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ABSTRACT 


The  external  nucleate  pool-boiling  heat-transfer 
coefficient  of  both  smooth  and  enhanced  horizontal  tubes 
in  R-114-oil  mixtures  (0  to  10  mass  percent  oil)  were 
measured  for  heat  fluxes  from  1  to  91  kW/m2  at  a  saturation 
temperature  of  2.2  C.  The  enhanced  tubes  tested  were  a 
Wieland  GEWA-T  finned  tube  containing  1.02  fins/mm,  a 
Hitachi  Thermoexce 1 -E  and  a  Hitachi  Thermoexce 1 -HE  tube. 
The  Thermoexce 1 -E  and  -HE  tubes  with  their  re-entrant 
cavity  designs  were  found  to  improve  the  heat- transfer 
coefficient  over  the  smooth  tube  value  at  a  constant  heat 
flux  by  a  factor  of  approx i mate  1 y 7  in  oil-free  R-114. 
while  the  GEWA-T  tube  improved  the  coefficient  by  a  factor 
of  about  4.  While  all  of  the  tubes  showed  a  generally 
decreasing  performance  with  the  presence  of  oil,  the  GEWA-T 
tube  resulted  in  the  minimum  reduction.  This  tube  showed 
up  to  20%  and  35%  reduction  (compared  to  the  oil-free  case) 
at  3  percent  and  10  percent  oil,  respectively.  The 
Thermoexce 1 -E  and  Thermoexce 1 -HE  tubes  showed  performance 
reductions  of  up  to  40%  and  60%,  respectively,  with  3 
percent  oil.  At  a  practical  heat  flux  of  about  30  kW/m: 
with  3  percent  oil,  the  GEWA-T.  Thermoexce 1 -E  and 
Thermoexce 1 -HE  tubes  outperformed  the  smooth  tube  by 
factors  of  4.8,  4.6  and  4.0,  respectively. 
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I.  INTRODUCTION 

A.  BACKGROUND 

In  the  design  of  refrigeration  systems  for  use  aboard 
Naval  vessels,  the  desire  is  to  make  the  system  as 
efficient  as  possible.  To  achieve  this  higher  efficiency, 
one  must  consider  the  use  of  a  high  pressure  refrigerant, 
since  this  gives  a  higher  energy  transfer  per  unit  volume 
of  vapor  (  i  .  e  .  .  Eu  =  h« , /v,  )  [Ref.  13,  and  results  in  less 
of  a  pressure  drop  along  the  vapor  path  for  a  given  heat 
load.  R-12,  which  is  a  high  pressure  refrigerant,  has  an 
energy  transfer  per  unit  volume  of  vapor  of  8570  kJ/m3 
[Ref.  23  and.  since  it  is  higher  pressure,  it  requires 
heavier  gauge  components.  R-ll.  on- the  other  hand,  is  a 
low  pressure  refrigerant  with  a  low  energy  transfer  per 
unit  volume  of  vapor  of  2086  kJ/m3  and  requires  lighter 
gauge  material . 

Table  1  provides  a  short  summary  of  some  of  the 
available  refrigerants  and  their  respective  energy  transfer 
per  unit  volume  of  vapor,  their  relative  toxicity  and  their 
required  system  pressure  for  optimum  operation.  As  can  be 
seen  from  Table  1,  R-12  and  R-22  possess  the  highest  values 
of  energy  transfer  per  unit  volume  of  vapor,  and  therefore 
would  provide  systems  of  higher  efficiency  than  the  other 
refrigerants  listed  since  pressure  losses  in  the  system 
would  be  less.  However,  since  both  R-12  and  R-22  are 


high-pressure  refrigerants,  their  respective  systems  would 
be  significantly  heavier  larger  than  those  of  the  other 
refrigerants.  However.  R-114  being  a  moderate-pressure 
refrigerant  with  an  energy  transfer  per  unit  vouiume  of 
vapor  of  3800  kJ/m3,  it  would  allow  for  a  relatively 
lighter  component  system  than  those  of  R-12  and  R-22. 


TABLE  1  [Ref.  33 
SUMMARY  OF  REFRIGERANTS 

Oper at i ona 1 
Sys  t em 

Re  f  r  i_ge r ant  k J  /m3 _ Toxicity _ _ _ Pressure 

R-ll  2086  moderate  (Group  5)  low 

R-12  8570  low  (Group  6)  high 

R-22  12480  moderate  (Group  5)  high 

R-113  1128  high  (Group  4)  low 

R-114  3800  moderate  (Group  6)  moderate 


Additionally,  R-114  belongs  to  the  refrigerant  group  (Group 
6)  having  the  lowest  toxioity  [ReF.  21.  and  it  has  also 
been  proven  fairly  stable  with  temperature  [Ref.  23. 
Therefore,  R-114  is  a  viable  alternative  to  the  higher 
pressure  refrigerants. 

B.  REVIEW  OF  REFRIGERANT-OIL  MIXTURE  BOILING  BEHAVIOR 

The  presence  of  oil  in  the  evaporators  of 
refrigeration  systems  is  a  common  occurrence.  This  is  due. 
in  part,  to  the  use  of  lubricants  in  the  hermetically 
sealed  refrigeration  compressors.  Over  a  period  of  time,  a 
considerable  amount  of  oil  is  introduced  into  the 
refrigerant.  subsequently  changing  the  physical 
properties  of  the  refrigerant.  Jensen  and  Jackman  [Ref.  33 


linearly  in  the  r e f r i ge r ant -o i 1  mixtures,  however,  the 


viscosity  and  the  surface  tension  do  not.  Henrici  and 
Hesse  [Ref.  43  experimentally  determined  in  1971  that  the 
surface  tension  of  mixtures  first  decreased  for  R-114-oil 
mixtures  of  up  to  2.5  percent  oil,  and  at  higher  oil 
concentrations,  the  surface  tension  increased.  This  non¬ 
linear  behavior  of  the  properties  of  the  mixtures  makes  it 
difficult  to  explain  the  changes  in  the  heat-transf er 
coefficient  for  these  mixtures. 

The  most  evident  example  of  the  addition  of  oil  to  the 
refrigerant  is  the  presence  of  foam  in  the  mixture  when 
boiling  takes  place.  With  the  oil  concentration  above  1 
percent,  there  is  a  substantial  amount  of  foam  generated 
from  the  nucleate  boiling  surface.  This  foam  is  generated 
because  the  refrigerant  in  the  ref r i gerant -o i 1  mixture  is 
more  volatile  than  the  oil  and  therefore  can  vaporize, 
creating  a  gas  bubble  surrounded  by  an  oil-rich  layer 
[Ref.  53.  Figure  1.1  shows  the  ideal  model  of  the  growth  of 
a  bubble  in  a  refrigerant-oil  mixture.  These  generated 
bubbles  are  at  lower  density  than  that  of  the  surrounding 
re f r i gerant -o i 1  mixture.  The  bubbles  ascend  to  the  liquid- 
vapor  interface  and  collect  at  the  interface  and  produce  a 
foam  layer.  As  reported  by  Henrici  and  Hesse  [Ref.  4],  the 
foaming  action  is  most  pronounced  in  refrigerant-oil 
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aft 
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mixtures  of 


1  to  10  percent  oil,  and  provide  a  the 


presence  of  oil  results  in  a  subsequent  decrease  in  the 
heat-transfer  performance  of  the  boiling  tube. 

Numerous  researchers  have  tried  to  explain  the  reasons 
for  the  decrease  in  the  heat- transfer  coefficient  in 
r ef r i gerant-o i 1  mixtures  of  increasing  oil  content.  Thome 
[Ref.  6],  after  an  extensive  literature  search,  found  that, 
in  1956,  Van  Wijk  et  al.  explained  the  effect  as  being  the 
result  of  the  evaporation  of  the  more- vo 1  at i 1 e  component 
leaving  an  oil-rich  layer  which  has  a  higher  resultant 
local  boiling  point.  This  increased  local  boiling  point 
requires  an  increase  in  the  amount  of  superheat  necessary 
to  continue  the  vaporization  and  bubble  growth.  The  result 
is  a  decrease  in  the  heat-transf er  coefficient. 

Stephan  and  Preusser  [Ref.  7]  determined  that  the  work 
required  to  form  bubbles  in  the  ref rigerant-oi 1  mixtures 
was  greater  than  in  an  equivalent  amount  of  pure 
refrigerant.  Therefore,  the  heat-transf er  coefficient  was 
lower  for  the  refrigerant-oil  mixture  than  for  the  pure 
refrigerant.  Chongrungr eong  and  Sauer  [Ref.  8]  reported 
that  the  rate  of  heat  diffusion,  which  is  governed  by  the 
thermal  properties  of  the  oil-rich  layer,  limits  the  bubble 
growth  and  that  the  effects  of  surface  tension  are 
negligible.  Thome  [Ref.  6]  agreed  with  the  above  factors 


and  concluded  that 


the  viscosity  variations  are  important 


in  explaining  the  decrease  in  the  heat-transf er 
coefficient  for  some  ref r i ger ant - o i 1  mixtures. 

Although  researchers  have  not  decided  on  any  one  set 
of  factors  that  determine  the  reasons  for  the  decrease  in 
the  heat-transf er  coefficient  as  oil  content  increases  in 
refrigerant-oil  mixtures,  all  agreed  that  the  physical  and 
thermal  properties  are  important  factors  in  explaining  the 
heat- transfer  behavior  of  r e f r i gerant -o i 1  mixtures. 

C.  REVIEW  OF  BOILING  PERFORMANCE  FROM  SMOOTH  TUBES 

The  U.  S.  Navy  attempts  to  reduce  the  size  of  its 
shipboard  refrigeration  systems  by  the  use  of  alternate 
refrigerants  (such  as  R-114).  At  the  same  time,  the  Navy 
desires  to  increase  the  performance  of  these  systems  by  the 
use  of  enhanced  surfaces.  However,  as  discussed  by 
Wanniarachchi  et  al.  [Ref.  9],  only  a  limited  amount  of 
data  is  available  regarding  the  boiling  of  refrigerants  in 
the  presence  of  oil. 

Stephan  [Ref.  10]  conducted  experiments  using  a 
smooth  horizontal  plate  as  a  boiling  surface  in  R-12 
and  R-22  refrigerants  and  reported  up  to  50  percent 
reductions  in  the  boiling  hea t - t r ansf er  coefficient 
with  a  refrigerant-oil  mixture  of  9  percent  oil  (by 
mass),  while  50  percent  oil  resulted  in  a  90 
percent  reduction  in  the  boiling  heat-transf er  coefficient. 


In  1971,  Henrici  and  Hesse  [Ref.  43  conducted 
experiments  involving  the  boiling  of  R-114  in  the  presence 
of  oil  from  a  smooth  copper  tube  and  reported  a  maximum 
decrease  of  20  percent  decrease  in  the  boiling  heat- 
transfer  coefficient  with  1  percent  oil  by  mass.  In  1981, 
Stephan  and  Mitrovic  CRef.  113  reported  that  for  the  GEWA-T 
tube  (manufactured  by  the  Wieland  Company)  in  R-114-oil 
mixtures,  the  boiling  heat-transfer  coefficient  of  the 
enhanced  surface  was  significantly  altered  by  the  presence 
of  the  oil.  In  1985,  Reilly  [Ref.  53  reported  a  drop  in 
the  heat- transfer  coefficient  of  0  to  35  percent  (depending 
on  the  heat  flux)  for  a  R-114-oil  mixture  of  10  percent  oil 
for  a  smooth  copper  tube. 

D.  REVIEW  OF  BOILING  PERFORMANCE  FROM  ENHANCED  SURFACES 

Aral  et  al.  [Ref.  123  tested  a  200-ton  R-12 
centrifugal  water  chiller.  The  chiller  which  used  an 
enhanced-surf ace  tube,  commercially  known  as  Thermoexce 1 -Em 
made  by  Hitachi,  was  28  percent  shorter  in  length  and  had 
an  approximate  improvement  in  the  overall  heat-transf er 
coefficient  of  50  percent.  Additional  research  by  Yilmaz 
and  Uestwater  [Ref.  133,  Marto  and  Lepere  [Ref.  143,  and 
Carnavos  [Ref.  153  on  commercially  available  enhanced 
surface  tubes,  in  various  refrigerants  other  than  R-114, 
indicated  that  a  porous -coated  surface  exhibited  the  best 
boiling  heat- transfer  coefficient  in  a  pure  refrigerant. 


Comparisons  by  Reilly  [Ref.  53  indicated  that  in  R-114, 
Union  Carbide’s  "High  Flux"  tube  with  a  porous-coated 
surface  resulted  in  an  improvement  of  at  least  a  factor  of 
seven  over  a  smooth  copper  tube  at  the  same  operational 
parameters.  Reilly’s  findings  complemented  the  research  of 
Yilmaz  and  Westwater  [Ref.  133.  Additionally,  Reilly 
reported  that  for  a  porous-coated  surface  in  an  R-114-oil 
mixtures  of  1  percent  oil,  the  oil  was  seen  to  degrade  the 
nucleate  pool-boiling  heat-transf er  performance  by  about 
20  percent. 

1.  Surfaces  with  Reentrant  Cavities 

Reentrant  cavities  have  been  investigated  for  a 
considerable  period  of  time.  Griffith  and  Wallis  [Ref.  153 
determined  that  the  reentrant  cavity  geometry  is  important 
in  the  generation  of  nucleate  boiling  sites  in  two  ways; 
(1)  the  diameter  of  the  reentrant  cavity  mouth  determines 
the  amount  of  superheat  required  to  initiate  boiling, 
and  (2)  the  cavity  shape  determines  the  stability  once 
the  boiling  has  begun.  Webb  [Ref.  163  stated  that  the 
key  to  the  high  performance  of  the  reentrant  grooved 
structures  can  be  attributed  to  three  factors:  (1) 
reentrant  cavity  within  a  critical  size  range,  (2) 
interconnected  cavities,  and  (3)  nucleation  sites  of  a  re¬ 
entrant  shape.  If  the  cavities  are  interconnected, 
adjacent  cavities  can  activate  each  other.  Webb  [Ref.  163 


also  stated  that  reentrant  cavities  provide  a  stable  vapor 
trap,  which  can  remain  active  at  low  values  of  superheat. 

The  GEUA-T  (1.02  fins/mm)  enhanced-sur f ace  tube, 
manufactured  by  Wi e 1 and-Wer ke  AG,  is  a  copper  tube  with 
spiral  fins  of  a  T-shaped  profile.  It  has  a  reentrant- type 
cavity  formed  by  cold  working  an  integral  finned  tube.  The 
two  adjacent  T-fins  form  a  spiral  cavity  with  an  opening 
of  0.25-mm  (see  Figure  1.2). 

The  Thermoexce 1 -E  and  Thermoexce 1 -HE  enhanced-sur face 
tubes  are  made  by  the  Hitachi  Company.  Like  the  GEWA-T 
tube,  the  Thermoexce 1 -E  (Figure  1.3)  and  Thermoexce 1 -HE 
(Figure  1.4)  tubes  are  also  formed  by  cold  working  integral 
fins?  however,  these  tubes  have  low  fins  which  have  a 
small  space  cut  out  at  the  fin  tips,  giving  the  appearance 
of  a  sawtooth.  These  sawtooth  fins  are  bent  parallel  to 
the  tube  axis,  to  a  horizontal  position,  forming  tunnels 
with  evenly  spaced  surface  pores.  This  cold  working 
procedure  results  in  a  high  area  density  of  reentrant 
nucleation  sites.  The  Thermoexce 1 -HE  tube  is  a  variation 
of  the  Thermoexce 1 -E  tube,  designed  for  improved 
performance  in  1 ow- hea t - f 1 ux  regions. 

Experimental  data  showing  the  effect  of  oil 
concentrations  on  the  hea t - t r ans f er  coefficient  of  the 
GEWA-T,  the  Thermoexce 1 -E  and  Thermoexce 1 -HE  tubes  in  R-114 


are  lacking. 


In  an  effort  to  gain  more  understanding  of 


the  pool-boiling  performance  of  these  tubes,  the  David  W. 


Taylor  Naval  Ship  Research  and  Development  Center  sponsored 
this  investigation.  The  smooth  copper  tube  and  the  GEWA-T 

tube  were  provided  by  the  Wieland  Company.  The  Thermoexce I  - 
E  and  Thermoexce 1 -HE  tubes  were  provided  by  the  Hitachi 
Company . 

E.  THESIS  OBJECTIVE 

The  objectives  of  this  thesis  are: 

1.  take  data  on  a  smooth  tube  with  the  use  of 

auxiliary  heaters  together  with  a  convection  shield 
to  enable  the  collection  of  more  data  with  a  minimum 
of  operator  effort, 

2.  take  boiling  data  on  a  smooth  copper  tube  in  R-114 

to  provide  a  baseline  for  follow-on  comparison 
for  enhanced  tubes  at  oil  concentrations  of  0,  1. 

2,  3,  6  and  10  percent  (by  mass)  at  a  boiling 

temperature  of  2.2  C  and  a  heat-flux  range  from  500 

W/m2  to  91  kW/m2 . 

3.  take  data  on  the  following  tubes  for  the  above- 
-mentioned  conditions: 

a.  GEWA-T  (1.02  fins/mm), 

b.  Thermoexce 1 -E,  and 


c.  Thermoexce 1 -HE . 


II.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 


A.  OVERALL  APPARATUS 

The  equipment  used  for  this  study  consisted  of  nine 
basic  components:  (a)  a  Pyrex-glass  tee  for  the  pool 
boiling  of  the  R-114  liquid;  (b)  a  Pyrex-glass  tee  for  the 
condensing  of  the  R-114  vapor;  (c)  an  R-114  liquid 
reservoir;  (d)  a  water-ethylene  glycol  mixture  sump;  (e)  an 
R-12  refrigeration  system  with  an  external  proportioning 
valve;  (f)  oil  reservoir;  (g)  a  vacuum  pump; 
(h)  instrumentation  and  data-acqu i s i t i on  system!  and  (i)  a 
convection  shield.  Figure  2.1  shows  schematically  the 
arrangement  of  the  various  components.  Complete  details  of 
the  design,  construction  and  operation  of  the  apparatus 
are  provided  by  Karasabun  [Ref.  17],  Reilly  [Ref.  5] 
and  Wanniarchchi  et  al.  [Ref.  13.  The  original  design  of 
Karasabun  was  modified  by  the  inclusion  of  the  convection 
shield  in  the  pool-boiling  section  as  shown  in  Figures  2.2 
and  2.3. 

The  operation  of  the  apparatus  consisted  of  the 
boiling  of  R-114  liquid  in  the  lower  glass  tee  (1), 
and  condensing  of  the  R-114  vapor  in  the  upper 
glass  tee  (2).  The  resultant  condensate  was  fed  to 
the  distribution  tube  within  the  boiling  section  by 
gravity.  A  water-ethylene  glycol  mixture  from  a  30-gallon 


vv 


sump  was  used  to  cool  the  copper  condensing  coil,  thereby 


condensing  the  R-114  vapor.  This  water -ethy 1 ene  glycol 
mixture  sump  was  maintained  between  -18  and  -14  C  by  a 
one-half  ton  R-12  air-conditioning  system.  The  mixture  was 
pumped  through  the  condensing  coil  and  controlled  manually 
with  valve  VC. 

Valve  VC  controls  the  amount  of  the  water-ethylene 
glycol  mixture  to  the  condensing  coil,  thereby  controlling 
the  pressure  within  the  test  loop.  For  example,  the  opening 
of  valve  VC  allows  more  of  the  water -ethy 1 ene  glycol 
mixture  to  pass  through  the  condensing  coil,  resulting  in 
additional  condensing  of  the  R-114,  and  subsequently 
lowering  the  temperature  and  pressure  within  the  entire 
system. 

As  the  subcooled  R-114  condensate  enters  the  boiling 
section  of  the  apparatus  through  the  distribution  tube,  it 
strikes  the  convection  shield.  This  shield  provides 
increased  mixing  of  the  subcooled  R-114  condensate  with  the 
remaining  R-114  (or  R-114/oil  mixture)  in  the  lower  portion 
of  the  boiling  section.  The  shield  deflects  the  incoming 
condensate  away  from  the  boiling  surface  and  allows  for 
further  mixing  of  the  condensate  and  the  existing  R-114 
liquid  by  means  of  the  auxiliary  cartridge  heaters.  This 
deflection  of  the  incoming  condensate  allowed  for  a  more 
even  temperature  distribution  of  the  liquid  R-114  and 
therefore  decreased  the  previously  reported  temperature 
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gradients  within  the  liquid  region.  This  additional  mixing 
of  the  liquid  R-114,  coupled  with  the  ability  to  use 
auxiliary  cartridge  heaters  in  the  liquid  and  allowed  for 
the  timely  achievement  of  steady-state  conditions  without 
the  previously  observed  thermal  gradients  in  the  liquid 
region  of  the  tests  done  without  the  aid  of  the  convection 
shield.  It  was  also  observed  throughout  the  operation  of 
the  apparatus  that  the  data  obtained  using  the  convection 
shield  had  a  higher  degree  of  repeatability  than  the  data 
obtained  without  the  shield. 

B.  BOILING-TUBE  CONSTRUCTION 

The  tubes  tested  during  the  course  of  this  study  were 
the  smooth  copper  tube,  the  GEWA-T  (1.02  fins/mm),  the 
Thermoexce 1 -E  and  the  Thermoexce 1 -HE  tubes.  All  of  these 
tubes  had  an  active  test  section  (i.e.,  heated)  of  length 
of  203.2-mm  (8  in),  with  an  outside  diameter  of  15.9  mm 
(5/8  in)  and  an  inside  diameter  of  12.7  mm  (1/2  in).  The 
remaining  114.3-mm  section  on  either  side  of  the  heated 
test  section  was  smooth  and  unheated  and  did  not  nucleate 
under  any  heat  flux  or  any  oil  concentration  condition. 
Figure  1.2  shows  the  schematic  of  the  GEWA-T  finned 
surface.  Karasabun  [Ref.  17]  described  the  analysis  used 
in  the  da ta - reduct i on  program,  allowing  for  the  two  end 
surfaces  to  be  treated  as  extended  fins  from  the  center 
section  and  subsequently  accounting  for  the  heat  losses. 


The  center  section  of  each  tube  tested  was  heated  by  a 
1000-W  240 -V  stainless  steel  cartridge  h6ater  inserted  in 
the  center  section  of  each  tube.  The  cartridge  heater  was 
inserted  in  the  boiling  tube  and  surrounded  by  a  copper 
sleeve  with  eight  1.3  mm  by  1.3  mm  thermocouple  channels. 
Figure  2.4  shows  the  layout  of  the  installed  tube  in  the 
boiling  section.  Figure  2.5  shows  the  details  of  the 
thermocouple  layout.  The  thermocouple  hot  junctions  were 
inserted  into  the  copper  sleeve  and  the  copper  sleeve  was 
then  tinned;  and  heat  was  maintained  to  keep  the  solder  in 
a  molten  state  as  the  sleeve  was  inserted  and  positioned  in 
the  center  section  of  the  boiling  tube.  Upon  completion  of 
the  thermocouple  installation,  calibration  of  the 
thermocouples  was  achieved.  Reilly  [Ref.  53  determined  in 
the  course  of  his  investigation  that  since  the  data- 
reduction  program  utilized  the  differences  between  the 
thermocouples  in  all  computations,  such  as  wall 
temperatures  minus  saturation  temperatures,  the  corrections 
provided  by  the  thermocouple  calibration  were  basically 
unnecessary  since  the  calibration  procedures  were  only 
necessary  for  items  dependent  on  absolute  temperature. 
Appendix  A  summarizes  the  procedure  used  by  Karasabun 
[Ref.  17]  and  Reilly  [Ref.  5]  to  calibrate  the 
thermocoup  1 es. 
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Figure  2.4  Sectional  View  of  Boiling  Tube 


Bottom 


C.  CONVECTION  SHIELD 

Initial  tests  on  the  smooth  copper  tube  revealed  an 
axial  thermal  gradient  in  the  liquid  over  the  length  of  the 
boiling  tube.  This  thermal  gradient  increased  with 
increasing  oil  concentration  and  with  increasing  heat  flux, 
with  the  thermal  gradient  ranging  from  0.5  K  to  1.5  K  over 
the  length  of  the  boiling  section.  This  thermal  gradient 
resulted  in  questionable  data  since  steady-state  conditions 
could  not  be  maintained  over  the  course  of  the  individual 
data  runs. 

During  the  early  stages  of  this  investigation,  it  was 
found  that  the  steady-state  conditions  were  fairly 
difficult  to  achieve.  This  was  especially  true  after  making 
a  step  change  in  the  heat  flux.  For  example,  if  the  heat 
flux  was  increased,  the  cooling  flow  rate  to  the  condenser 
must  be  increased  to  maintain  the  system  pressure.  Since 
the  system  showed  about  a  30-second  time  lag  between  the 
change  of  the  valve  position  and  the  subsequent  system 
response,  the  pressure  was  often  seen  to  undergo  cyclic 
variations.  Uith  considerable  experience  of  the  operator, 
it  was  possible  to  achieve  steady-state  conditions  within  5 
to  10  minutes,  as  demonstrated  by  Reilly  [Ref.  5]. 

In  order  to  minimize  the  efforts  necessary  by  the 
operator,  two  600-W  auxiliary  cartridge  heaters  were 
installed  within  the  boiling  tee.  With  this  arrangement. 


steady-state  conditions  could  be  maintained  very  easily  by 


maintaining  the  total  heat  load  in  the  boiling  tee  at  a 
constant  value.  For  example,  if  the  heat  duty  for  the 
boiling  tube  was  increased,  the  heat  duty  on  the  auxiliary 
heaters  must  be  decreased  by  the  same  magnitude.  In  this 
manner,  steady-state  conditions  could  be  maintained  with 
minimum  attention.  However,  this  arrangement  required  that 
the  heat  duty  on  the  auxiliary  heaters  be  a  maximum  when 
the  heat  duty  on  the  boiling  tube  is  at  a  minimum.  As 
determined  by  the  early  data  runs  (see  Figure  2.6),  the 
convective  effects  created  by  the  bubbles  of  the  auxiliary 
heaters  created  an  artificially  enhanced  hea t - t r ans f er 
coefficient  at  low  heat-flux  settings. 

In  order  to  alleviate  the  problem  discussed 
above,  while  maintaining  the  advantages  offered  by 
the  auxiliary  heaters,  the  present  investigator 
considered  the  use  of  a  convection  shield.  As 
shown  in  Figure  2.2,  this  shield  was  positioned 
around  the  test  tube,  thus  isolating  it  from  the 
convective  bubble  effects  created  by  the  auxiliary  heaters. 

The  convection  shield  also  provided  increased 
mixing  of  the  refrigerant  and  condensate  without 
interfering  with  the  boiling  phenomenon  of  the  boiling 
tube  or  the  circulation  of  the  o i 1  - r ef r i gerant 
mixture  in  the  boiling  section  by  both  the  boiling 
tube  and  the  auxiliary  heaters. 


Figure  2.6  Convective  effects  of  Auxiliary  Heaters 
with  and  without  Convective  Shield 


D. 


DATA  ACQUISITION  AND  REDUCTION 


The  data-acqu i s i t i on  system  used  for  this  research  was 
a  Hewlett-Packard  9826A  computer  and  a  Hewlett-Packard 
3497A  Automatic  Data  Acquisition  Unit.  The  HP-9826A 
computer  was  used  to  control  the  HP-3497A,  to  provide  hard- 


*  copy 

i 


storage  capacity 


and  to  analyze  the  data  obtained 


during  this  investigation. 

The  data-acqu i s i t i on  system  provided  the  means  by 
which  the  researcher  was  able  to  review  system  performance 
and  stability  prior  to  the  taking  of  any  data.  System 
parameters  were  monitored  on  the  video  screen  following  the 
scanning  of  all  of  the  channels. by  the  HP-3497A.  Following 
the  automatic  scanning  of  all  of  the  channels,  the  computer 
would  display  the  desired  saturation  temperature  (2.2  °C). 
the  observed  saturation  temperature,  the  liquid  temperature 
at  both  ends  of  the  boiling  section  and  their  average,  the 
vapor  temperature,  and  the  temperature  of  the  water- 
ethylene  glycol  mixture  sump  (-15  to  -18  °C).  Once  the 
steady-state  conditions  were  met  for  the  required  satura¬ 
tion  temperature  within  the  prescribed  variance  (iO.l  K) 
of  the  desired  saturation  temperature,  the  data-acqui s i t i on 
system  was  initiated  for  the  individual  data  run. 
Since  there  was  a  time  lag  in  the  request  for  the 
taking  of  data  and  the  actual  printing  of  the  data  due  to 
the  sampling  of  each  channel  twenty  times  before  recording 


the  data,  the  da ta - acqu i s i t i on  system  would  print  the  data. 


al lowing  the  researcher  to  review  the  results  prior  to 
storing  the  data.  If,  during  the  sampling  of  the  data,  a 
transient  trend  was  observed  due  to  the  cycling  of  the 
R-12  refrigerant  system  or  the  temperature  difference 
between  the  thermocouples  measuring  the  liquid  R-114  was 
not  within  the  specified  temperature  range,  the  researcher 
had  the  capability  to  reject  the  affected  data  set.  This 
capability  allowed  the  researcher  to  apply  a  cursory 
analysis  of  each  data  point  and  ascertain  if  the 
temperatures  measured  during  that  individual  data  set  were 
within  the  prescribed  limits  of  variance. 

The  information  scanned  on  each  channel  of  the  HP- 
3497A  were  thermal  emf’s  provided  by  copper-constantan 
(type-T)  thermocouples  made  of  0 . 25 -mra-d i ame tc r  (30  gauge) 
wire.  Two  additional  measurements  were  provided  to  the  HP- 
3497A  by  the  power  -  sens i ng  device,  as  described  by 
Karasabun  [Ref.  17],  which  converted  the  AC  voltage  dialed 
into  the  variac  of  the  main  boiling  tube  heater  into  a  DC 
signal,  allowing  the  HP-3497A  to  accurately  scan  the 
converted  value  and  then  allow  the  HP-9826A  to 
automatically  apply  this  value  to  the  analysis  calculation. 

Upon  the  initiation  of  the  command  to  take  data,  the 
data-acqu i s i t i on  system  would  sample  each  channel  and  then 
compute,  according  to  the  step-wise  procedure  outlined  by 
Karasabun  [Ref.  17]  and  print  the  results.  Appendix  B 
contains  a  summary  of  the  data-acquisition  program  and  a 


listing  of  the  da ta - r educ t i on  program  utilized  during  the 
course  of  this  research.  Appendix  C  contains  a 
representative  sample  of  the  data  printout  received  for 


each  data  run. 


I  I  I_._  EXPER  I  MENTAL  P R OCEDURES 

A.  INSTALLATION  OF  TUBE  IN  BOILING  APPARATUS 

The  external  surfaces  of  each  of  the  boiling  tubes 
were  cleaned  with  Nitol  (2  percent  nitric  acid  and  98 
percent  ethyl  alcohol)  solution  and  then  rinsed  with 
acetone  to  ensure  tube  cleanliness.  This  treatment  was  also 
effective  in  the  removal  of  oil  and  grease  residues  from 
the  tube  surfaces. 

After  air  drying,  the  individual  tube  to  be  tested  was 
installed  in  the  boiling  section  of  the  apparatus.  The 
system  was  sealed  and  then  evacuated  to  approximately  29 
inHg  by  means  of  .  a  portable  mechanical  vacuum  pump  (6). 
shown  in  Figure  2.1.  The  apparatus  remained  under  the 
evacuated  condition  for  a  period  of  about  one  hour  to  check 
if  there  were  any  major  leaks.  System  pressure  was  measured 
by  means  of  a  Marsh  pressure  gauge  (30  inHg  to  150  psi 
range  with  12.5  inHg  and  £0.5  psi  accuracy, 
respectively).  If  there  was  not  a  noticeable  drop  in  the 
vacuum  reading  as  observed  on  the  pressure  gauge,  the 
system  was  charged  with  R-114  vapor,  which  was  condensed  by 
the  wa t e r - e t hy I ene  glycol  mixture  and  then  the  absolute 
pressure  was  raised  to  approximately  0.186  MPa  (27  psi). 


Upon  reaching  this  pressure. 


an  Automatic  Halogen  Leak 


Detector,  TIF  5000,  with  a  sensitivity  of  3  ppm 
concentration,  was  used  to  check  for  system  leaks. 

Upon  successful  completion  of  the  leak-detection  tests 
and  taking  the  necessary  steps  to  fix  any  leaks  found,  the 
R-114  liquid  level  was  adjusted  to  be  20  mm  above  the  top 
of  the  boiling  tube  surface.  Liquid-level  adjustments  were 
made  by  either  boiling  off  the  excess  refrigerant  to  the 
R-114  reservoir  (3),  or  drawing  the  required  amount  of 
refrigerant  liquid  into  the  system  from  the  R-114 
reservoir.  The  required  liquid  level  for  boiling  was 
maintained  constant  for  al 1  tubes  tested  by  means  of  a 
permanent  scribe  mark  on  the  inside  of  the  oil  intake 
endbe 1 1  of  the  boiling  section.  At  this  initial  refrigerant 
level,  the  system  was  now  ready  for  boiling  tests  in  pure 
refrigerant. 

B.  GENERAL  SYSTEM  OPERATION 

Table  2  provides  a  summary  of  all  of  the  83  data  runs 
made  during  this  investigation.  Each  data  run  is  annotated 
to  indicate  the  conditions  involved  in  the  respective  run. 
The  data  runs  were  numbered  sequentially  with  an  indicator 
code  for  each  tube  preceding  the  number  of  the  data  run. 
Since  this  research  is  a  follow  up  of  previous  research, 
the  initial  data  run  was  indicated  as  number  168.  All  data 
runs  were  undertaken  at  a  saturation  temperature  of  2.2°C. 


A 

listing  of 

the 

boiling  tube  indicator  codes 

for  Table  2 

i  s 

as  foil ows : 

WH 

- -  W i e 1  and 

Hard  Copper  Smooth 

Tube 

GTB 

- -  W i e 1  and 

Gewa-T  (1.02  fins 

per  mm) 

TXE 

--  Hitachi 

Thermoexce 1 -E 

THE 

- -  Hi tach i 

Thermoexcel-HE. 

Initially,  all  data  runs  were  taken  over  a  range  of 
nine  different  heat-flux  settings,  ranging  from  2  kW/m2  to 
91  kW/m2 .  Additional  heat  fluxes  of  0.5,  0.8,  1.2  kW/m2 
were  later  included  for  al 1  tubes  with  enhanced  surfaces  to 
study  the  onset  of  nucleate  boiling  in  the  o i 1  -  ref r i gerant 
mixtures  and  to  study  the  effects  of  hysteresis  on  these 
tubes.  The  initiation  of  all  data  runs  was  carried  out  in 
the  same  manner  to  ensure  consistency  in  the  obtained  data. 
Each  run  was  started  by  an  initial  system  "warm-up"  period 
allowing  the  R-12  refrigeration  system  to  obtain  a  minimum 
coolant  starting  temperature  of  about  -15°C.  This  was 
accomplished  by  lighting  off  of  the  R-12  system  in 
conjunction  with  the  effluent  recirculation/cooling  pump 
(8)  and  slight  opening  of  the  control  valve  VC.  This 
procedure  allowed  the  system  to  gradually  decrease  system 
temperature  and  pressure,  while  minimizing  the  adverse 
effects  of  shocking  the  boiling  tube.  The  system  was 
operated  in  this  condition  until  a  stable  condition  was 
established  with  the  saturation  temperature  of  2.2°C.  This 
operation  was  most  critical  when  the  testing  sequence 
involved  increasing  heat-flux  conditions  since 
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TABLE  2 

Summary  of  Data  Runs 


Run 

No . 

% 

Oi  1 

Heat  No . of  Data 
Flux  Points 

Aux i 1 iary 
Heater 

Convect i on 
Shield 

WH168 

0 

DEC 

20 

off 

wi thout 

WH 169 

0 

INC 

22 

off 

w i thout 

WH170 

1 

DEC 

20 

off 

w i thout 

WH171 

1 

INC 

22 

off 

wi thout 

WH172 

2 

DEC 

20 

off 

wl thout 

WH173 

2 

INC 

20 

off 

w i thout 

WH174 

3 

DEC 

20 

off 

wi thout 

WH175 

3 

I  NC 

22 

off 

wi thout 

UH176 

3 

DEC 

20 

off 

wi thout 

WH  1 77 

6 

DEC 

20 

off 

without 

WH178 

6 

INC 

21 

off 

w i thout 

WH  1 79 

10 

DEC 

18 

off 

w i thout 

WH180 

10 

INC 

13 

off 

w i thout 

WH  181 

10 

INC 

12 

off 

w i thout 

GTB182 

0 

DEC 

24 

on 

with 

GTB183 

0 

INC 

23 

on 

with 

GTB184 

0 

DEC 

22 

on 

with 

GTB185 

1 

DEC 

20 

on 

with 

GTB186 

1 

INC 

22 

on 

with 

GTB187 

2 

DEC 

34 

on 

with 

GTB188 

2 

INC 

18 

on 

with 

GTB189 

3 

INC 

22 

on 

with 

GTB190 

3 

DEC 

20 

on 

with 

GTB191 

6 

DEC 

27 

on 

with 

GTB192 

6 

I  NC 

20 

on 

with 

GTB193 

10 

INC 

20 

on 

with 

GTB194 

10 

DEC 

21 

on 

with 

WH  195 

0 

I  NC 

22 

on 

with 

WH196 

0 

DEC 

20 

on 

with 

WH  197 

0 

DEC 

17 

on 

with 

WH  1 98 

0 

INC 

22 

on 

with 

WH199 

1 

INC 

20 

on 

with 

WH200 

1 

DEC 

20 

on 

with 

WH201 

2 

INC 

20 

on 

with 

WH202 

2 

DEC 

20 

on 

with 

WH203 

3 

I  NC 

20 

on 

with 

WH204 

3 

INC 

20 

on 

with 

WH205 

3 

DEC 

20 

on 

with 

WH206 

3 

INC 

21 

on 

with 

WH207 

6 

DEC 

20 

on 

with 

WH208 

6 

INC 

20 

on 

with 

WH209 

10 

INC 

20 

on 

with 
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WH210 

TXE210 

TXE211 

TXE212 

TXE213 

TXE21A 

TXE215 

TXE216 

TXE217 

TXE218 

TXE219 

TXE220 

TXE221 

TXE222 

TXE223 

TXE224 

TXE225 

TXE226 

TXE227 

TXE228 

TXE229 

TXE230 

THE231 

THE232 

THE233 

THE234 

THE235 

THE236 

THE237 

THE238 

THE240 

THE241 

THE242 

THE243 

THE244 

THE245 

THE246 

THE247 

THE248 

THE249 


TABLE  2 

Summary  of 

Data  Runs 

( cont ' d ) 

*  Heat 

No . of  Data 

Aux i 1 i ary 

Convect i on 

Oil  Flux 

Points 

Heater 

Shield 

f  a  i  lure 


to  strictly  adhere  to  this  would  result  in 
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premature  nucleate  boiling  and  the  subsequent  purging  of 
the  data.  This  initial  warm-up  period  required  a  minimum 
of  one  hour  to  ensure  adequate  system  performance. 

Upon  reaching  stable  conditions,  the  initial  heat  flux 
for  the  boiling  tube  was  set  by  means  of  the  resistor 
variac.  Also,  the  heat  flux  for  the  auxiliary  cartridge 
heaters  was  set  by  a  second  variac.  The  heat  flux  of  the 
auxiliary  cartridge  heaters  was  limited  to  settings  that 
allowed  the  heating  of  the  liquid  R-114  without  creating 
excess  turbulence  on  the  boiling  tube.  As  discussed 
earlier,  the  boiling  tube  was  shielded  from  the  effects  of 
the  auxiliary  heater  by  means  of  the  convection  shield.  As 
the  oil  concentration  of  the  refrigerant-oil  mixture 
increased,  the  heat-flux  setting  for  the  auxiliary  heater 
was  additionally  limited  to  lower  heat  fluxes.  This 
precluded  the  generated  oil  foam  from  carrying  over  the 
sides  of  the  convection  shield  and  interfering  with  the 
operation  of  the  boiling  tube. 

After  setting  of  the  desired  heat  flux  for  the  data 
run,  the  system  was  monitored  for  the  vapor  temperature, 
sump  temperature  and  the  liquid  refrigerant  temperature 
(two  readings;  left  side  and  right  side  of  the  boiling 
tube).  Control  valve  VC  was  adjusted  to  allow  the  system  to 
reach  steady-state  conditions.  Additionally,  the  liquid 
refrigerant  temperatures  were  compared  to  ensure  minimum 


TV* 


m 


Swra 


variance.  After  fulfilling  the  requirements  of  the  steady- 
state  conditions,  the  data-acqui si t ion  system  was  initiated 
to  take  data. 

The  data-acqu i s i t i on  system  would  scan  each  assigned 
channel,  compute  the  heat-transf er  coefficient  and  print 
the  results  for  operator  viewing.  If  the  data  were  within 
the  prescribed  limits,  they  were  stored  on  a  floppy  disk 
for  permanent  records.  A  representative  sample  of  a  data 
printout  is  contained  in  Appendix  B. 

After  each  run  was  completed,  oil  was  added  to 
the  boiling  section  by  means  of  a  graduated  cylinder 
through  valve  VI.  Table  3  contains  the  listing  of 


the 

requi red 

vo 1 ume 

of 

oi  1 

added  to 

obtain 

the 

proper 

mass 

percent 

of 

oi  1 

for  each  series  of 

data 

runs . 

Mixing 

of  the 

o  i  1 

and 

refrigerant 

1 i qu i d 

was 

accomp 1 i shed 

by  means 

of 

the 

aux i 1 i ary 

heater , 

to 

ensure  a  homogeneous  mixture.  Foaming  of  the  oil- 
refrigerant  mixtures  was  observed  over  the  entire  heat- 
flux  range  of  each  data  run,  although  carryover  of 
the  foam  over  the  convection  shield  was  limited  to  the 
highest  heat  flux.  As  the  oil  concentration  increased, 
carryover  increased,  but  was  limited  to  the  vertical 
section  of  the  apparatus.  Carryover  oil  was  not 
observed  in  the  condenser  section  of  the  apparatus. 
Figures  3.1  to  3.3  demonstrate  the  foaming  action  of 


the 


o i 1  -  ref r i gerant  mixtures 


at 


a 


heat  flux  of  37 


>A 

PS 


kW/m*  for  oil  concentrations  of  0,  6  and  10  percent. 

Data  runs  were  made  for  each  boiling  tube  tested  with 
pure  R-114  and  followed  by  runs  of  oil-refrigerant 
mixtures  of  increasing  oil  content.  0 i 1  - r e f r i gerant 

mixtures  of  1,  2,  3,  6  and  10  percent  oil  by  mass  were 

used  for  al 1  tubes. 

TABLE  3 

OIL  CONCENTRATIONS  VALUES 


Percent 
Oi  1 


0 

1 

2 

3 

6 

10 


Total  Volume  (cc)  Step  Change  Ccc) 
of  Oil _ in  Vo  1  ume  of  Oil 


26. 

8 

26. 

8 

54. 

2 

27. 

4 

82. 

2 

28. 

,  0 

169. 

,  6 

87. 

.  4 

295. 

,  2 

125. 

,  6 

M 


37  kW/m 


91  kW/m 


Figure  3.2 


Foaming  Action  for  R-114  with 
6  Percent  Oil 
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91  kW/m‘ 


Figure  3.3 


Foaming  Action  for  R-114  with 
10  Percent  Oil 
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Figure  4.1  Performance  Comparison  for  Pure  R-114  Boiling 
from  a  Smooth  Tube. 


Smooth  Tube. 


Figure  4.3  Performance  Comparison  for  R-114  with  10  Percent 
Oil  Boiling  from  a  Smooth  Tube. 


-2.2°C.  This  data  set  represents  a  series  of  data 
runs,  starting  with  a  heat  flux  of  100  kW/m2  and 
decreasing  the  heat  flux  to  a  minimum  of  1  kW/m2 .  The 
three  sets  of  curves  are  nearly  identical  with  the 
deviation  between  the  Reilly  data  and  the  most  recent  data 
being  attributed  to  the  pressure  effect  and  the  possible 
aging  of  the  tinned  interface  between  the  thermocouple 
sleeve  and  the  inner  surface  of  the  smooth  copper  tube. 
The  maximum  discrepancy  of  the  computed  wall  superheat 
between  Reilly's  data  and  the  present  data  is  about  twenty 
percent  at  the  maximum  heat-flux  setting.  This  figure  also 
shows  the  typical  uncertainty  of  about  10  percent  computed 
for  the  wall  superheat  values  during  this  investigation. 

B.  BOILING  PERFORMANCE  OF  THE  SMOOTH  TUBE 

Figure  4.4  shows  the  performance  of  the  smooth  copper 
tube  in  pure  R-114  and  ref ri gerant-oi 1  mixtures  of  3  and 
10  mass  percent  oil.  The  behavior  of  the  smooth  copper 
tube  in  pure  R-114  liquid  demonstrates  typical 
nucleate  pool-boiling  performance.  As  can  be  seen  in 
Figure  4.4.  the  line  between  points  A  and  B  represents 
a  region  of  constant  slope  (on  a  log-log  scale)  for 
the  variation  of  heat  flux  with  the  tube  wal 1  superheat 
(Two-Tsat).  This  section  of  the  curve  represents  natural 
convection  (i.e.,  no  bubbles  were  generated). 
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Upon  reaching  point  B  (see  Figure  4.4),  the  boiling 
tube  wall  superheat  begins  to  decrease  even  though  the 
heat  flux  is  continually  increased.  At  this  point,  the 
tube  begins  to  demonstrate  the  characteristics  of 
incipient  nucleate  boiling.  This  region  is  known  as  the 
mixed  boiling  region,  where  the  transition  from  natural 
convection  to  nucleate  boiling  occurs.  This  event  is 
characterized  by  the  increasing  number  of  active 
nucleation  sites  with  increasing  heat  flux.  It  should  be 
noted  that  during  all  data  runs,  the  unheated  ends  of  the 
tubes  did  not  show  any  nucleation  sites.  Reilly  CRef.  51 
also  stated  that  the  unenhanced  ends  of  the  boiling  tube, 
at  all  observed  heat  fluxes,  underwent  only  natural 
convection  despite  a  small  amount  of  heat  axially  conducted 
along  the  tube  walls.  It  was  observed  that  this  transition 
from  natural  convection  to  nucleate  boiling  occurs  very 
quickly,  usually  within  a  few  seconds  after  the  activation 
of  the  first  nucleation  site. 

At  the  heat  flux  represented  by  point  C  in  Figure 
4.4,  all  the  available  activation  sites  of  the  boiling  tube 
appear  to  be  active  and  the  wa 1 1  superheat  begins  to 
increase  with  increasing  heat  flux.  From  point  C  to 
point  D,  where  the  heat  flux  is  increased  to  a  maximum 
of  91  kW/m2 ,  new  activation  sites  are  generated  (maximum 
heat  flux  was  limited  by  cartridge  heater  rated  output). 


Figure  4.4  Heat-Transfer  Performance  for  Smooth  Tube  in 
0,  3  and  |0  Percent  Oil  . 


As  the  oil  concentration  was  increased,  there  was  a 
significant  increase  in  the  amount  of  foam  generated. 

After  obtaining  a  maximum  heat  flux  at  point  D,  the 
heat  flux  was  decreased  and  is  represented  as  the  region 
between  points  D  and  E  on  Figure  4.4.  This  region  of  the 
curve  represents  decreasing  wall  superheat;  however,  more 
activation  sites  that  were  created  at  the  higher  heat-flux 
settings  remain  active  to  lower  heat  fluxes.  This 
increased  range  over  which  the  boiling  tube  continues  to 
nucleate  results  in  an  increase  in  the  heat- transfer 
performance  since  the  wall  superheat  temperature  for  the 
nucleate  boiling  at  any  particular  heat  flux  is  less  than 
the  tube  wall  superheat  for  the  natural  convection  heating 
at  the  same  heat  flux. 

Stephan  [Ref.  9]  reported  that  the  effects  of  adding 
oil  introduces  a  mass -d i f f us i on  resistance  to  the  liquid 
and  therefore  lowers  the  heat- transfer  coefficient.  Reilly 
[Ref.  5],  in  his  discussion  of  his  smooth-tube  data  at  a 
saturation  temperature  of  -2.2°C  (Figure  4.5),  stated  that 
the  3  mass  percent  oil  and  10  mass  percent  oil 
concentrations  demonstrated  higher  wall  superheat 
temperature  than  that  observed  in  the  pure  R-114  liquid. 
However,  Reilly  observed  that  the  wall  superheat  for  the  3 
percent  oil  concentration  was  higher  than  that  for  the 
10  percent  oil  concentration.  He  attributed  this  to  the 


perce 


non-linear  physical  property  characteristics  of  the 


refrigerant- 

oil  mi xtures 

.  As  can 

be 

seen 

i  n 

Figure  4.4, 

the  general 

trend  of 

the  data 

i  s 

the 

same 

with  the 

except  ion 

of  the  3 

percent 

o  i  1 

concentration  wall 

superheat , 

which  near  1 y 

over  1  ays 

the 

data 

for 

pure  R-114 

liquid.  This  deviation  from  the  Reilly  data  may  be  due  in 
part  to  the  difference  in  the  saturation  temperatures.  The 
addition  of  oil  to  the  pure  R-114  liquid,  up  to  a  maximum 
of  10  percent,  should  delay  the  transition  from  natural 
convection  to  nucleate  pool  boiling  on  the  tube.  With  a 
r ef r i ger ant -o i 1  mixture  of  10  mass  percent  oil,  it  was 
observed  that  the  surface  nucleation  sites  did  appear  to 
spread  more  slowly  as  the  heat  flux  was  increased.  This 
condition  continued  until  the  tube  was  fully  engulfed  in 
the  nucleate  boiling  mode.  This  observation  was  also  noted 
by  Reilly  [Ref.  51  and  it  agrees  with  the  contention  of 
Chongrungr eong  and  Sauer  [Ref.  8]  and  Thome  [Ref.  6]  that, 
with  the  exception  of  the  surface  tension,  the  non-linear 
variations  in  the  physical  properties  of  the  refrigerant- 
oil  mixtures  explain  the  variation  in  the  boiling  heat 
transfer  of  the  tube. 

Figure  4.6  shows  the  heat - tr ans f er  coefficient  of  the 
smooth  tube  in  re f r i ge rant - o i 1  mixtures,  plotted  as  a 
function  of  heat  flux.  Figure  4.7  shows  the  same  plot  for 
Reilly's  data  for  a  saturation  temperature  of  -2.2°C 


q  <W/m2> 

Boiling  Heat-Transfer  Coefficient  for  Smooth  Tube 


[Ref.  5].  As  can  be  seen,  the  two  graphs  are  nearly 
identical,  representing  good  repeatability  of  the  data. 
Also,  the  effect  of  the  addition  of  up  to  6  percent  oil 
(Figure  4.7)  consistently  represents  an  approximate 
decrease  of  10  percent  in  the  heat- transfer  coefficient. 
As  the  oil  concentration  was  increased  from  6  to  10 
percent,  the  decrease  in  the  heat-transfer  coefficient  was 
a  maximum  of  about  35  percent.  These  decreases  in  the 
heat-transf er  coefficients  were  dependent  on  the  comparison 
of  the  pure  R-114  and  the  10  percent  oil  concentration 
at  equal  heat-flux  settings. 

C.  BOILING  PERFORMANCE  OF  THE  GEWA-T  TUBE 

Figure  4.8  shows  the  nucleate  pool-boiling  performance 
of  the  GEWA-T  (1.02  fins/mm)  tube  in  R-114-oil  mixtures. 
The  magnitude  of  the  tube-wall  superheat  was  considerably 
less  than  the  values  obtained  for  the  smooth  tube  discussed 
in  Section  IV. B.  This  lower  wall  superheat  can  be 
attributed  to  the  specially  formed  fins  of  the  GEWA-T  tube. 
The  performance  of  the  GEWA-T  tube  is  similar  to  that  of 
the  smooth  tube  with  the  curves  of  the  increasing  heat  flux 


to  the 

point 

of 

incipient 

nuc 1 ea  te 

bo  i  1  i  ng 

(point 

B)  . 

paralleling 

the 

performance 

of  the 

smooth 

copper 

tube 

(Figure 

4.4) 

.  The 

same  condition  exists  throughout 

t  he 

entire  range  of  decreasing  heat  fluxes. 
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The  region  bounded  by  points  A  and  B  shows  the 
characteristic  natural -convection  heating  behavior.  The 
onset  of  incipient  nucleate  boiling  at  point  B  is  clearly 
defined  for  the  pure  R-114  liquid.  However,  as  the  oil 
concentration  was  increased,  the  point  of  incipient 
nucleate  boiling  becomes  less  clearly  defined  and  the 
region  of  mixed  boiling  (point  B  to  C)  increased  with 
increased  heat  flux.  It  was  also  observed  that  in  the 
region  from  A  to  B,  the  3  percent  oil  mixture  shoved  the 
smallest  wall  superheat  and  the  pure  refrigerant  showed  the 
highest  wall  superheat.  The  reason  for  this  behavior  is 
not  fully  known  at  present,  but  may  be  due  to  convection 
effects  within  the  channels  of  the  GEWA-T  tube  which  are 
enhanced  by  the  foaming  action  in  the  presence  of  oil. 

The  region  bounder1  by  the  points  C  and  D  for  the 
increasing  heat-flux  condition  is  indicative  of  fully 
developed  nucleate  pool  boiling.  The  region  between  points 
D  and  E  represents  nucleate  pool  boiling  for  decreasing 
heat  flux.  Figure  4.8  shows  that,  as  the  oil 
concentration  was  increased,  the  wall  superheat  increased 
at  a  given  heat  flux.  This  is  contrary  to  the  results  for 
the  condition  of  increasing  heat  flux  in  the  region  of 
na tura 1 -convect i on  heating. 

When  comparing  the  performance  of  the  smooth  tube 
(Figure  4.4)  to  that  of  the  GEWA-Tube  (Figure  4.8),  the 
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Figure  4.8  Heat-Transfer  Performance  for  GEWA-T  Surface. 


maximum  adverse  effect  of  the  increased  oil  content  for  the 


smooth  tube,  occurs  at  the  higher  heat-flux  settings.  For 
the  GEWA-T  tube,  the  maximum  decrease  in  the  tube 
performance  for  increased  oil  content,  occurs  in  the  lower 
and  intermediate  heat-flux  range. 

Figure  4.9  shows  the  heat - t r ans f e r  coefficient  as  a 
function  of  heat  flux  over  the  range  of  R-114-oil 
mixtures  (0,  1,  2,  3,  6  and  10  percent  oil).  Data  observed 
over  the  range  of  R-114-oil  mixtures  of  1,  2,  3  and  6 
percent  oil  are  nearly  coincident.  The  10  percent  oil 
concentration  data  represent  an  approximate  30  percent 
reduction  of  the  heat-transfer  coefficient  when  compared  to 
the  heat-transfer  coefficient  obtained  during  the  data  for 
pure  R-114.  "Figure  4.10  shows  the  direct  comparison  of  the 
GEWA-T  tube  with  the  smooth  tube  for  R-114-oil  mixtures  of 
0,  3  and  10  percent  oil.  The  GEWA-T  tube  represents  an 


overa  1  1 

enhancement  of  4.2 

over 

the  base  1 i ne 

smoot  h 

copper 

tube  at 

a  heat  flux  of  40 

kW/m* 

with 

pure 

R-114 

and  an 

overa  1  1 

enhancement  of  3.8 

with 

10 

percent 

oil. 

These 

r e  su 1 ts 

compare  favorably 

with 

t  he 

r esu 1 ts 

o  f 

Y  i  1  maz  , 

Hwalek  and  Westwater  [Ref.  181  for  a  12.3-mm  OD  GEWA-T  tube 
in  p-xylene,  in  which  they  obtained  an  enhancement  of  5.3 
in  the  pool-boiling  hea t - t r ans f e r  coefficient  over  their 
baseline  smooth  copper  tube.  Additionally,  Yilmaz  and 
Westwater  [Ref.  123  tested  a  similar  surface  in  isopropyl 


alcohol  in  which  an  enhancement  of  2.0  was  observed.  Marto 


and  Lepere  CRef.  133  tested  a  17.9-mm  OD  GEWA-T  surface  in 
both  R-113  and  FC-72  and  reported  enhancements  of  2.8  and 
2.5.  respectively. 

Figure  4.11  shows  more  clearly  the  degradation 
that  oil  causes  in  the  boiling  hea t - t r ans f er 
performance  for  the  GEWA-T  tube.  This  figure  plots 
the  hea t - t r ans f e r  coefficient  of  the  GEWA-T  tube  in 
the  oil  mixture  relative  to  the  heat- transfer 
coefficient  in  pure  R-114  as  a  function  of  oil 
concentrat i on .  From  the  shape  of  the  curves,  the 
oil-caused  degradation  does  not  show  a  clear  trend 
with  the  oil  concentration  at  all  practical  heat 


fluxes. 

Over 

the 

range 

o  f 

oil 

concentrations,  except 

at  a 

heat 

flux 

setting 

of 

37 

kW/m2  ,  there  is  only 

a  small  degradation  (up  to  a  maximum  of  35%)  in  the 
boiling  performance.  This  is  probably  due  to  the 
large  channels  and  reentrant  cavities  in  the  tube, 
al lowing  the  "pumping"  action  of  the  bubbles  to  continue 
to  remove  the  oil-rich  liquid  throughout  the  entire 
heat-flux  range  without  a  significant  reduction  in 
performance.  The  only  significant  decrease  in 
performance  occurred  at  heat  fluxes  greater  than  5 
kW/m2  with  10  percent  oil. 
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Figure  4.9  Boiling  Heat -Transfer  Coefficient  for  GEWA-T  Surface. 


Boiling  Heat  Transfer  Coefficient  Comparison  for  GEWA  T 
Surface  and  Smooth  Surface  . 


Figure  4.11  Relative  Effect  of  Oil  on  GEWA-T  Boiling  Heat-Transfer 
Performance . 


D. 


BOILING  PERFORMANCE  OF  THE  THERMOEXCEL-E  TUBE 


Figure  4.12  shows  the  nucleate  pool-boiling 
performance  of  the  Thermoexce 1 -E  tube  in  R-114-oil 
mixtures.  As  can  be  seen  in  the  figure,  the  region  of 
natura 1 -convect i on  heating  for  the  condition  of  increasing 
heat  flux,  points  A  to  B,  for  the  pure  R-114  is  very  small, 
with  the  point  of  incipient  nucleate  boiling  occurring  at  a 
reasonably  low  heat  flux  (1.5  kW/m2 ) .  At  this  point  of 
incipient  nucleate  boiling,  the  tube  performs  in  the  mixed 
boiling  region  as  a  limited  number  of  nucleation  sites 
become  active.  From  point  C  to  C',  the  tube  appears  to 
operate  in  more  of  a  natura 1  - convect i on  heating  mode  as  no 
new  nucleation  sites  are  activated,  although  previously 
activated  sites  continue  to  nucleate.  This  phenomenon  is 


probab 1 y 

due  to 

the 

i nab i 1 i ty 

of  the 

cartridge  heaters 

to 

pr  oduce 

uni  form 

heat  flux  at 

these 

1  ower 

sett ings . 

At 

point  C ’ 

with 

a 

hea  t - f lux 

setting 

of  5 

kW/m2 ,  more 

nuc 1  eat i on 

sites 

are 

generated 

and  as 

the 

heat  flux 

i  s 

increased 

to  60  kU/m2 

,  the  tube 

becomes 

fully 

nuc 1 eated . 

For  the  3  percent  oil  concentration,  the  point  of 
incipient  nucleate  boiling,  point  B',  decreases  as  the  non¬ 
linear  physical  properties  of  the  R-114-oil  mixture  allow 
the  onset  of  the  mixed  boiling  region  to  occur.  Also,  the 
region  of  the  mixed  boiling,  point  B'  to  D'  ,  is 
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far  larger  than  that  for  pure  R-114,  as  the  generation  of 
new  nucleation  sites  is  retarded  by  the  increased  oil 


content . 

For  the  10  percent  oil  concentration,  on  the  other 
hand,  the  point  of  incipient  nucleate  boiling,  point  B’*, 
is  delayed  considerably.  At  approximately  14  kW/m* ,  the 
transition  from  natural -convective  heating  to  nucleate 
boiling  occurs.  The  exact  reason  for  the  difference  between 
the  3  percent  and  10  percent  behavior  is  not  known,  however 
it  is  possibly  due  to  the  change  in  the  physical  properties 
of  the  R-114-oil  mixture. 

The  region  of  nucleate  pool  boiling  for  the  condition 
of  decreasing  heat  flux  represented  by  points  D  to  E.  D' 
to  E'  and  D’ '  to  E' ’  for  0,  3  and  10  percent  oil 
concentrations,  respectively,  for  the  Thermoexce 1 -E  tube, 
represent  classic  nucleate  pool-boiling  characteristics. 
As  the  oil  concentration  increases,  the  wall  superheat  also 
increases  throughout  the  range  of  heat-flux  settings.  At 
the  lower  heat  fluxes,  the  largest  increase  in  the  wall 
superheat  occurred  at  3  percent  oil  concentration  when 
compared  to  the  pure  R-114,  while  the  10  percent  oil 
concentration  was  observed  to  undergo  continual  degradation 
in  tube  performance  over  the  entire  range  of  heat-flux 


kW/m2  ,  the  R-114-oil  mixture  of  10  percent  oil  was  observed 
to  undergo  a  reduction  in  tube  performance  by  a  factor  of 
3.  For  the  practical  heat-flux  range  of  30  to  40  kW/m2 , 
there  is  an  increase  in  wall  superheat  by  factors  of  1.4 
and  2.6  for  R-114-oil  mixtures  of  3  and  10  percent, 
respective  1 y . 

Figure  4.13  shows  the  heat-transf er  coefficient  of  the 
Thermoexce 1 -E  tube  in  refrigerant-oil  mixtures,  plotted  as 
a  function  of  heat  flux.  As  can  be  seen  in  the  figure,  the 
heat-transfer  coefficient  remains  nearly  constant 
throughout  the  range  of  heat-flux  settings  for  pure  R-114, 
while  there  are  slight  reductions  for  the  1,  2  and  3 
percent  oil  concentrations.  Additionally,  it  can  be  seen 
that  at  these  lower  oil  concentrations,  the  heat - t r ans f er 
coefficients  are  nearly  coincident.  As  the  R-114-oil 
mixture  oil  concentration  is  increased  to  6  percent,  the 
tube  performance  increases  with  increasing  heat  flux  until 
reaching  a  heat-flux  setting  of  60  kW/m2 ,  at  which  time 
tube  performance  begins  to  deteriorate.  This  is  contrary 
to  the  tube  performance  for  0,  2,  and  3  percent  oil 
concentrations  in  which  the  heat-transf er  coefficient  was 
observed  to  increase  slightly  over  the  range  of  heat 
fluxes.  For  the  10  percent  oil  concentration,  the  heat- 
transfer  coefficient  decreased  as  the  heat  flux  was 


i ncr eased . 
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Transfer  Coefficient  for  Thermoexcel-E  Surface 


Figure  4.14  shows  a  direct  comparison  of  the 
Thermoexce 1 -E  tube  to  the  smooth  copper  tube  for  R-114-oil 
mixtures  of  0,  3  and  10  percent.  The  Thermoexce 1 -E  tube 
presents  enhancements  of  6.6  and  4.1  over  the  performance 
of  the  smooth  copper  tube  in  pure  R-114  and  an  R-114-oil 
mixture  of  10  percent  oil,  respectively,  at  a  heat-flux 
setting  of  40  kW/m2 .  The  results  for  the  pure  R-114  compare 
favorably  with  the  results  obtained  by  Yilmaz,  Palen  and 
Taborek  [Ref.  19]  for  the  Thermoexce 1 -E  in  p-xylene. 

Figure  4.15  shows  the  degradation  that  oil  causes  in 
the  boiling  heat-transf er  performance  for  the  Thermoexce 1 -E 
tube.  This  figure  plots  the  ratio  of  the  heat - 1 r ans f er 
coefficient  of  the  Thermoexce 1 -E  tube  with  oil  to  the  value 
without  oil  as  a  function  of  oil  concentration.  The  shapes 
of  the  curves  continue  to  demonstrate  the  non-linear 
variation  of  the  physical  properties  of  refrigerant-oil 
mixtures.  However,  the  significant  degradation  of  the  tube 
performance  is  probably  due  to  the  inability  of  the 
generated  bubbles  to  scavenge  the  interior  tunnels  of  the 
Thermoexce 1 -E  tube.  This  is  in  contrast  to  the  GEWA-T 
tube.  where  more  scavenging  is  possible  due  to  its  larger 
channe 1  size. 

E.  BOILING  PERFORMANCE  OF  THE  THERMOEXCEL-HE  TUBE 

Figure  4.16  shows  the  nucleate  pool-boiling 

performance  of  the  Thermoexce 1 -HE  tube  in  R-114-oil 
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Boiling  Heat-Transfer  Coefficient  for  Thermoexcel-E  vs 
Smooth  Surface . 


Figure  4.15  Relative  Effect  of  Oil  On  TTiermoexcel-E  Boiling 
Heat-Transfer  Performance. 


mixtures.  In  the  region  between  points  A  and  B/B', 
natural -convection  heating  occurs  for  increasing  heat  flux. 
Although  similar  in  performance  to  the  Thermoexce 1 -E  tube, 
the  Thermoexce 1 -HE  tube  presented  a  more  linear 
characteristic  than  that  of  the  Thermoexce 1 -E  tube, 
representing  a  more  uniform  heat  distribution.  At  point  B 
on  the  figure,  the  point  of  incipient  nucleate  boiling  for 
pure  R-114  liquid  occurs  at  approximately  2  kW/m2 .  For 
R-114-oil  mixture  of  3  and  10  percent  oil  concentrations, 
the  point  of  incipient  nucleate  boiling,  point  B'.  occurs 
at  approximately  8  kW/m2 .  Also  in  this  region,  the  pure 
R-114  and  the  R-114-oil  mixture  of  3  percent  oil 
concentration  are  nearly  coincident  in  the  range  of  heat 
flux  from  500  W/m2  to  2  kW/m2 . 

From  points  B  to  C,  the  Thermoexce I -HE  tube  in  the 
pure  R-114  liquid  achieved  full  nucleate  boiling  at 
approximately  37  kW/m2 ,  with  intermediate  data  points  in 
the  region  of  mixed  boiling.  These  intermediate  data 
points  indicate  that  the  Thermoexce 1 -HE  tube  was  not  able 
to  activate  all  nucleation  sites  simultaneously  for 
increasing  heat-flux  conditions.  The  R-114-oil  mixture  of 
10  percent  oil  was  able  to  achieve  activation  of  all 
nucleation  sites  upon  reaching  22  kW/m2 ,  without  any  data 
points  contained  in  the  mixed  boiling  region.  With  3 
percent  oil,  the  Thermoexce 1 -HE  tube  was  able  to  achieve 


f  u  1  I 


nucleate  boiling  at  37  kW/m2 . 


The  Thermoexce 1 -HE 


tube  in  the  R-114-oil  mixtures  appeared  to  activate 
all  nucleation  sites  simultaneously. 


The 

region  from  points 

C  to  D 

and 

from  C*  to  D' 

represent 

the  nucleate  pool 

bo i 1 i ng 

for 

increasing  heat 

flux.  The  performance  of  the  Thermoexce 1 -HE  tube  in  this 
region  appears  to  decrease  as  the  oil  content  was 
increased.  The  10  percent  oil  concentration  increased  the 
wall  superheat  by  a  factor  of  approximately  2  above  the 
wall  superheat  for  the  pure  Rrll4  liquid  for  heat-flux 
settings  of  40  kW/m2 . 

The  regions  between  the  points  D  and  E  and  D*  and  E' 
represent  nucleate  pool-boiling  for  decreasing  heat-flux 
conditions.  As  can  be  seen  from  the  figure,  the  greatest 
reduction  in  performance  occurs  in  the  lower  heat-flux 
settings  between  the  pure  R-114  and  the  3  percent  oil 
concentration.  This  is  similar  to  the  results  for  the 
Thermoexce 1 -E  tube,  where  the  tube  performance  decreased 
as  oil  content  increased.  There  is  a  significant  decrease 
in  the  Thermoexce 1 -HE  tube  performance  over  the  range  of  - 
heat-flux  settings  as  the  oil  concentration  is  increased  to 
10  percent. 

Figure  4.17  shows  the  heat-transfer  coefficient  of  the 
Thermoexce 1 -HE  tube  in  ref rigerant-oi 1  mixtures,  plotted 
as  a  function  of  heat  flux.  As  can  be  seen,  the 
performance  is  very  similar  to  that  of  the  Thermoexce 1 -E 


tube  (Figure  4.13). 


As  was  observed  for  the  Thermoexce 1 -E 


tube,  the  Thermoexce 1 -HE  tube  undergoes  a  slight  reduction 
in  the  heat - t rans f er  coefficient  as  the  oil  concentration 


was  increased  from  1 

to 

3  percent.  As 

the 

o  i  1 

concentrat i on 

is  increased 

to  6 

and  10  percent 

oil. 

t  he 

performance 

reduction  is 

qu  i  te 

noticeable.  With 

the 

10 

percent  o  i  1 

concentration, 

tube 

performance  decreases 

at 

both  the  lower  and  higher 

heat 

-flux  settings. 

Th  i  s 

i  s 

contrary  to  the  findings  for  the  Thermoexce 1 -E  tube  where 
the  performance  decreased  over  the  entire  range  of  heat- 
flux  settings,  from  low  heat  flux  to  high  heat-flux 
settings. 

Figure  4.18  shows  a  comparison  of  the  Thermoexce 1 -HE 
tube  to  the  smooth  copper  tube  at  R-114-oil  mixtures  of  O, 
3  and  10  percent  oil  concentrations.  The  Ther moex ce 1 -HE 
tube  displayed  an  enhancement  of  6.1,  4.6  and  4.1  at  a  heat 
flux  setting  of  40  kW/m2  and  R-114-oil  mixtures  of  0.  3  and 
10  percent  oil,  respectively,  over  the  same  parameters  for 
the  smooth  copper  tube.  These  results  are  nearly  identical 
to  those  found  for  the  Thermoexce 1 -E  tube. 

Figure  4.19  shows  the  degradation  that  the  oil  has  on 
the  boiling  heat-transfer  performance  of  the  Thermoexce 1 -HE 
tube.  At  the  lower  heat  fluxes,  the  effect  of  the  oil 
concentration  is  slightly  more  pronounced  than  that 
observed  for  the  Thermoexce 1 -E  tube  throughout  the  range  of 
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Transfer  Coefficient  for  Thermoexcel-lfE 


Boiling  Heat-Transfer  Coefficient  Comparison  for 
Thermoexcel-HE  Surface  vs.  Smooth  Surface  . 


oil  concentrations,  although 


more  severe  at  the  lower  oil 


concentrations.  At  •  the  higher  heat  fluxes,  the  trend 
appears  to  reverse. 

F.  COMPARISON  TO  REILLY'S  DATA  FOR  THE  HIGH  FLUX  TUBE 

Figures  4.20,  4.21  and  4.22  provide  a  comparison  of 
the  smooth,  High  Flux,  GEWA-T,  Thermoexce 1 -E  and 
Thermoexce 1 -HE  tubes  for  0,  3  and  10  percent  oil 

concentrations,  respectively.  For  pure  R-114  liquid,  the 
High  F 1 ux , Thermoexce 1 -E  and  Thermoexce 1 -HE  tubes  operating 
at  a  heat  flux  of  approximately  20  kW/nt* ,  show  nearly 
identical  performance.  As  the  heat  flux  was  increased,  the 
performance  of  the  High  Flux  tube  continued  to 

increase  while  the  performance  of  Thermoexce 1 -E  and 

Thermoexce 1 -HE  tubes  remained  nearly  constant.  At  heat- 
flux  settings  lower  than  20  kW/m2 ,  the  Thermoexce 1 -E  and 
Thermoexce 1 -HE  tubes  showed  higher  heat-transfer 
coefficients  than  that  for  the  High  Flux  tube.  As  the  oil 
concentration  was  increased  to  3  percent,  the  performance 
of  all  five  tubes  decreased,  as  can  be  seen  in  Figure  4.21. 
Additionally,  the  GEWA-T  tube's  performance  appeared  to 
decrease  the  least  since,  as  discussed  earlier  in  section 
IV. C,  the  finned  surface  of  the  GEWA-T  tube  is  not 
adversely  effected  by  the  small  increase  in  the  oil  content 
due  to  its  large  reentrant  cavities.  From  this 

observation,  it  can  be  stated  that  the  performance  of  the 


Comparison  of  Boiling  Heat-Transfer  Coefficients  for 
Smooth,  High  Flux  [Ref.  5]  ,  GHWA-T,  Thermoexcel-E 
and  Thermoexcel -HE  Surfaces  for  3  Percent  Oil  Content 


tubes  with  porous  coatings  or  small  reentrant  cavities  is 


adversely  affected  by  the  increase  in  oil  content  due  to 
the  generation  of  an  oil-rich  film  in  the  nucleation  sites 
impeding  the  tube's  performance.  This  is  further  evidenced 
by  Figure  4.22,  which  depicts  the  performance  of  all  five 
tubes  in  R-114-oil  mixtures  of  10  percent  oil.  The 
performance  of  all  the  tubes  is  reduced  significantly, 
while  the  performance  of  the  High  Flux  tube  is  severely 
decreased  in  the  range  of  heat-flux  settings  from  22  kW/m2 
to  the  maximum  setting  of  91  kW/m2 .  At  approximately  45 
kW/m2 ,  the  performance  of  the  High  Flux  tube,  the  GEWA-T. 
the  Thermoexce  1  .-E  and  the  Ther  moexce  1  -ME  tubes  are 
approximately  equal,  although  the  performance  of  the 
Thermoexce 1 -E  and  the  Thermoexce 1 -HE  tubes  remain  nearly 
constant  over  the  range  of  heat-flux  settings  of  14  kW/m2 
to  50  kW/m2 .  Although  the  performance  of  the  High  Flux 
tube  exceeds  that  of  the  GEWA-T,  the  Thermoexce 1 -E  and  the 
Thermoexce 1 -HE  tubes  over  this  range,  its  performance 
begins  to  severely  degrade  at  the  midpoint  of  the  heat-flux 
range  indicated.  In  the  heat-flux  range  above  45  kW/m2 , 
the  GEWA-T  tube  outperformed  all  the  other  tubes.  This  is 
attributed  to  the  ability  of  bubbles  generated  in  its 
larger  channels  to  continually  scavenge  the  oil-rich  layer 
from  the  tube,  and  therefore  remove  the  insulating  oil  from 
the  tube's  channels  and  exterior  surface. 
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V .  CONCLUSIONS  AN D  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  In  pure  R-114,  at  a  heat  flux  of  30  kW/m2  the 

pool-boiling  heat-transfer  coefficient  of  the 
Thermoexce 1 -E ,  Thermoexce 1 -HE  and  GEWA-T  tubes 
were  approximately  6.6,  6.1  and  4.2  times, 

respectively,  larger  than  the  baseline  smooth  tube. 

2.  In  R-114-oil  mixtures  of  3  percent  oil,  at  a 

heat  flux  of  30  kW/m2 ,  the  pool-boiling  heat- 
transfer  coefficient  of  the  Thermoexce 1 -E,  Thermo- 
excel-HE  and  GEWA-T  tubes  were  approximately  4.8, 
4.6  and  4.0  times,  respectively,  larger  than  the 
smooth  tube  under  the  same  conditions. 

3.  In  R-114-oil  mixtures  of  10  percent  oil,  at  a  heat 

flux  of  30  kW/m2 ,  the  performance  of  the  Thermoexcel- 
E  'and  Thermoexce 1 -HE  tubes  were  nearly  identical, 
with  an  enhancement  of  4.1  over  the  pool-boiling 
heat- transfer  coefficient  of  the  smooth  tube.  The 
GEWA-T  tube,  under  the  same  conditions,  produced  an 
enhancement  of  3.8  over  the  smooth  tube. 

4.  As  oil  content  in  the  refrigerant-oil  mixture 

increased,  the  performance  of  enhanced  surfaces 
decreased;  however,  the  amount  of  this  decrease  was 
directly  related  to  the  size  of  the  subsurface 
channels  (i.e.,  the  larger  the  reentrant  cavity  or 
channel,  the  less  adverse  effect  of  the  increased  oil 
content ) . 

5.  As  the  oil  content  increases  in  the  channels,  the 
greater  the  adverse  effect  of  the  oil-rich  layer  has 
on  the  surface,  insulating  the  surface  from  the 
refrigerant  and  therefore  decreasing  the  performance 
of  the  surface. 

6.  The  GEWA-T  surface  in  refrigerant-oil  mixtures,  with 
its  larger  subsurface  channels,  was  able  to 
scavenge  the  oil-rich  layer  from  its  cavities  by 
the  bubbling  action  generated  within  the  channels. 

7.  The  performance  of  the  Thermoexce 1 -E  and  Thermoexcel- 
H£  tubes,  with  their  smaller  pore  size  and  smaller 
tunnels,  showed  significant  performance  decreases 
as  oil  content  was  increased.  The  smaller  pore  and 
tunnel  size  under  these  conditions,  were  unable  to 
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effectively  disperse  the  oil-rich  film  from  the  sub 
-surface  cavities  and  therefore  allowing  the  oil-rich 
film  to  insulate  the  tube’s  surface. 


B.  RECOMMENDATIONS 

1.  The  single-tube  testing  should  be  expanded  to  include 
multi-tube  bundles.  This  would  provide  more  realistic 
data  with  regard  to  the  effects  of  tubes  operating  in 
close  proximity  to  each  other. 

2.  A  uniform  heat  flux  is  necessary  within  the  boiling 
tube  and  within  the  boiling  section  itself.  There¬ 
fore,  the  present  cartridge  heaters  need  to  be 
replaced  with  more  precise  heaters,  capable  of  pro¬ 
ducing  a  uniform  heat  flux  in  the  axial  direction. 

3.  The  testing  of  surfaces  should  be  expanded  to  study 
the  effects  of  various  refrigerant  oils. 

4.  Conduct  experiments  to  establish  the  physical  proper¬ 
ties  of  r ef r i gerant-o i  1  mixtures  to  enable  further 
understanding  into  the  boiling  mechanisms  involved. 

5.  There  needs  to  be  additional  studies  undertaken  on 
the  effects  of  additional  boiling  temperatures  on  the 
mixture  properties  of  the  r ef r i ger ant -o i 1  mixtures. 


APPENDIX  A 

THERMOCOUPLE  CALIBRATION  [Ref.  5] 

Karasabun  [Ref.  17]  described  the  thermocouple 
calibration  procedure  in  great  detail.  Reilly  [Ref.  5] 
stated  that  the  da t a - r educ t i on  program  utilized  differences 
in  thermocouple  readings  in  all  computations,  such  as  the 
amount  of  the  tube  wall  superheat  (Two  -  Tsat).  Reilly 
stated  that  the  calibration  of  the  thermocouples  would  be 
necessarv  if  the  thermocouples  were  measuring  quantities 
such  as  absolute  temperature  to  determine  actual  fluid 
properties.  The  procedure  utilized  by  Karasabun  and  Reilly 
for  the  initial  instrumentation  for  the  experimental 
apparatus  is  as  contained  in  Appendix  A  of  Reference  5. 

To  determine  the  calibration  for  the  thermocouples, 
two  thermocouples  were  made,  one  from  the  beginning  of  the 
wire  rol I  and  the  second  from  the  end.  The  seventh-order 
polynomial  equation  provided  by  the  wire  manufacturer  was 
corrected  by  adding  to  it  a  second-order  curve  fit  of  the 
variation  of  the  manufacturer's  commuted  temperature  for  a 
given  emf  from  a  known  set  of  reference  temperatures  as 
measured  by  a  Hewlett-Packard  2804A  quartz  thermometer.  The 
thermometer  had  a  temperature  resolution  of  ±0.0001  K  and 
an  accuracy  of  ±0.03  K. 


The  manufacturer's  emf  to  temperature  conversion 
equat ion  is: 

T  =  a0  ♦  a,  E  +  a-  E2  *  a3  E3  +  a4  E4  +  (A.l) 

a3  Es  +  a4  E*  +  a?  E7 

where 

E  =  thermocouple  reading  in  volts 

T  =  temperature 

a„  =  0.100860910 

a,  =  25727.94369 

a,  =  -767345.8295 

a3  =  78025595.81 

a*  =  -9247486589 

a,  =  6.97688x10*  1 

a*  *  -2. 66192x10* J 

a-,  =  3. 94078x10*  4 

Figure  A. 1  [Ref.  5]  shows  the  difference  in  the  quartz 
thermocouple  readings  and  the  actual  thermocouple  readings 
(discrepancy)  plotted  as  a  function  of  temperature.  As  can 
be  seen  from  the  figure,  the  two  thermocouples  agreed  to 
within  0.05  K  of  each  other.  The  manufacturer's 

calibration  equation  needed  an  increase  of  approximately 
0.1  K  to  more  accurately  convert  the  emf's  to  the  true 
temperature.  The  second-order  polynomial  correction  was  as 
f o 1  lows; 

DCP  =  b0  4  b,T  +  b2Tz  (A.  2) 


where 


DCP 


discrepancy  (K) 


b0  =  8. 6268968x10- 2 

b,  =  3. 7619902x10- s 

b,  =  -5. 0689259x10' s 
T2  =  thermocouple  readings  from 
equa t i on  A . 1  in  C 

Therefore.  the  temperatures  computed  by  the  da ta- r educt i on 
program  (DRP5)  were  emf's  converted  directly  to 
temperatures  by  equation  A.l  and  corrected  by  equation  A. 2. 
This  calibration  was.  thus  estimated  to  result  in 
temperature  readings  within  10.05  K  of  the  true 
temperature . 
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APPENDIX  B 


DATA-REDUCT I  ON  PROGRAM  CREF.  5] 


% 


X 


S 

: 

$ 


The  data-reduct i on  program  below  consists  of  the 
following  sections: 

Main  Program  -  Menu  of  subprogram  options 
Sub  Main  -  Takes  or  reprocesses  data 

Sub  plot  -  Plots  data  on  log-log  scale 

Sub  Poly  -  Computes  least  squares  curve  fit 

Sub  Plin  -  Plots  data  on  1 i near  -  1 i near  scale 

Sub  Stats  -  Computes  average  and  standard 

deviation  of  data 

Subprogram  Main  consists  of  the  following  steps: 

1.  Creates  data  file  for  data  and  plotting  file. 

2.  Annotates  tube  type. 

3.  Monitors  heat  flux  or  saturation  temperature  to 
establish  steady-state  conditions. 

4.  Scans  HP3497  channels  and  stores  in  data  file. 

5.  Converts  raw  emf's  to  temperatures,  current  and 
voltages. 

6.  Computes  the  heat- transfer  rate  for  the  tube’s 
internal  cartridge  heater. 

7.  Computes  the  average  wall  temperature  of  the 

tube,  the  average  wall  superheat  and 


the  film 


temperature, 


ft 

i 


8.  Computes  the  physical  properties  of  the 
R-114  using  the  given  correlation  temperature. 

9.  Computes  the  natural -convection  heat- transfer 
coefficient  of  R-114  for  the  non-boiling 
tube  ends. 

10.  Computes  the  heat  loss  from  the  non-boiling 
tube  ends. 

11.  Calculates  the  corrected  heat  fluxes  from 
the  testing  tube  to  the  R-114. 

12.  Calculates  the  boiling  heat- transfer  coefficient 
of  the  R-114  from  the  tube  being  tested. 

13.  Prints  the  data  and  stores  the  heat  flux  and  tube 
wall  superheat  values  in  the  plot  file. 

The  following  is  the  listing  of  the  complete  data- reduct i on 
program  (DRP5)  written  in  Basic  3.0  for  the  Hewlett-Packard 


9826  computer. 


1000!  FILE  NAME :  0RP5 
1003!  DATE:  October  19.  1984 

1006!  REVISED:  April  25.  1986 
1009! 

1012  COM  /Idp/  Ido 
1 01 S  PRINTER  IS  1 
1018  CALL  Select 

1021  INPUT  ‘WANT  TO  SELECT  ANOTHER  OPTION  ( 1 -Y.O-N)?* .Isei 

1024  IF  Is# 1“ 1  THEN  GOTO  1018 

1027  8EEP 

1030  SEEP 

1033  PRINTER  IS  1 

1036  PRINT  “DATA  COLLECTION/REPROCESSING  COMPLETED* 

1 039  END 

1042  SUB  Main 

1045  COM  /Idp/  Idp 

1048  COM  /Cc /  C(7).Ical 

1051  COM  /Mil/  D2.Di.0o.L.Lu.Kcu 

1054  DIM  Ent'(  l2>.T(12),D1a(1U  .D2a(  1 1 )  .Dial  1 1 )  ,Doa(  11 )  .La(  1 1 )  ,Lua(  1 1 )  ,Kcua(  1 1  > , 
Et( 19) , TnS< 4) C 151 

1057  DATA  0.10086091.25727.94369.-767345.3295.78025595.81 
1060  DATA  -9247486589. 6. 97688E-H1, -2. 66192E+13, 3. 94078E+14 
1063  READ  C(»> 

10661  DATA  “Smooth" . "High  Flux*  ,  “Thern»o#x#l-E*  ,  “Thermoexel-HE’ 

1069  DATA  Saooth.High  Flux , Turbo-B .High  Flux  Mod.Turbo-B  Mod 
1072  READ  Tn®<*> 

1075  PRINTER  IS  701 
1078  BEEP 

1081  IF  Idp-4  THEN  PRINTER  IS  1 
1084  IF  Ido-4  THEN  GOTO  1774 

1087  INPUT  “ENTER  MONTH.  DATE  AND  TIME  <MM:DD:HH:MM:SS) * .Date® 

1090  OUTPUT  709: “TD“ :Date® 

1093  OUTPUT  709: *TD* 

1096  ENTER  709;Dat#S 
1099  PRINT 

1102  PRINT  *  Month,  date  and  time  :*;Date® 

1105  PRINT 

1108  PRINT  USING  "1  OX, ““NOTE:  Program  name  :  DRP4**“ 

1111  BEEP 

1114  INPUT  “ENTER  DISK  NUMBER". Dn 

1117  PRINT  USING  “16X.““Disk  number  -  “*,ZZ*:Dn 

1 1 20  BEEP 

1123  INPUT  “ENTER  INPUT  MODE  (0-30S4A. !-FILE)*,Im 
1126  BEEP 

1129  INPUT  “SELECT  HEATING  MODE  (0-ELEC:  1 -MATER) *, Ihm 
1132  BEEP 

1135  INPUT  “ENTER  THERMOCOUPLE  TYPE  (0-NEM. I -OLD) *  .leal 
1138  IF  Ia-0  THEN 
1141  BEEP 

1144  INPUT  “GIVE  A  NAME  FOR  THE  RAH  DATA  FILE*.D2  file® 

1147  PRINT  USING  “16X.““Ne«  file  name:  “* .  1 4A“ : 02_f i le® 

1150  Sizel-20 

T 1  S3  CREATE  BOAT  02  file®. Size! 

1156  ASSIGN  §File2  TO  Q2_file® 

1159! 

1162!  DUMMY  FILE  UNTIL  Nrun  KNOHN 
1165  D1  f ileS-‘DUMMY“ 
t 1 63  CREATE  BOAT  D1  fileS.Sizel 
1171  ASSIGN  4File1  TO  Dl.file® 

1174  OUTPUT  SFilel :Date® 


1177  IF  Ihm-0  THEN 
1180  BEEP 

1183  INPUT  “ENTER  NUMBER  OF  DEFECTIVE  TCS  ( O-DEFAULT) *.Idtc 
1186  IF  Idtc-0  THEN 
1189  Ldtcl-0 
1192  Ldtc2-0 

1 195  PRINT  USING  “16X,**No  defective  TCs  exist* “* 

1 1 98  END  IF 

1201  IF  Idtc-1  THEN 

1204  BEEP 

1207  INPUT  “ENTER  DEFECTIVE  TC  LOCATION* .Ldtct 

1210  PRINT  USING  “16X.““TC  is  defective  at  location  **,D*;Ldtc1 

1213  Ldtc2-0 

1216  END  IF 

1219  IF  Idtc-2  THEN 

1222  BEEP 

122S  INPUT  “ENTER  DEFECTIVE  TC  LOCATIONS* .Ldtct ,Ldtc2 

1228  PRINT  USING  *16X,"*TC  are  detective  at  locations  ** ,D,4X,D* iLdtcl ,Ldtc2 

1231  END  IF 

1234  IF  Idtc>2  THEN 

1237  BEEP 

1240  PRINTER  IS  1 

1243  BEEP 

1246  PRINT  ‘INVALID  ENTRY* 

1249  PRINTER  IS  701 
1252  GOTO  1180 
1255  END  IF 
1258  END  IF 

1261  OUTPUT  SFilef iLdtc! .Ldtc2 
1264!  Ib-1  option 
1267  ELSE 
1270  3EEP 

1273  INPUT  “GIVE  THE  NAME  OF  THE  EXISTING 'DATA  FILE* ,D2_f ileS 
1276  PRINT  USING  “16X.*“01d  file  nape:  “ * , 1 4A* ;D2_f 1 le* 

1279  ASSIGN  3File2  TO  D2_fileS 
1282  ENTER  3File2:Nrun 
1285  ENTER  $Fiie2:DoldS 

1288  PRINT  USING  “16X.“*This  data  set  taken  on  :  “*, I4A*;DoldS 
1291  ENTER  »Fiie2:Ldtc1 .Ldtc2 
1294  IF  Ldtc! >0  QR  Ldtc2>0  THEN 

1297  PRINT  USING  “16X. "“Thermocouples  were  defective  at  locat ions: “ * ,2(3D, 4X) “ ; 

Ldtcl ,Ldtc2 

1300  END  IF 

1303  ENTER  *File2;Itt 

1306  END  IF 

1309!  IF  Im-0  AND  Ihm-1  THEN  1595 
1312  BEEP 

1315  INPUT  “WANT  TO  CREATE  A  PLOT  FILE?  ( 0-N. 1 -Y) * , Iplot 
1318  IF  Ip lot* 1  THEN 
1321  BEEP 

1324  INPUT  “GIVE  NAME  FOR  PLOT  FILE* ,P_f i leS 
1327  CREATE  80AT  P  file*. 4 
1330  ASSIGN  SPlot  TO  P_file$ 

1333  END  IF 
1336  IF  Ihm-1  THEN 
1339  BEEP 

1342  INPUT  “WANT  TO  CREATE  Uo  FILE?  < 0-N, 1 -Y ) * . Iuf 
1345  IF  Iuf - 1  THEN 
1348  BEEP 

1351  INPUT  “ENTER  Uo  FILE  NAME".UfiIeS 


99 


13S4 

CREATE  BOAT 

Uf ileS. 4 

1357 

ASSIGN  SUfile  TO  UfileS 

1360 

END  IF 

1363 

BEEP 

1366 

INPUT  "UANT 

TO  CREATE  Re  FILE?  (0-N. 1 -Y) * . Ire 

1369 

IF  Ire-1  THEN 

1372 

BEEP 

1375 

INPUT  "ENTER 

Re  FILE 

NAME", Ref ileS 

1378 

CREATE  BDAT 

RefileS. 

10 

1381 

ASSIGN  defile  TO  RefileS 

1384 

END  IF 

1387 

END  IF 

1390 

PRINTER  IS  1 

1393 

IF  Ia-0  THEN 

1396 

BEEP 

1399 

PRINT  USING 

"4X, ""Select  tube  number*"" 

1402 

IF  Ihm-0  THEN 

1405 

PRINT  USING 

"SX. *"0 

Smooth  4  inch  Ref**" 

1408 

PRINT  USING 

*SX.  *"1 

Smooth  4  inch  Cu  (Preas/Slide) 

141  1 

PRINT  USING 

"SX. *"2 

Soft  Solder  4  inch  Cu""* 

1414 

PRINT  USING 

"SX.*"3 

Soft  Solder  4  inch  HIGH  FLUX"" 

1417 

PRINT  USING 

"6X.""4 

Uieland  Hard  8  inch*"* 

1420 

PRINT  USING 

"6X.""5 

HIGH  FLUX  3  inch""* 

1423 

PRINT  USING 

"6X.""6 

GEUA-K  19  Fina/in""" 

1426 

PRINT  USING 

"6X.""7 

GEMA-K  26  Fina/ln"** 

1429 

PRINT  USING 

*6X. "*8 

GEWA-T  19  Fma/in"* 

1432 

PRINT  USING 

"6X."*9 

GEHA-T  26  Fma/in""" 

1435 

PRINT  USING 

"6X.*"10 

THERMOEXCEL-E"*" 

1438 

PRINT  USING 

"6X,*"11 

THERMOEXCEL-HE*"* 

1441 

ELSE 

1444 

PRINT  USING 

"SX, *"0 

Smooth  tube*"* 

1447 

PRINT  USING 

"6X.""1 

High  Flux""* 

1450 

PRINT  USING 

"SX. ""2 

Turbo-B"*" 

1453 

PRINT  USING 

"6X.""3 

High  Flux  Mod"*" 

1456 

PRINT  USING 

"6X, **4 

Turbo-B  Mod*** 

1459  END  IF 

1462  INPUT  Itt 

1465  OUTPUT  OF i 1*1; Itt 

1468  END  IF 

1471  PRINTER  IS  701 

1474  IF  Itt<10  THEN  PRINT  USING  *16X,‘"Tube  Number:  *",D":Itt 

1477  IF  Itt>9  THEN  PRINT  USING  "16X,""Tube  Number:  **,DD*:Itt 

1480  IF  Ihm-1  THEN  PRINT  USING  "16X,"*Tube  Type:  — , ISA’ : TnS( Itt) 

1 483  BEEP 

1486  INPUT  "ENTER  OUTPUT  VERSION  < O-LONG. I -SHORT, 2-NONE) " , Iov 
1 489  BEEP 

1492  INPUT  "SELECT  (0-LIQ, 1 -VAP,2-(LIQ+VAP)/2) " . Ilqv 
1495! 

1498!  DIMENSIONS  FO  TESTED  TUBES 

1S01 !  ELECTRIC  HEATED  MODE 

1504!  D1-0lameter  at  thermocouple  positions 

1507  OATA  .01 11 125.. 0111 125, .0111125. .0129540.. 012446.. 0129540. .0100965 
1510  DATA  .0100965. .01  157. .01157, .01157. .01157 
1513  READ  Dial*) 

1516  DT -Dial Itt) 

1519! 

1522!  D2”0laaeter  of  test  section  to  the  base  of  fins 

1525  DATA  .015875.  .015875. .015875. .015824. .015875. .01S824. .01270 

1528  DATA  .0127, .0138,  .0138. .0138. .0138 

1531  READ  02a ( *) 


100 


i  .TLk  ".A 


1534! 

1537!  Di'tnside  diameter  of  unennanced  ends 

1540  D«TA  .0127, .0127. .0127, .0132. .0127. .0132. .01 11 125,. 01 11 12S. .01 18. .01 18. .01 
18. .01  Hi 

1543  READ  Qia( * ) 

1 S46 ! 

1549!  Du“Outside  diameter  of  unenhanced  ends 

1S52  DATA  .015875. .015875, .015875. .015824. .01 5875.. 015824.. 01270 . .01270 , . 01 331  . 
. 01331  .  .01331  .. 01331 
1555  READ  Doa<*) 

1558! 

1SS1 !  L'Length  of  enhanced  surface 

1564  DATA  .1016.. 1016,. 1016,. 10 t 6,. 2032.. 2032.. 2032,. 2032.. 2032,. 2032.. 2032,. 20 
32 

1567  READ  La(  *) 

1570! 

1573!  Lu-Length  of  unenhanced  surface  at  the  ends 

1576  DATA  .Q2S4, .02S4, .0254. .0254, .0762. .0762. .0762, .0762. .0762, .0762. .0762. .07 
62 

1579  READ  Lua ( * ) 

1582! 

1585!  Kcu-Thermal  Conductivity  of  tube 

1588  DATA  398, 344, 344. 4S, 344, 45. 344, 944, 398, 398, 398. 398 

1591  READ  Kcua( * ) 

1594  IF  Ihm-1  THEN 
1597! 

1600!  Data  statements  for  water  heating  mode 
1603! 

1606  DATA  0.015875.0.015875.0.0169.0.0138.0.0169.0.0.0,0.0 
1609  READ  D2a( *) 

1612  DATA  0.0127.0.0127.0.0145,0.0127,0.0145.0.0.0.0.0 
1615  READ  Dial*) 

1613  DATA  0. 015875, 0.01S87S.0. 0169, 0. 015875. Q. 0169, 0.0, 0,0,0 
1621  READ  Doa ( * ) 

1624  DATA  0.3048.0.3048,0.3048.0.3048,0.3048.0.0,0.0.0 
1627  READ  La(*> 

1630  DATA  0.0254.0.0254,0.0254.0.0254,0.0254,0.0.0.0.0 
1633  READ  lua(*) 

1636  DATA  398. 4S.398, 45,398.0 ,0 ,0 ,0 .0 
1633  READ  Kcua<*) 

1642  END  IF 
1645  D2-02a(Itt) 

1648  Di-Oia(Itt) 

1651  Do-Ooa(Itt) 

1654  L-ta(Itt) 

1657  Lu-Lua(Itt) 

1660  Kcu-Kcua( Itt) 

1663  Xn-.8 

1666  Fr-.3 

1669  IF  Itt-0  THEN  Cf-t.70E+9 

1672  IF  ltt>0  THEN  CF-3.7037EH0 

1675  A-PI*<Do‘2-Di‘2)/4 

1673  P'PI*0o 

1681  IF  Ihm-I  THEN 
1 684  BEE P 

1687  INPUT  ’TUBE  INITIATION  NODE.  (I-HOT  WATER. 2-STEAM. 3-COLD  UATER)’.Itim 
1690  IF  Itim-1  THEN  PRINT  USING  "16X.*"Tube  Initiate:  Hot  Water1*"” 

1693  IF  Itim-2  THEN  PRINT  USING  "16X,""Tube  Initiate:  Steam""" 

1696  IF  Itim-3  THEN  PRINT  USING  "16X.”"Tube  Initiate:  Cold  Water""’ 

1699  INPUT  "TEMP/VEL  MODE:  (0-T-CONST.V-DEC: ! -T-DEC, V-C0NST;  2-T-INC. '/-CONST) * . 

I  tv 


1702  IF  Itv-0  THEM  PRINT  USING  * 16X. H ‘Temp/Vel  Mode:  Constant/Decreasing 

1705  IF  Itv-1  THEN  PRINT  USING  *1SX.  "Temo/Vel  Mode:  Decreasing/Constant 

1708  IF  I  tv-2  THEN  PRINT  USING  ‘ 16X. " “Temp/Ve 1  Mode:  Increasing/Constant 

1711  INPUT  “WANT  TO  RUN  WILSON  PLOT?  ( 1 -Y.0-N) *,IuuI 

1714  IF  Ihn-1  AND  Iuil-0  THEN 

1717  IF  Itt-0  THEN  Ci-.032 

1720  IF  Itt-I  OR  Itt-3  THEM  Ci-.OSS 

1723  IF  Itt-2  OR  Itt-4  THEN  Ci-.062 

1 726  BEEP 

1729  INPUT  'ENTER  Cl  (DEF:  WH-.032,HF-.059,TB-.062>*,Ci 
1732  PRINT  USING  '16X.  "Sieder-Tate  *" 

173S  PRINT  USING  '16X,"  Constant  -  *'.Z.4D*:Ci 
1 738  ENO  IF 

1741  END  IF 

1744  IF  Ihn-1  AND  Im-1  AND  Ini 1-1  THEN 
1747  IF  Itt-0  THEN  Ci-,032 
1750  IF  Itt-1  OR  Itt-3  THEN  Ci-.0S9 
1753  IF  Itt-2  OR  Itt-4  THEN  Ci-.062 
1756  ASSIGN  @Fiie2  TO  * 

1759  CALL  Wilson(Cf.Ci) 

1762  ASSIGN  8Me2  TO  D2_fileS 

1765  ENTER  3File2:Nrun,Dold$.Ldtc1 ,Ldtc2. Itt 

1763  END  IF 
1771  Nsub-0 

1774  IF  I dp-4  THEN  Ihn-1 
1777  IF  Ihn-1  THEN  Nsub-8 
1780  Ntc-G 

1783  IF  Ihn-0  THEN  Ntc-12 

1786  J-1 

1789  Sx-0 

1792  Sy-0 

1795  Sxs-0 

1798  Sxy-0 

1801  Repeat:  ! 

1804  IF  In-0  THEN 
1807  Dtld-2.22 

1810  Ido-2 

1313  ON  KEY  0.15  RECOVER  1801 
1816  PRINTER  IS  1 

1819  PRINT  USING  '4X,  "SELECT  OPTION*" 

1822  PRINT  USING  ‘6X.  "0-TAKE  DATA’" 

132S  IF  Ihn-0  THEN  PRINT  USING  *SX."T-SET  HEAT  FLUX*" 

1828  IF  Ihn-1  THEN  PRINT  USING  *6X. "1-SET  WATER  FLOW  RATE"* 

1831  PRINT  USING  '6X.  "2-SET  Tsat"* 

1834  PRINT  USING  '4X.**N0TE:  KEY  0  -  ESCAPE”* 

1837  BEEP 

1840  INPUT  Ido 

1843  IF  Ido >2  THEN  Ido-2 

1846  IF  Ido-0  THEN  2296 

1849! 

1852!  LOOP  TO  SET  HEAT  FLUX  OR  FLOWMETER  SETTING 

1855  IF  Ido-1  THEN 

1358  IF  Ihn-0  THEM 

1861  OUTPUT  709; 'AR  AF12  AL13  VR5" 

1864  BEEP 

1867  INPUT  'ENTER  DESIRED  Odp'.Dodp 

1870  PRINT  USING  '4X.  "DESIRED  Qdp  ACTUAL  Qdp"' 

1873  Err-1000 

187S  FOR  1-1  TO  2 

1879  OUTPUT  709: 'AS  SA' 


1882  Su«-0 

1885  FOR  Ji-t  TO  5 

1888  ENTER  709 :£ 

1891  Su«-Suir*£ 

1894  NEXT  Ji 

1897  IF  1-1  THEN  Volt-Sun*5 
1900  IF  1-2  THEN  Amp-£ 

1903  NEXT  I 

190S  Aqdp-Ool  t*A*p/(PI*02*L) 

1909  IF  ABSi Aqdp-Dqdp) >£rr  THEN 

1912  IF  AqdP>Dqdp  THEN 

1915  BEEP  4000.. 2 

1918  BEEP  4000.. 2 

1921  BEEP  4000.. 2 

1924  ELSE 

1927  BEEP  2S0..2 

1930  BEEP  250.. 2 

1933  BEEP  2S0..2 

1936  END  IF 

1939  PRINT  USING  *4X.MZ.3DE,2X,MZ.3DE" :Dqdp,Aqdp 

1942  WAIT  2 

1945  GOTO  1876 

1948  ELSE 

1951  BEEP 

1954  PRINT  USING  *4X,MZ.3DE.2X,NZ.3DE';Dqdp,Aqdp 

1957  Err-500 

1960  WAIT  2 

1963  GOTO  1876 

1966  END  IF 

1969  ELSE 

1972  BEEP 

1975  INPUT  'ENTER  FLOWMETER  SETTING' ,F«a 

1978  GOTO  1819 

1981  END  IF 

1984  END  IF 

1987! 

1990!  LOOP  TO  SET  Tsat 
1993  IF  Ido-2  THEN 
1996  IF  lied t—  1  THEN  2011 
1999  BEEP 

2002  INPUT  'ENTER  DESIRED  Tsat'.Dtid 

2005!  PRINT  USING  *4X.**  DTsat  ATsat  Rate  Tv  Rate"' 

2008  Ilcdt-1 
2011  Qldl -0 
2014  01d2-0 

2017  Nn-1 
2020  Nrs-Nn  MOD  IS 
2023  Nn-Nn+1 
2026  IF  Nrs-t  THEN 

2029  PRINT  USING  *4X."  Tsat  Tidl  TId2  Tv  Tsump  Tinlet  Tpile 


Tout"* 

2032  END  IF 

2035  IF  Ih»-Q  THEN  OUTPUT  709:*AR  AF8  AL M  VR5 " 

2038  IF  Ihm-1  THEN  OUTPUT  709: *AR  AFO  ALS  VR5 " 

2Q41  FOR  I-l  TO  6 

2044  IF  Ih«-0  AND  I>4  THEN  2125 

2047  Sun-0 

2050  OUTPUT  709: 'AS  SA' 

2053  FOR  Ji-1  TO  20 
2056  ENTER  709:£Uq 
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2059  Sum-Sum*El  iq 
2062  NEXT  Ji 
2065  Eliq-Sum/20 
2068  Tld-FNTvsvlEl iq) 

2071  IF  1-1  THEN  TIdl-TId 

2074  IF  1-2  THEN  Tld2-TId 

2077  IF  1-3  THEN  Tv-Tld 

2080  IF  1-4  THEN  Tsump-Tld 

2083  IF  1-5  THEN  Tiniet-TId 

2086  IF  1-6  THEN  Tout- T Id 

2089  NEXT  I 

2092  IF  Iha- I  THEN 

209S  OUTPUT  709: *AR  AF20  AL20  VRS‘ 

2098  OUTPUT  709:* AS  SA“ 

2101  Su»-0 

2104  FOR  Kic-1  TO  20 

2107  ENTER  709:£ 

21 1 0  Sum-Sum*E 

2113  NEXT  Kit 

2116  Eat ( 7 ) -ABS ( Sun/20 ) 

2119  To i ie-Eat ( 7) /3.96E-4 
2122  END  IF 

2125  Atld-(Tld!>TId2)*.S 

2128  IF  ABSl Atld-Dtld) > .2  THEN 

2131  IF  Atld>DtId  THEN 

2134  BEEP  4000.. 2 

2137  BEEP  4000.. 2 

2140  8EEP  4000.. 2 

2143  ELSE 

2146  BEEP  250.. 2 

2149  BEEP  250.. 2 

2152  BEEP  250.. 2 

2155  END  IF 

2158  Errl-Atld-Qldl 

2161  Oldl-Atld 

2164  Err2*Tv-Qld2 

2167  QId2-Tv 

2170  IF  Iha-0  THEN  PRINT  USING  "4X .StMDD . DD.2X ) * :Dt Id. TIdl . Tld2 . Tv. Tsump 

2173  IF  Iha-1  AND  Idp-0  THEN  PRINT  USING  ’4X. 7 ( HDD .DD.2X) " ;Dt id. TIdl . TId2. Tv . Ts 

ump.Tirrlet.Tpile 

21 76*  IF  Iha-1  AND  Idp-4  THEN  PRINT  USING  "4X .5 ( HDD. DD.2X ) ,3(N3D .DD ,2X) * ;Dtld, TI 

d 1 , T I  d2 . T v , T sump .Tin 

2179  WAIT  2 

2182  GOTO  2020 

2185  ELSE 

2188  IF  ABS( At ld-Qt Id) > . 1  THEN 

2191  IF  Atld>0tld  THEN 

2194  BEEP  3000.. 2 

2197  BEEP  3000.. 2 

2200  ELSE 

2203  BEEP  800.. 2 

2206  BEEP  800.. 2 

2209  END  IF 

2212  Err  1 -Atld-QIdl 

2215  Oldl-Atld 

2218  Err2-Tv-01d2 

2221  Qld2-Tv 

2224  IF  Iha-0  THEN  PRINT  USING  *4X.5(MDD.DD.2X) “ ;Dt ld.Tldl ,Tld2. Tv. Tsump 
2227*  IF  Iha-1  THEN  PRINT  USING  *4X .5 ( HDD . DD . 2X ) ,3( M3D . DD . 1 X ) * ; Dt Id. T Idl . Tld2 . Tv 
. Tsump. T inlet. Tptle, 


2230 
2233 
223S 
2239 
2242 
2245 
2248 
225! 
22S4 
2257 
,  Tp  i  li 
2260 
2253 
2266 
2269 
2272 
2275! 
2278 
2281 
2284 
2287 
2290! 
2293! 
2296 
2299 
2302 
2305 
2308 
2311 
2314 
2317 
2320 
2323 
2326 
2329 
2332 
2335 
2338 
2341 
2344 
234  7 
2350 
2353 
2356 
2359 
2362 
2365 
2368 
2371 
2374 
2377 
2380 
2383 
2386 
2389 
2392 
2395 
2398 
2401 
2404 
2407 


up  T  T 

GOTO  2020 

ELSE 

8EEP 

Errl -Atld-Qldl 
Oldt "At  Id 
Err2-Tv-01d2 
01d2-Tv 

IF  Ih»-0  THEN  PRINT  USING  *4X . 5 (HDD . DD ,2X) ‘ : D t Id, Tld 1 , Tld2 . Tv , Tsump 
IF  Ihm-1  THEN  PRINT  USING  "4X ,8(MDD .QD.2X) * ;Dt Id, Tldl , Tld2. Tv, Tsump, T inlet 
i.  Tout 
WAIT  2 
GOTO  2020 
END  IF 
END  IF 
END  IF 

ERROR  TRAP  FOR  Ido  OUT  OF  BOUNDS 

IF  Ido>2  THEN 

BEEP 

GOTO  1819 
END  IF 

TAKE  DATA  IF  Im-0  LOOP 
IF  Ikol-1  THEN  2308 
BEEP 

INPUT  ’ENTER  BULK  OIL  Z*,Bop 
Ilcol*1 

IF  Ihm-0  THEN  OUTPUT  709: *AR  AFO  ALII  VRS* 

IF  Ihm-1  THEN  OUTPUT  709; ‘AR  AFO  AL5  VRS“ 

IF  Ihm-0  THEN  Ntc-12 
FOR  1-1  TO  Ntc 
OUTPUT  709: “AS  SA" 

Sum-0 

FOR  Ji-1  TO  20 
ENTER  709 :E 
Sum-Sum+£ 

IF  I-(9-Nsub)  OR  I-( 1 O-Nsub)  THEN  Etc Ji-t ) -E 

NEXT  Ji 

Kdl-0 

IF  I-(9-Nsub)  OR  I-( 1 O-Nsub)  THEN 

Eave-Sum/20 

Sum-0 . 

FOR  Jk-0  TO  19 

IF  ABSCEtUk  )-Eave)<S.QE-S  THEN 
Sum-Sum+Et( Jk ) 

ELSE 

Kdl-Kdl*1 
END  IF 
NEXT  Jk 

IF  I-( 9-Nsub f  OR  I-( 1 0-Nsub)  THEN  PRINT  USING  *4X.-*Kdi  -  **,DD*;Kdi 

IF  <d  I  > 1 0  THEN 

BEEP 

BEEP 

PRINT  USING  “4X,“’Too  much  scattering  m  data  -  repeat  data  set“““ 

GOTO  1816 
END  IF 
END  IF 

Erat ( I ) -Sum/ ( 20-Kdi ) 

NEXT  I 

IF  Ihm-1  THEN 

OUTPUT  709: 'AR  AF20  AL20  VRS" 


rnmm 


2590 
2593 
2596 
2599 
2602 
2605 
2608 
261  1 
2614 
2617 
2620 
2623 
2626 

2629 

2632 
2635 
2638 
2641 
2644 
2647 
2650 
2653 
2656 
2659 
2662 
2665! 
2668! 
2671 
2674 
2677 

2630 

2633 
2686 
2689 
2692 
2695 
2698 
2701 
2704 
2707 
2710! 
2713! 
2716 
2719 
2722 
2725 
2728 
2731 
2734 
2737! 
2740! 
2743 
E-10) 
2746! 
2749 
2752 
2755 
2758 
276! 
2764 


NEXT  I 

Tw«Twa/(8-Idtc) 

ENO  IF 

Tld*T(9-Nsub) 

Tld2"T ( 1 O-Nsub ) 

Tlda-<  Tld*Tld2)-*.S 
Tw«T ( 1 1 -Nsub ) 

IF  Itt<3  AND  Iha-0  THEN 
TId2— 99.99 
Tv-( T( 10) +T( 1 1 ) )/2 
END  IF 

Tsuap-T ( 12-Nsub) 

IF  Iha-0  THEN  2635 
Tinlet**T( 13-Nsub ) 

Tout-T ( 1 4-Nsub ) 

IF  Iha-0  THEN 
Aap-ABS( Ir ) 

Voit-ABS(Vr)*2S 
Q-Volt*flap 
END  IF 

IF  Itt-0  AND  Iha-0  THEN 
Kcu-FNKcu(  Tw ) 

ELSE 

Kcu-Kcua( Itt) 

END  IF 

FOURIER  CONDUCTION  EQUATION  WITH  CONTACT  RESISTANCE  NEGLECTED 
IF  Iha-0  THEN  T«-T*-Q*(_QG(  D2/D1 )/(2*PI*Kcu*L) 

IF  Ilgw-0  THEN  Tsat-Tlda 
IF  Ilqw-1  THEN  Taat-< Tlda+Tw)*.5 
IF  Ilqw-2  THEN  Tsat-Tw 
IF  Iha-1  THEN 
Tawg-Tinlet 

Grad-37 .9853+-.  1 04388*Tawo 
Tdrop-ABSlEat  ( 7 1  1*1  .  E-*-6/  ( 1  0*Grad) 

Tavgc-Txnlet-Tdrop*.5 
IF  A8S( Tawg-Tawgc) >.0!  THEN 
Tawg-<  Tawg*Tawgc)*.5 
GOTO  2689 
END  IF 

COMPUTE  HATER  PROPERTIES 
IF  Iha-1  THEN 
Kw-FNKu( Tawg ) 

Mu**a*FNMuw<  Tawg) 

Cpw*FNCpw( Tawg) 

P  rw*FNP  rvi  <  Tawg ) 

Rho*-FNRho*(  Tawg) 

Twx-Tawg 

Compute  MOOT 

Mdot-3.9657E-3>Fms*<3.6!955E-3-Fras*<8.82006E-6-Fms*< 1 .23688E-7-Fms*4.31897 

)  ) 

Ndot»Hdot*< 1  . 0365-Txnlet'M  1 . 96644E-3-T iniet+5 .252E-S) ) /1 .0037 
Kdt-0 

Q-Hdo  t *Cpa*  Td  rop 

Lratd-Tdrop/LQGt  (Tmlet-Tsat)/(T tnle t-Tdrop-Tsat )  ) 

Uo-Q/(PI*Do*L*Lmtd) 

Rw*0o*LQG( Do/Di  1712. *Kcu) 

T«*  Tsat+Fr*Lmtd 


27S7  V*-Hdot/(Rhow-PI*Oi ‘2/4) 

2770  Raw*Rhov*Vv*Q i/Muua 

2773  H  t*Ci*Kia/D  l -Raw ' . 8*Prw  * ( 1 /3 . ) *( Muwa/FNMuw( Twi ) ) * . 1 4 
2776  Twic*Tavq-Q/(PI*Oo*L*Hi ) 

2779  IF  PBS ( Twi -Twtc) > . 0 1  THEN 

2782  Twi*( Twi+T*ic)*.S 

2785  GOTO  2773 

2788  END  IF 

273  T  Tm-<  Twi  +  Tvic)*.S 

2794  Ho-1/< 1/Uo-Do/(Dt*Hi)-Rw) 

2797  END  IF 

2800  END  IF 

2803  IF  Iha-1  THEN 

2806  Thatab*Q/(Ho*PI*Oo*L) 

2809  Tw»Tsat*Thatab 

2812  ELSE 

2815  Th«tab“T«-Tsat 

2813  END  IF 

2821  IF  Th«tab<0  THEN 
2824  BEEP 

2827  INPUT  *TUALL<TSAT  (O-CONTINUE.  1 -END) ‘.lav 

2830  IF  Iav-0  THEN  GOTO  1804 

2833  IF  lav- I  THEN  3130 

2836  END  IF 

2839! 

2842!  COMPUTE  VARIOUS  PROPERTIES 
284S  Tf ila-( T*+Tsat)*.5 
2848  Rho-FNRholTf ilnJ 
2851  Mu-FNHu ( Tf i la ) 

2854  K-FNK(Tf iln) 

2857  Cp-FNCp(Tfilii) 

2860  Bata-FNBeta< Tf i la) 

2863  Hfq-FNHfg< Tsat) 

2866  Ni-Mu/Rho 
2869  AIpha-K/fRho-Cp) 

2872  Pr- Ni/Alpha 
2875  Psat-FNPsat< Tsat) 

2878! 

2881!  COMPUTE  NATURAL -CONVECTIVE  HEAT-TRANSFER  COEFFICIENT 
2884!  FOR  UNENHANCED  END(S) 

2887  Hbar-190 

2890  Fa-(Hbar-*P/(Kcu-A) ) " .5*Lu 
2893  Tanh-FNTanh(Fa) 

2896  Thata-Thatab*Tanh/Fa 

2899  Xx-(9.81 *8ata*Thatab*0o '3*Tanh/(F9-*Ni*Alpha) ) * . 166667 

2902  Yy-d-M  .559/Pr)  *  (9/16)  )  *  (8/27 ) 

2905  Hbarc«K/Do*< .6*.387*Xx/Yy ) '2 

2908  IF  ABSdHbar-Hbarc) /Hbarc)  >  .  001  THEN 

2911  Hbar-(Hbar-Hbarc)*.5 

2914  GOTO  2890 

2917  END  IF 

2920! 

2923!  COMPUTE  HEAT  LOSS  RATE  THROUGH  UNENHANCED  ENO(S) 

2926  Ql-(Hbar*P*Kcu*A) * .5*Th»tab*Tanh 
2929  Qc-Q-2*Q1 
2932  Aa-PI*02*L 
2935! 

2938!  COMPUTE  ACTUAL  HEAT  FLUX  AND  BOILING  COEFFICIENT 

2941  Qdp-Qc/As 

2944  Htuba-Qdp/Thetab 
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2947  Cst‘«<Cp*Thetab/Hfg>/<Qdp/(Mu*Hfg>*<  . 0 1 4/ (  9. 31  *Rho )  *  .5)  ‘  <  1  /3. ) *Pr  ‘  1 .  7 ) 

2950! 

2953!  RECORD  TIME  OF  DATA  TAKING 
29S6  IF  In-0  THEN 
2959  OUTPUT  709: 'TD* 

2962  ENTER  709: ToldS 

2965  END  IF 

2968! 

2971 !  OUTPUT  DATA  TO  PRINTER 
2974  PRINTER  IS  701 
2977  IF  Iov-0  THEN 
2980  PRINT 

2983  PRINT  USING  *10X.**Data  Sat  Nuaber  ■  **  .DDD.2X,  "Bulk  Oil  %  -  ".DD.O.SX.l 
4A* : J.Bop. ToldS 
2986  IF  Iha-0  THEN 

2989  PRINT  USING  *10X,"TC  No:  1  2  3  4  5  6  7 

0-  ■  ■ 

2992  PRINT  USING  *10X."Teap  :"  ,8( 1X.MDD.DD) * ;T(l),T(2),T(3),T(4),Tl5).T(6),T< 
7) .1(8) 

2995  PRINT  USING  *10X."  Twa  Tliqd  TIigd2  Tvapr  Psat  Tsuao*" 

2998  PRINT  USING  *1 0X.21MD0 .00. 1 X) , 1 X .MOD. DD . IX, 2< 1 X.NDD.DD) .2X.MDD. D* ; Tw, Tld , T 
ld2, Tv. Psat, Tsuao 

3001  PRINT  USING  *10X."  Thetab  Htube  Qdo*" 

3004  PRINT  USING  ‘1  OX, MOD. 3D. 1 X.MZ.3DE. 1X.MZ.3DE* ; The tab, Htube. Qdp 
3007  ELSE 

3010  PRINT  USING  *10X,"  Fa*  V*  Tsat  Tinl  Tdrop  Thatab  q  Uo 

Ho"* 

3013*  PRINT  USING  *1  OX, 412D.DD. IX ) . Z.3D , JX.DD.DD. 1X.3(MZ.3DE. IX) * ;F«s, Vw.Tsat , Ti 

niet , Tdrop. Thatab.Qd 

3016  END  IF 

3019  END  IF 

3022  IF  Iov-1  THEN 

3025  IF  J-1  THEN 

3028  PRINT 

3031  IF  Iha-0  THEN 

3034  PRINT  USING  *10X."  RUN  No  0iI7.  Tsat  Htuba  Qdp  Thatab*** 

3037  ELSE 

3040  PRINT  USING  "10X."  FMS  OILX  TSAT  HTUBE  QDP  THETAB*** 

3043  END  IF 

3046  END  IF 

3049  IF  Ihm-0  THEN 

3052  PRINT  USING  *12X, 3D. 4X,D0.2X. MOD. DD,3( IX, MZ.3DE) *:J. Bop, Tsat, Htube, Qdp. The 
tab 

3055  ELSE 

3053  PRINT  USING  "12X.3D.4X.DD.2X ,M0D.DD,3( IX.MZ.3DE) * :Fn»s,8op, Tsat, Htube, Qdp, T 
he  tab 

306 1  END  IF 
3064  END  IF 
3067  IF  Ia-0  THEN 
3070  BEEP 

3073  INPUT  "OK  TO  STORE  THIS  DATA  SET  (1-Y,0-N)?*.0k 
3076  END  IF 

3079  IF  Ok  —  1  OR  Im-t  THEN  J-J-t-1 
3082  IF  Ok  - 1  ANO  Ia-0  THEN 

3085  IF  Iha-0  THEN  OUTPUT  W i la  I  :  Boo, ToldS, Eraf ( * ) . Vr. Ir 
3088  IF  Iha-1  THEN  OUTPUT  SFtlel ;8op.ToidS.Emf (*> ,Fms 

3091  END  IF 

3094  IF  luf-l  THEN  OUTPUT  »Ufile:Vw.Uo 

3097  IF  Ire-1  THEN  OUTPUT  §Ref iie:Fms.Re« 

3100  IF  (Ira-!  OR  Ok-I)  ANO  Iplot-1  THEN  OUTPUT  *Plot:Qdp, Thetab 
3103  IF  Ia-0  THEN 


■i* 

'el 

*2 

V. 

'»* 


I 


4$ 

i 


3106  BEEP 

3109  INPUT  “HILL  THERE  BE  ANOTHER  RUN  < 1-Y.Q-N)?* ,Go_on 
3112  Nrun*J 

3115  IF  Go_on-0  THEN  3130 
3118  IF  Go  onOO  THEN  Repeat 
3121  ELSE 

3124  IF  J<Nrun+1  THEN  Repeat 
3127  END  IF 
3130  IF  Ib-0  THEN 
3133  BEEP 

3136  PRINT  USING  *10X.**NOTE:  ".ZZ.**  data  rune  Mere  stored  in  file  **,10A*;J 
1 ,D2_file* 

3139  ASSIGN  JFilel  TO  * 

3142  OUTPUT  #File2:Nrun-1 
3145  ASSIGN  SFilel  TO  Dt  file* 

3148  ENTER  9File1 :DateS.Ldtc1 ,Ldtc2. Itt 
3151  OUTPUT  9Fiie2:Date$,Ldtc1 ,Ldtc2.Itt 
3154  FOR  1-1  TO  Nrun-1 
3157  IF  Ih«-0  THEN 

3160  ENTER  SFilel :Bop, ToldS.Emf (♦) .Vr.Ir 
3163  OUTPUT  «File2:Bop.ToldS,Eaf(*) , Vr.Ir 
3166  ELSE 

3169  ENTER  SFilel :Bop. ToldS.Eaf (*) .Frw 
3172  OUTPUT  #File2;Bop,ToldS,Eaf (•) .Fms 
3175  END  IF 
31 78  NEXT  I 

3181  ASSIGN  «F;le1  TO  * 

3184  PURGE  "DUMMY" 

3187  END  IF 
3190  BEEP 
3193  PRINT 

3196-  IF  Iplot-1  THEN  PRINT  USING  " 1 0X . ""NOTE:  **,ZZ.”  X-Y  pairs  uere  stored  i 
plot  data  f i le  ,  1 
3199  ASSIGN  3File2  TO  * 

3202  ASSIGN  aPlot  TO  * 

3205  IF  Inf— 1  THEN  ASSIGN  »Ufile  TO  * 

3208  IF  Ire-1  THEN  ASSIGN  SRefile  TO  * 

3211  CALL  Stats 
3214  BEEP 

3217  INPUT  'LIKE  TO  PLOT  DATA  (1 -Y ,0-N) ,0k 
3220  IF  Olr-1  THEN  CALL  Plot 
3223  SU8END 
3226! 

3229!  CURVE  FITS  OF  PROPERTY  FUNCTIONS 
3232  DEF  FNKcu(T) 

3235!  OFHC  COPPER  250  TO  300  K 

3238  Tk-T+273. 15  !C  TO  K 

3241  <-434-.  11 2*  Tic 

3244  RETURN  < 

3247  FNEND 
3250  DEF  FNMu(T) 

3253!  170  TO  360  <  CURVE  FIT  OF  VISCOUS ITY 

3256  Tk-T+273. 15  !C  TO  K 

3259  Mu-EXP ( -4. 4636+1  1 01 1  .47/Tlc )  )*1 .0E-3 

3262  RETURN  Mu 

3265  FNEND 

3268  DEF  FNCp(T) 

3271!  180  TO  400  <  CURVE  FIT  OF  Cp 

3274  Tk-T+273. 15  !C  TO  < 

3277  Cp-. 40188+1 .65007E-3*Tk+! .5!494E-S*Tk ,2-6.67853E-10*Tk *3 


3280  -  Cp*Cp*1 000 
3233  RETURN  Co 
3236  FNENO 
3289  DEF  FNRho(T) 

3292  Tlc-T+273. 15  !C  TO  K 

3295  X—  t  —  <  t  .8*Tlc/753.95)  !K  TO  R 

3298  Ro-36.32+61 . 146414*X* ( 1 /3)+16.41801S*X+1 7.476838*X' .5+1 . 1 19828*X*2 

3301  Ro-Ro/. 062423 

3304  RETURN  Ro 

3307  FNENO 

3310  DEF  FNPr(T) 

3313  Pr-FNCp(T)*FNHu(T)/FNK(T) 

3316  RETURN  Pr 
3319  FNENO 
3322  DEF  FNK(T) 

3325!  T<360  K  WITH  T  IN  C 
3328  K* . 071 - . 000261 +T 
3331  RETURN  X 
3334  FNENO 
3337  DEF  FNTanh(X) 

3340  P-EXP(X) 

3343  Q-1/P 

3346  Tanh-(P-0)/(P+Q) 

3349  RETURN  Tanh 
3352  FNEND 
3355  DEF  FNTvsvIV) 

3358  COM  /Cc/  C(7),Ical 
3361  T-C(O) 

3364  FOR  1-1  TO  7 
3367  T-T+C(I)*V‘I 
3370  NEXT  I 
3373  IF  Ical-1  THEN 

3376  T-T-6. 7422934E-2+T*(  9. 0277043£-3-T+( -8.325991 7E-5) > 

3379  ELSE 

3382  T-T+8. 626897E-2+T*  <  3 . 76 1 99E-3- T*5 . 0689259E-5 ) 

3385  END  IF 
3388  RETURN  T 
3391  FNEND 
3394  DEF  FNBeta(T) 

3397  Rop-FNRho( T+. 1 ) 

3400  Ro«-FNRho(  T  - .  1) 

3403  Be ta--2/ ( Rop+Roe ) * ( Rop-Ron ) / . 2 

3406  RETURN  Beta 

3409  FNEND 

3412  DEF  FNHfglT) 

3415  Hfg-1 . 374 1344E+5~T*<  3. 309436 1E+2+T+1 .2165143) 

3418  RETURN  Hfg 

3421  FNEND 

3424  DEF  FNPsat(Tc) 

3427!  0  TO  80  deg  F  CURVE  FIT  OF  Pgat 
3430  Tf-1.8*Tc+32 

3433  Pa-5.94552S+Tf+( . 1S3S2082+Tf+< 1 .4840963E-3+Tf+9.S150671E-6) ) 

3436  Pg-Pa-14.7 

3439  IF  Pg>0  THEN  !  +-PSIG.—  m  Ha 

3442  Psat-Pg 
3445  ELSE 

3448  Psat-Pg*29.92/14.7 

3451  END  IF 
3454  RETURN  Psat 
3457  FNEND 


3460  DEF  FNHs»ooth<X.Boo. I sat) 

3463  DIM  6(5) .BIS) ,C(S) .0(5) 

3466  DATA  .20526. .2S322. .319048. .55322. .79909. 1 .00258 
3463  DATA  .74515. . 72992. .731 89, . 71 22S . .68472 . .641 97 
3472  DATA  .41092. . 17726. .25142. .54806. .31916. 1 .0845 
3475  DATA  .71403. . 72913. . 72S6S. .696691  . .665867 . .61 889 
3478  READ  A(*).B(*),C(*).D(*) 

3481  IF  Bop<6  THEN  I-Bop 
3484  IF  Bop-6  THEN  1-4 
3487  IF  Bop- 10  THEN  1-5 
3490  IF  Isat- I  THEN 
3493  Hm-EXP<A(I)+B(I)*LOG(X)) 

3496  ELSE 

3499  Hm-EXP<C(I)*0<I)*L0G<X) ) 

3502  ENO  IF 
350S  RETURN  H* 

3508  FNEND 

3511  DEF  FNPoly(X) 

3514  COM  /Cply/  A(  10, 10)  ,C(  10)  .B(4 )  .Nop ,  Ipmt.Qpo, Ilog, If n, I jom.Njoin 
3517  Xl-X 
3S20  Poly-B(0) 

3S23  FOR  1-1  TO  Nop 

3526  IF  Ilog-1  THEN  XI-LOG(X) 

3529  Poiy-Poly+8( I)*X1‘ I 
3532  NEXT  I 

3535  IF  Ilog-1  THEN  Poly-EXP (Poly) 

3538  RETURN  Poly 
3541  FNEND 
3544  SUB  Poly 

3547  DIM  R(  10)  ,S(  10) .Sy( 12) .Sx( t2) .Xx( 100) .Yy( 100) 

3550  COM  /Cply/  A( 10.10) ,C( 10) ,8(4) .N.Iprnt , Opo . Ilog.Ifn.Ijoin.Njoin 
3S53  COM  /Xxyy/  Xp(  25 ) .  Yp  ( 25 ) 

3556  FOR  1-0  TO  4 
3559  B( I) -0 
3562  NEXT  I 
3565  BEEP 

3568  INPUT  ‘SELECT  (0-FILE. 1 -KEYBOARD. 2-PROGRAM) ‘.Im 
3571  Ia-Im+1 
3574  BEEP 

3577  INPUT  ‘ENTER  NUMBER  OF  X-Y  PAIRS*. No 
3580  IF  Im- 1  THEN 
3583  BEEP 

3586  INPUT  ‘ENTER  DATA  FILE  NAME* ,D  film* 

3589  BEEP 

3532  INPUT  ‘LIKE  TO  EXCLUDE  DATA  PAIRS  <  t-Y,0-N)?Mmd 
3535  IF  Imd-1  THEN 
3598  BEEP 

3601  INPUT  ‘ENTER  NUMBER  OF  PAIRS  TO  BE  EXCLUDED* , Ipmx 
3604  END  IF 

3607  ASSIGN  3Film  TO  D  film* 

3610  ELSE 
3613  BEEP 

3616  INPUT  ‘WANT  TO  CREATE  A  DATA  FILE  ( 1 -Y.0-N)?* . Yms 
3619  IF  Yma-1  THEN 
3622  BEEP 

3625  INPUT  'GIVE  A  NAME  FOR  DATA  FILE*,D_f lie* 

3628  CREATE  BDAT  D  file*. 5 
3631  ASSIGN  *Film  TO  0  film* 

3634  ENO  IF 
3637  END  IF 


crow 


v  \  •  wvvvwv; 


iwwswww 


3640 

3643 

3646 

3649 

3652 
3655 

3653 
3661 
3664 
3667 
3670 
3673 
3676 
3679 
3682 
3685 
3688 
3631 
3694 
3697 
3700 
3703 
3706 
3709 
3712 
37  IS 
3718 
3721 
3724 
3727 
3730 
3733 
3736 
3739 
3742 
3745 
3748 
3751 
3754 
3757 
3760 
3763 
3766 
3769 
3772 
3775 
3778 
3781 
3784 
3787 
3790 
3793 
3796 
3799 
3802 
3805 
3808 
3811 
3814 
3817 
3320 


THEN 


BEEP 

INPUT  “ENTER  THE  ORDER  OF  POLYNOMIAL* ,N 
FOR  1*0  TO  N*2 
Sy l I ) *0 
Sx(  D-0 
NEXT  I 

IF  led* I  AND  In-1 
FOR  1*1  TO  Ipex 
ENTER  *Fil«;X.Y 
NEXT  I 
END  IF 

FOR  1*1  TO  Np 
IF  Ip-1  THEN 

Opo-2  THEN  ENTER  9FiU:X.Y 
0po<2  THEN  ENTER  9FiU:Y.X 
Opo-i  then  r-r/x 
Ilog-1  THEN 
Opo-2  THEN  Xt-X/Y 


THEN 

THEN 


Y-LOG(Xt) 
Y-L0G1 Y) 


THEN 


IF 
IF 
IF 
IF 
IF 

X-LOG(X) 

IF  Opo-2 
IF  0po<2 
END  IF 
END  IF 
IF  Ip-2 
BEEP 

INPUT  “ENTER  NEXT  X-Y  PAIR’.X.Y 
IF  Yes-1  THEN  OUTPUT  9File;X.Y 
END  IF 

IF  Ip<3  THEN 
Xx(I)-X 
Yy  (  I ) - Y 
ELSE 

X-Xp(I-l) 

Y-Yp(I-l) 

END  IF 
R(0)-Y 

Sy(0)-Sy(0)>Y 

S(D-X 

Sx(i)-Sxtn+x 
FOR  J-1  TO  N 
R( J)*R( J-1 )*X 
Sy ( J)-Sy ( J) +R( J) 

NEXT  J 

FOR  J-2  TO  N*2 
SU)-S(J-1  )*X 
Sx( J)-Sx( J)+S( J) 

NEXT  J 
NEXT  I 

IF  Yes-1  ANO  la-2  THEN 
BEEP 

PRINT  USING  *12X,QD. ""  X-Y  pairs 

END  IF 

Sx(0)-Np 

FOR  1-0  TO  N 

CCI)-Sy(I) 

FOR  J-0  TO  N 
A( I. J)-Sx< I*J) 

NEXT  J 
NEXT  I 

FOR  I-fl  TO  N-1 
CALL  Dividel I) 


'xere  stored  in 
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3823  CALL  Subtraction 

382S  NEXT  I 

3829  B(N)-C(N)/A(N.N) 

3832  FOR  I- 0  TO  N-l 
3835  B( N-1 -I ) »C ( N-1  -I ) 

3838  FOR  J-0  TO  I 

3841  8 1 N- 1  - 1 )  -B  ( N- 1  - 1 )  -A  ( N- 1  - 1 ,N- J ) *B ( N- J ) 

3844  NEXT  J 

3847  8  <  N- 1  - 1 )  -B  ( N- 1  - 1 )  /A  <  N- 1  - 1 ,  N- 1  - 1 ) 

3850  NEXT  I 
3853  '.PRINTER  IS  701 
3856  IPRINT  B<*) 

3859  ! PRINTER  IS  70S 
3862  IF  Iprnt-0  THEN 

3865  PRINT  USING  *12X,* ’EXPONENT  COEFFICIENT**’ 

3868  FOR  1-0  TO  N 

3871  PRINT  USING  " 1 5X , DO . SX , MO . 7DE " ;  I ,  B  ( I ) 

3874  NEXT  I 
3877  PRINT  *  * 

3880  PRINT  USING  *12X.”DATA  POINT  X  Y  Y ( CALCULI 

EPANCY”* 

3883  FOR  1-1  TO  Np 
3886  Yc-B(O) 

3889  FOR  J-1  TO  N 
3892  Yc-Yc+B( J)*Xx(  I) * J 
3895  NEXT  J 
3898  D-Yy< I) -Yc 

3901  PRINT  USING  *1SX,3D.4X.4(MD.5DE, IX) ’ ; I ,Xx( I) . Yy ( I) , Yc.O 

3904  NEXT  I 

3907  END  IF 

3910  ASSIGN  SFile  TO  * 

3913  SUBEND 
3916  SUB  DividslM) 

3919  COM  /CpLy/  A<  10. 10 )  ,CC  1 0)  ,B(4)  ,N.Iprnt,Opo.  Ilog,  If n.I jom.Njoin 
3922  FOR  I-M  TO  N 
3925  Ao-A(I.M) 

3928  FOR  J-M  TO  N 

3931  A< I.J)-A( I, J)/Ao 

3934  NEXT  J 

3937  C(I)-C(I)/Ao 

3940  NEXT  I 

3943  SUBENO 

3946  SUB  Subtract(K) 

3949  COM  /Cply/  A( 10, 10) ,C( IQ1 ,S(4) ,N,Iprnt,Opo.Ilog,Ifn.I Join.Njoin 

3952  FOR  I-K  TO  N 

3955  FOR  J-K- 1  TO  N 

3958  A(I.J)«A(K-1 ,J)-Al I.J) 

396 1  NEXT  J 

3964  Ct  D-CCK-l  )-CC  I) 

3967  NEXT  I 
3970  SU8END 
3973  SUB  PUn 

3976  COM  /Cply/  A< 1 0 . 1 0 ) . C( 1 0 ) .81 4) ,N, Iprnt ,0po . Ilog, If n. I Join.N join 
3979  COM  /Xxyy /  Xx<25) . Yy ( 25 ) 

3982  PRINTER  IS  705 
3985  BEEP 

3988  INPUT  "WANT  TO  PLOT  Uo  vs  Vw7  (1 -Y , 0-N ) * , Iuo 
3991  IF  Iuo-0  THEN 
3994  BEEP 

3997  INPUT  ’SELECT  (0-h/h07.  gams  tubs.  1  -MHF) /Ms*)  *  , Irt 


Y( CALCULATED)  DISCR 


4000  BEEP 

4003  INPUT  "SELECT  h/h  RATIO  ( 1 -FILE. 0 -COMPUTED) " . Ihrat 
4006  IF  Ihrat-0  THEN 
4009  BEEP 

4012  INPUT  "WHICH  Tsat  ( 1-6. 7,0— 2.2) * , Isat 

4015  END  IF 

4018  Xain-0 

4021  Xaax-10 

4024  Xstep-2 

4027  IF  Irt-G  THEN 

4030  Yain-0 

4033  Yaax- 1.4 

4036  Ystep-. 2 

4039  ELSE 

4042  Ymin-0 

4045  Yaax- 15 

4048  Ystep-5 

4051  END  IF 

4054  ELSE 

4057  Opo-2 

4060  Yain-0 

4063  Yaax- 12 

4066  Ystep-3 

4069  Xain-0 

4072  Xnax-4 

4075  Xstep- 1 

4078  END  IF 

4081  IF  Ihrat- 1  THEN 

4084  Yam-0 

4087  Yaax-18 

4090  Ystep-3 

4093  Xain-0 

4096  Xaax-9 

4099  Xstep-2 

4102  END  IF 

4105  BEEP 

4108  PRINT  "INlSPT : IF  2300.2200,3300,6800;* 

4111  PRINT  "SC  0 , 1 00 . 0 , 1 00 ; TL  2.0;“ 

4114  Sf  x- 1 00/(Xnax-Xam) 

4117  Sfy-1 00/( Yaax-Yain) 

4120  PRINT  "P D  0.0  PD" 

4123  FOR  Xa-Xain  TO  Xaax  STEP  Xstep 

4126  X-(Xa-X»in)*Sfx 

4129  PRINT  *PA":X,*,0;  XT;* 

4132  NEXT  Xa 

4135  PRINT  "PA  100.0;PU;" 

4138  PRINT  *PU  PA  0.0  PD* 

4141  FOR  Ya-Ynun  TO  Yaax  STEP  Ystep 
4144  Y-( Ya-Yrain) *Sfy 

4147  PRINT  *PA  0.*;Y,*YT* 

4150  NEXT  Ya 

4153  PRINT  "PA  0.100  TL  0  2* 

4156  FOR  Xa-Xain  TO  Xaax  STEP  Xstep 

4159  X-(  Xa-Xn»in)*Sf  x 

4162  PRINT  "PA* ;X, * , 100:  XT" 

4165  NEXT  Xa 

4168  PRINT  “PA  100.700  PU  PA  100.0  PD" 

4171  FOR  Ya-Yain  TO  Yaax  STEP  Ystep 

4174  Y-( Ya-Yain)*Sfy 

4177  PRINT  *PD  PA  100.".Y,*YT* 

4180  NEXT  Ya 


4183  PRINT  'PA  100. 100  PU* 

4136  PRINT  'PA  0,-2  SR  1 .5.2“ 

4189  FOR  Xa-Xain  TO  Xaax  STEP  Xstep 
4192  X-<Xa-Xainl*Sfx 

4144  PRINT  “Pa"*y  *  fl • * 

4198  IF  Iuo-0  THEN  PRINT  *CP  -2.-1 ;LB* :Xa! ** 

4201  IF  Iuo-1  THEN  PRINT  'CP  -1 .S. -1 ;LB" ;Xa; ** 

4204  NEXT  Xa 

4207  PRINT  'PU  PA  0.0“ 

4210  FOR  Ya-Yam  TO  Yaax  STEP  Ystep 
4213  IF  A8S(Ya)<1.£-5  THEN  Ya*0 
4216  Y-(  Ya-Yatn)*Sfy 

4219  PRINT  “PA  0.“;Y.““ 

4222  IF  Iuo-0  THEN  PRINT  “CP  -4. -.2S;LB* : Ya: ** 

4225  IF  Iuo-1  THEN  PRINT  “CP  -3,-.2S;LB*;Ya:*“ 

4228  NEXT  Ya 

4231  XlabalS-’Oll  Percent* 

4234  IF  Iuo-0  THEN 

4237  IF  Irt-0  THEN 

4240  YlabelS“*h/HQZ* 

4243  ELSE 

4246  YlabelS-'H/hsaooth' 

4249  END  IF 

4252  PRINT  “SR  1.S.2:PU  PA  50,-10  CP“:-LEN(XlabelS)/2:“0:L8“:XlabeiS:*“ 

42S5  PRINT  “PA  -11,50  CP  0 , * ; -LEN( Y labels) /2*5/6:  ‘DI  0, 1 ;LB* ;  YiabelS; ** 

42S8  PRINT  ‘CP  0.0“ 

4261  ELSE 

4264  PRINT  “SP0:SP2* 

4267  PRINT  “SR  1.2.2.4JPU  PA  -8,3S:DI  0.1;LBU:PR  I,0.S:LBo:PR  -1.0.5;LB  (IcW/m 

4270  PRINT  'PR  - 1 , 0 . 5 : SR  1.1.5:LB2:SR  1.5.2:PR  .S..5:LB.:PR  .S.O’.LBK)' 

4273  PRINT  'PA  42.-10:01  1,0;LBV;PR  .4.-1;LB«:PR  1 , .S;LB(m/s) ' 

4276  •  PRINT  “SPO:SP1“ 

4279  END  IF 
4282  Ipn-0 
4285  BEEP 

4288  INPUT  'UANT  TO  PLOT  OATA  FROM  A  FILE  (  1  •  Y .  0-N)  ?“  ,Qlrp 

4291  Icn-0 

4294  IF  Olrp-1  THEN 

4297  BEEP 

4300  INPUT  'ENTER  THE  NAME  OF  THE  DATA  FILE'.D  file* 

4303!  IF  Iuo-0  THEN 
4306  BEEP 

4309  INPUT  'SELECT  (0 -LINEAR.  1 -L0G(X. Y) “ . Hog 
4312!  END  IF 

4315  ASSIGN  SFile  TO  0  files 
4318  BEEP 

4321  INPUT  “ENTER  THE  BEGINNING  RUN  NUMBER “,Md 
4324  BEEP 

4327  INPUT  “ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED* .Npairs 
4330  IF  Iuo-0  ANO  Ihrat-0  THEN 
4333  BEEP 

4336  INPUT  “ENTER  DESIRED  HEAT  FLUX" .Q 

4339  END  IF 

4342  BEEP 

4345  PRINTER  IS  1 

4348  PRINT  USING  '4X, ““Select  a  symbol:*** 

4351  PRINT  USING  '4X,'“1  Star  2  Plus  Sign'” 

4354  PRINT  USING  '4X.'*3  Circle  4  Square”' 

4357  PRINT  USING  '4X.”S  Rombus”" 


I 
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4360  PRINT  USING  “4X.*“S  Rignt-side-up  triangle'1** 
4363  PRINT  USING  *4X,**7  Up-3ide-down  triangle"** 
4366  INPUT  Sya 
4363  PRINTER  IS  70S 
4372  PRINT  “PU  01* 

4375  IF  Sya-1  THEN  PRINT  *SM** 

4378  IF  Sy*-2  THEN  PRINT  *SM+“ 

4381  IF  Sya-3  THEN  PRINT  *SMo" 

4384  Nn-4 

4387  IF  Ilog" !  THEN  Nn- I 
4390  IF  Md>1  THEN 
4393  FOR  1-1  TO  (ttd-t) 

4396  ENTER  3File:Xa.Ya 
4399  NEXT  I 
4402  END  IF 
440S  IF  Ihrat-0  THEN 
4408  QI-Q 

4411  IF  Ilog-I  THEN  Q-LOG(O) 

4414  END  IF 
4417  FOR  1*1  TO  Npairs 
4420  IF  Iuo-0  AND  Ihrat-0  THEN 
4423  ENTER  3File;Xa.B<*> 

4426  Ya-8(0) 

4429  FOR  K-1  TO  Nn 
4432  Ya-Ya+8(K)*Q"K 
4435  NEXT  < 

4438  END  IF 

4441  IF  Iuo-1  OR  Ihrat-1  THEN 
4444  ENTER  $File:Xa.Ya 
4447  IF  Iuo-1  THEN  Ya-Ya/1000 
4450  END  IF 

4453  IF  Iuo-0  AND  Ihrat-0  THEN 
4456  IF  Ilog-1  THEN  Ya-EXP<Yal 
4459  IF  Ilog-0  THEN  Ya-Ql/Ya 
4462  IF  Irt-0  THEN 
4465  IF  Xa-0  THEN 
4468  Yo-Ya 
4471  Ya-1 
4474  ELSE 
4477  Ya-Ya/Yo 
4480  END  IF 
4483  ELSE 

4486  H«a-FNHsmooth(Q,Xa. Isat) 

4489  Ya-Ya/Hsa 
4492  END  IF 
4495  END  IF 
4498  Xx( 1-1 )-Xa 
4501  Yy ( 1-1 ) -Ya 
4504  X-<  Xa-Xain) *Sf x 

4507  Y-(Ya-Yain)*Sfy 
4510  IF  Sya>3  THEN  PRINT  *SN* 

4513  IF  Sya<4  THEN  PRINT  *SR  1.4, 2. 4* 

4516  PRINT  *PA“ ,X. Y. "* 

4519  IF  Sya>3  THEN  PRINT  *SR  1 .2.1.5" 

4522  IF  Sym-4  THEN  PRINT  *002.4.99.0,-8.-4.0.0,3.4 
4525  IF  Sya-5  THEN  PRINT  *UC3.0 , 99 . -3, -6 . -3.6.3 .6. 
4528  IF  Sya-6  THEN  PRINT  *UC0 .5 .3 . 99 .3 . -8 . -6 . 0 .3 .3: * 
4S31  IF  Sy«-7  THEN  PRINT  *UC0 . -5.3 .99, -3.8. 6,0 . -3. -8 
4534  NEXT  I 
4537  BEEP 

4540  ASSIGN  3Fi le  TO  * 


to- 


4543  END  IF 
4546  PRINT  *PU  SM* 

4549  BEEP 

4552  INPUT  “WANT  TO  PLOT  A  POLYNOMIAL  ( I *Y , Q-N) 7  ,Okp 
4555  IF  Qkp-1  THEN 

4553  BEEP 

4561  PRINTER  IS  1 

4564  PRINT  USING  *4X. "Select  line  type:"* 

4567  PRINT  USING  *SX,"0  Solid  line"* 

4570  PRINT  USING  *6X."1  Dashed"* 

4573  PRINT  USING  *SX."2.,,S  Longer  line  -  dash"* 

4576  INPUT  Ipn 
4579  PRINTER  IS  705 
4582  BEEP 

4585  INPUT  'SELECT  <0-LlN£AR,t-L0G<X.Y))*.Ilog 

4S88  Iprnt-1 

4591  CALL  Poly 

4594  IF  Iuo-1  THEN 

4597  BEEP 

4600  INPUT  ’DESIRE  TO  SET  X  Lower  and  Upper  Limit  l 1 -Y ,0-N>7* , Ixl im 

4603  IF  Ixlim-0  THEN 

4606  XI 1-0 

4609  Xul-7 

4612  END  IF 

4615  IF  Ixlim-1  THEN 

4618  BEEP 

4621  INPUT  ’ENTER  X  Lower  Limit*, XI 1 
4624  BEEP 

4627  INPUT  'ENTER  X  Upper  Limit *,Xul 
4630  END  IF 
4633  END  IF 

4636  FOR  Xa-Xll  TO  Xul  STEP  Xstep/25 
4639  Icn-Icn-1 
4642  Ya-FNPoly(Xa) 

4645  IF  Iuo-1  THEN  Ya-Ya/1000 

4648  Y-( Ya-Ymin)*Sfy 

4651  X-(Xa-Xain)*Sf x 

4654  IF  Y<0  THEN  Y-0 

4657  IF  Y>1 00  THEN  GOTO  4637 

4660  Pu-0 

4663  IF  Ipn- I  THEN  Idf-Icn  MOD  2 

4666  IF  Ipn-2  THEN  Idf-Icn  MOD  4 

4669  IF  Ipn-3  THEN  Idf-Icn  MOD  8 

4672  IF  Ipn-4  THEN  Idf-Icn  MOD  16 

4675  IF  Ipn-5  THEN  Idf-Icn  MOD  32 

4678  IF  Idf-1  THEN  Pu-1 

4681  IF  Pu-0  THEN  PRINT  ’PA*.X, Y, *PD* 

4684  IF  Pu-1  THEN  PRINT  *PA* ,X, Y, *PU* 

4687  NEXT  Xa 
4690  PRINT  *PU* 

4693  GOTO  428S 
4696  END  IF 
4699  BEEP 

4702  INPUT  'WANT  TO  QUIT  ( 1 -Y ,0-N>7* , Iqui t 
4705  IF  Iquit-1  THEN  4711 
4708  GOTO  4285 
4711  PRINT  *PU  SPO’ 

4714  SUBEND 
4717  SUB  Stats 
4720  PRINTER  IS  701 
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4723 

J-0 

4728 

K-0 

4729 

BEEP  „  ,  ,  , 

4732 

IF  Ip  lot" I  THEN  ASSIGN  SFzle  TO  P.filaS 

473S 

BEEP 

4738 

INPUT  'LAST  RUN  No7 ( 0-QUIT) * ,Nn 

i 

4741 

IF  Nn-0  THEN  4849 

i 

4744 

Nn-Nn-J 

f 

4747 

S*-0 

1 

4750 

Sy-0 

4753 

St-0 

47S6 

Sxs-0 

4759 

Sys-0 

1 

4782 

Szs-0 

; 

4765 

FOR  1-1  TO  Nn 

4768 

J-J+1 

4771 

ENTER  3Fil®:Q.T 

4774 

H-Q/T 

4777 

Sx-Sx+Q 

4780 

Sxs-Sxs-K)*2 

4783 

Sy-Sy-T 

4786 

Sya-Sys+T*2 

i 

4789 

Sz-Sz+H 

4792 

Szs-Szs+H'2 

4795 

NEXT  I 

j 

4798 

Qav«-Sx/Nn 

I 

4801 

Tavq-Sy/Nn 

4804 

Havq-Sz/Nn 

% 

4807 

Sd«vq-SQR(A8S< (Nn*Sxs-Sx'2)/(Nn-*<Nn-1 ) ) ) ) 

- 

4810 

Sd«vt-SOR( ABS( (Nn*Sys~Sy‘2)/(Nn*(Nn-1 ) ) ) ) 

V 

• 

4813 

Sdevh-SOR(  A8S(  (Nn-*Szs-Sz‘2)  /(Nn*(Nn-1 )  ))  ) 

J 

4816 

Sh- 1 00*Sdevh/Have 

L* 

4819 

Sq- 1 00*Sdevq/Qave 

M 

4822 

St-1 00*Sdevt/Tave 

1 

4825 

IF,  K-l  THEN  4843 

1 

4828 

PRINT 

4831 

PRINT  USING  *1 IX,  "DATA  FILE:  ”  , I4A* ;FileS 

i 

4834 

PRINT 

SdevQ  Thetab  Sd«vT" 

1 

X 

4337 

m 

PRINT  USING  "1 IX.  "RUN  Htube  SdevH  Qdp 

'} 

4840 

K-1 

3D .  IX. 3D . 2D ' : J . Have , Sh , Qave 

U 

4843 

PRINT  USING  “1 1X,DD,2(2X,D.3DE. 1X.3D.2D) .2X.DD. 

8 

,Sq, Tave.St 

I 

4846 

GOTO  4735 

3 

4849 

ASSIGN  $File1  TO  * 

S 

4852 

PRINTER  IS  1 

r* 

4855 

SUBEND 

L 

4858 

SUB  Coet 

,  If  n, Ijoin.Njoin 

1 

4861 

4864 

4867 

CON  /Cp iy /  A(I0,10),CC10),B(4) .N. Iprnt.Opo. I log 
BEEP 

INPUT  'GIVE  A  NANE  FOR  CROSS-PLOT  FILE'.CpfS 

I 

4870 

CREATE  8DAT  CpfS.2 

V 

4873 

ASSIGN  SFile  TO  CptS 

& 

4876 

BEEP 

4879 

INPUT  'SELECT  ( 0 -LINEAR. 1 -L0G(X. Y) ) * , Ilog 

4882 

BEE? 

•v 

4885 

INPUT  'ENTER  OIL  PERCENT  (- 1 -STOP )*, Bop 

4888 

IF  Bop<0  THEN  4900 

s 

4891 

CALL  Poly 

1 

4894 

OUTPUT  4File:8op.B<*> 

p 

>' 

V 

4897 

GOTO  4882 

V 

* 

[ 

119 

4900 

4903 

4906 

4909 

4912 

4915! 

4918 

4921 

4924 

4927 

4930 

4933 

4936 

4939 

4942 

4945 

4948 

4951 

4954 

49S7 

4960 

4963 

4966 

4969 

4972 

4975 

4978 

4981 

4984 

4987 

4990 

4993 

4996 

4999 

5002 

5005 

5008 

501 1 

5014 

5017 

5020 

S023 

5026 

5029 

5032! 

5035 

5038 

5041 

5044 

5047 

5050 

5053 

5056 

5059 

5062 

5065 

5068 

5071  ! 

5074! 

5077 


ASSIGN  sFile  TO  * 

SUBEND 

SUB  Mi  IsoniCf .Ci ) 

CON  /Mil/  02. Di. Do. L.lu.Xcu 
DIM  Emt (12) 

WLISON  PLOT  SUBROUTINE  DETERMINE  CF  AND  CI 
BEEP 

INPUT  "ENTER  DATA  FILE  NAME". FiUS 
BEEP 

PRINTER  IS  1 

PRINT  USING  "4X, * “Select  option:'"* 

PRINT  USING  “4X . *'  0  Vary  Cf  and  Ci*"" 

PRINT  USING  "4X , 1  Fix  Cf  Vary  Ci*'* 

PRINT  USING  *4X."  2  Vary  Cf  Fix  Ci"' 

INPUT  'ENTER  OPTION' , let ix 

PRINTER  IS  701 

IF  Icfix-0  THEN  4960 

IF  Icf ix>0  THEN  BEEP 

IF  Icf ix- I  THEN  INPUT  'ENTER  Cf'.Cst 

IF  Icfix-2  THEN  INPUT  'ENTER  CI'.Ci 

PRINTER  IS  1 

INPUT  'Uant  To  Vary  Coet f P( 1 -V , Q-N ) * , Iccoef 
IF  Iccoef  —  1  THEN  INPUT  "ENTER  COEFF'.R 
PRINTER  IS  701 

IF  Icfix-0  OR  Icfix-2  THEN  Cfa-.004 

IF  Icf  ix- 1  THEN  Cfa-Csf 

Cia-Ci 

Xn-.3 

Fr-.3 

Jj-0 

Rr-3. 

IF  Iccoef- 1  THEN  Rr-R 
PRINTER  IS  1 
PRINT  Do.Dt.Kcu 
ASSIGN  4FUa  TO  File* 

ENTER  *Fi  iel.Nrun, Dates, Ldtcl ,Ldtc2,Itt 
Rw-0o*L0G ( Do/0 i ) / ( 2*Kcu ) 

Sx-0 

Sy-0 

Sxy-0 

Sx2-0 

Sy2-0 

FOR  1-1  TO  Nrun 

ENTER  4F  i la: Bop. To  ids. Emt ( * ) ,  Fms 
CONVERT  EMF'S  TO  TEMPERATURE 
FOR  J-1  TO  5 
T(  J)-FNTv3v(Emt ( J)  ) 

NEXT  J 

Tsat-(T( 1 )+T(2) )*.S 
Tavg-T (5) 

Grad-37.3853*. !04388*Tavg 
Tdrop-tait  (  7 )  *1  .  E+6/  ( 1  0  .  *6 rad) 

T  avgc- T ( 5 ) - Td  rop* . 5 
IF  A8S l Tavg-Tavgc ) > . 01  THEN 
Tavg-<  Tavg  +  favgci *.S 
GOTO  5050 
END  IF 


Compute  properties  of  water 
Km-FNKw ( Tavg ) 


5080  Muua*FNMuw*  Tavg ) 

5083  Cpw*FNCpu i Tavg ) 

5086  Pru*FNPru( Tavg ) 

5089  Rho«”FNRhow( Tavg) 

5092! 

5095!  Compute  properties  of  Freon-114 

5098  Lmtd” Tdrop/LQG( (T(S)-Tsat)/(T(S)-Tdrop-Tsat)  ) 

5 ! 0  T  IF  Jj-0  THEN 
5104  Tu“Tsat+Fr*Lmtd 
S107  Tbetab*Tu-Tsat 
5110  Jj-1 
5113  ENO  IF 
5116  Tf * < Tu+Tsa t ) * . S 
5119  Rho-FNRho(Tf ) 

5122  Mu-FNMu(Tf) 

5125  K-FNK(Tf) 

5128  Cp-FNCo(Tf) 

5131  Beta-FNBeta( Tf ) 

5134  Hf g-FNHf g( Tsat ) 

5137  Ni“Mu/Rho 
5140  Alpha-K/(Rho*Cp) 

5143  Pr-Ni/Alpha 
5146! 

5149!  Analysis 
5152!  COMPUTE  MOOT 
5155  A“PI*( Do '2*0 i ‘2) /4 
5158  P-PI*0o 

5161  Mdo  t«  3 . 9657E- 3  +Fms*  < 3.61 9S5E-3-Fms* ( 8 . 82006E-6-F ms* ( 1 .23688E-7-Fms*4.31897 
E- 10))) 

5164  Q“Mdot*Cpw*Tdrop 

5167!  COMPUTE  NATURAL -CONVECTIVE  HEAT-TRANSFER  COEFFICIENT 
5170!  FOR  UNENHANCED  END ( 3 ) 

5173  Hbar-190 

5176  Fe*(Hbar*P/(Kcu*A) ) * .5*Lu 
5179  Tanh-FNTanh(Fe) 

5182  Theta-Thetab*Tanh/Fe 

5185  Xx”(9.81*8eta*Thetab*Oo'3*Tanh/ (Fe*Ni ♦Alpha)  )  * .  166667 
5188  Yy«< 1 ♦( .559/Pr) ‘ (9/16)) • (8/27) 

5191  Hbarc-</0o*< .6+.387*Xx/Yy ) *2 

5194  IF  ABS ( ( Hbar-Hbarc) /Hbar ) > . 001  THEN 

5197  Hbar* ( Hbar^Hba re ) * . 5 

5200  GOTO  5176 

5203  END  IF 

5206! 

5209!  COMPUTE  HEAT  LOSS  RATE  THROUGH  UNENHANCED  ENDS 

5212  01 *( Hbar*P*Kcu*A ) * . 5*Thetab*Tanh 

5215  Qc*0-2*Q1 

5218  As-PI*02*L 

5221 !  COMPUTE  ACTUAL  HEAT  FLUX 

5224  QdP”Qc/As 

5227  IF  let ix*0  OR  Ictix>1  THEN  Csf « 1 /Cf * ( 1 . /Rr ) 

5230  Thetao-Csf /Ce*Hfg*(Qdp/(Mu*Hfg)*( . 0 1 4/ ( 9 . 81 *Rho ) ) * .S) ‘ ( 1 /Rr ) *Pr *  1 . 7 

5233  Ho“Qdp/Thetab 

5236  Omega *Ho/Cf 

5239  Uo-0/(PI*0o*L4Lmtd) 

5242  Vw-Mdot/(Rhou*PI*0i '2/4) 

5245  Re«”Rhow*Vw*0 i /Muua 
5248  Twi*Tw*Q*Rg/(PI*0o*L) 

5251  Gan»a*<u/D  i *Reu ' . 3*Pru ' ( 1/3. )*( Muua/FNMuwl Tw  i )  )  ‘ .  1 4 

5254!  PRINTER  IS  1 

5257  Yu” < 1 . /Uo-Ru ) *0mega 
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5260  Xe-Oeega*Oo/ < Gama*0 1 ) 

5263  Sx-Sx*Xe 

5266  Sy-Sy+Ye 

5269  Sxy-Sxy*Ye*Xe 

5272  Sx2-Sx2+Xe*Xw  * 

5275  Sy2"Sy2+,Yw*Yu 

5278  NEXT  I 

5281  ASSIGN  SFile  TO  * 

5284  M-<  Sx*Sy-Nrun*Sxy ) / (Sx*Sx-Nrun*Sx2) 

5287  C-(Sy-Sx*M)/Nrun 

5290  IF  Ictix-0  OR  Icfix-3  OR  Icftx-4  THEN 

5293  Cic-1/N 

5296  Cfc-1/C 

5299  END  IF 

5302  IF  Icfix-1  THEN 

5305  Cic-1/N 

5308  Cfc-Cf 

5311  END  IF 

5314  IF  Icfix-2  THEN 

5317  Cic-Ci 

5320  Cfc-1/C 

5323  END  IF 

5326  IF  A0S< (Ci-Cic)/Cic)>.001  OR  ABS< <Cf-Cfc)/Cfc)>.001  THEN 
5329  Ci-(Ci*Cic)*.S 
5332  Cf-(Cf-K:fc)*.5 
5335  PRINTER  IS  1 

5338  PRINT  USING  *2X.”  Csf  -  ”  .NZ.3DE.2X.  ”  Ci  -  ”,NZ.3DE*:Csf  ,Ci 

5341  PRINTER  IS  701 

5344  GOTO  5002 

5347  END  IF 

5350  PRINT 

5353  PRINTER  IS  701 

5356  PRINT  USING  "23X.”  Cf  Ci*” 

5359  PRINT  USING  *3X "ASSUMED  ”  ,MZ:3DE.3X.MZ.3DE" :Cfa.Cia 

5362  PRINT  USING  *8X. "'CALCULATED  ” .MZ.3DE.3X.MZ.3DE* ;Caf .Ci 

5365  PRINT 

5368  Su*2-Sy2-2*M*Sxv-2'K>$y*M*2*Sx2*2*M*C*Sx+Nrun*C’2 
5371  PRINT  USING  "10X.”Su«  of  Squares  -  ** .Z.3DE" :Sun2 
5374  PRINT  USING  "10X. "Coefficient  -  "* .D.DDD" :Rr 
5377  SUBEND 
5380  DEF  FNMue(T) 

5383  A-247 .3/ ( T+1 33. 1 5 ) 

5386  Mu-2.4E-5*10*A 

5389  RETURN  Mu 
5392  FNENO 
5395  DEF  FNCpe(T) 

5398  Coe-4.21 1208S8-T*(2.26826E-3-T*< 4.42361E-5+2. 71428E-7*T) ) 

5401  RETURN  Coe* 1000 

5404  FNEND 

5407  DEF  FNRhoe(T) 

5410  Ro-999 . 52946+ T* ( .01269-T*<S.482513E-3-T*1 .234147E-5) ) 

5413  RETURN  Ro 
5416  FNEND 
5419  DEF  FNPre(T) 

5422  P rw-FNCoe ( T ) *FNMue ( T ) /FNKe  <  T ) 

5425  RETURN  Pre 

5423  FNEND 

5431  DEF  FNKe(T) 

5434  X-(T*273. 151/273. 15 

5437  Ke--.92247+X*< 2.9395-X*(  1 ,3007-X*( . 52577- . 07344*X ) ) ) 


RETURN  Kw 
FNEND 
SUB  Plot 

COM  /Cply/  A(t0,10),C(!0).B(4) ,Nop . Iprnt ,Qpo .Ilog.Ifn.Ijoin.Njoin 
DIM  Bs<3) 

INTEGER  Ii 
PRINTER  IS  1 
BEEP 

INPUT  "SELECT  HER TING  MODE  (O-ELECTRIC,  f -WATER) ", Iha 

Idv-Q 

8EEP 

INPUT  "LIKE  DEFAULT  VALUES  FOR  PLOT  (1  -Y,0-N)?",ldv 

Opo-O 

BEEP 

PRINT  USING  *4X. "'Select  Option:*"" 

PRINT  USING  "SX.""0  q  varsus  delta-T""" 

PRINT  USING  *SX.""1  h  varsus  dalta-T""" 

PRINT  USING  *SX , *"2  h  varsus  q"*" 

INPUT  Opo 
BEEP 

INPUT  "SELECT  UNITS  (0-SI, t -ENGLISH) ",  Ian 
PRINTER  IS  705 
IF  Idv<M  THEN 
BEEP 

INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  X-AXIS", Cx 
BEEP 

INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  Y-AXIS", Cy 
BEEP 

INPUT  "ENTER  MIN  X-VALUE  (MULTIPLE  OF  tO>".X*in 
BEEP 

INPUT  "ENTER  MIN  Y-VALUE  (MULTIPLE  OF  tO>*,Yniin 
ELSE 

IF  0po-0  THEN 

IF  Iha-1  THEN 

Cy-t 

Cx-2 

Xnun- 1 

Ymin-1 0000 

ELSE 

Cy-3 

Cx-3 

X«in« . 1 

Yrain- ! 00 

END  IF 

END  IF 

IF  Opo- 1  THEN 

IF  Ih»-I  THEN 

Cy-2 

Cx-2 

Xain- . 1 

Yn»in-I00 

ELSE 

Cy-3 

Cx-2 

Xnun- 1 000 
Ymin-1 00 
END  IF 
END  IF 

IF  Qpo-2  THEM 
IF  Ih«-1  THEN 
Cy-2 


Cx*1 

Xmm-iOOOO 
Yain*1 00 
ELSE 
Cv-3 
Cx-3 
Xain-.l 
Yarn- 100 
END  IF 
END  IF 
END  IF 
BEEP 

PRINT  *IN:SP1;IP  2300. 2200. 8300. 6800:* 

PRINT  ‘SC  0,100.0. 100’.  TL  2.0:“ 

Sfx- 100/Cx 
Sfy-100/Cy 

INPUT  "WANT  to  BY-PASS  CAGE?  ( t-Y.O-N) Ibyp 
IF  Ibyp* 1  THEN  SO 49 
PRINT  *PU  0.0  PD' 

Nn-9 

FOR  I- I  TO  Cx-1 
Xat-Xain*10*  1 1-1 ) 

IF  I -Cx-1  THEN  Nn-1 
FOR  J*1  TO  Nn 

IF  J-1  THEN  PRINT  *TL  2  0* 

IF  J-2  THEN  PRINT  ‘TL  1  O' 

Xa-Xat*J 

X-LGT ( Xa/Xa i  n ) *Sf  x 
PRINT  *PA*;X,',0:  XT;* 

NEXT  J 
NEXT  I 

PRINT  'PA  100.0;PU:* 

PRINT  *PU  PA  0.0  PD* 

Nn-9 

FOR  I-t  TO  Cy-H 
Yat-Yain*10*(I-1) 

IF  I-Cy+1  THEN  Nn-1 
FOR  J-1  TO  Nn 

IF  J-1  THEN  PRINT  *TL  2  0* 

IF  J-2  THEN  PRINT  'TL  1  O' 

Ya-Yat*J 

Y-LGT(Ya/Y*vn)*Sfy 
PRINT  "PA  0. * ; Y, *YT* 

NEXT  J 
NEXT  I 

PRINT  'PA  0,100  TL  0  2* 

Nn-9 

FOR  1-1  TO  Cx*1 
Xat-Xain*1 0*  ( 1-1 ) 

IF  I-Cx-H  THEN  Nn-1 

FOR  J-1  TO  Nn 

IF  J-1  THEN  PRINT  'TL  0  2* 

IF  J>1  THEN  PRINT  *TL  0  1* 

Xa-Xat*J 

X-LGT (Xa/Xa in )*Sfx 
PRINT  'PA';X,'.100:  XT' 

NEXT  J 
NEXT  I 

PRINT  'PA  100,100  PU  PA  100.0  PD' 
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5303  Nn-9 

5306  FOR  1-1  TO  Cy-1 
5809  Yat-Y*in*l0‘  (1-1 ) 

5812  IF  I-Cy+1  THEN  Nn-I 

5815  FOR  J*1  TO  Nn 

5318  IF  J-1  THEN  PRINT  *TL  0  2* 

5821  IF  J>1  THEN  PRINT  *TL  0  1* 

5824  Ta-Yat*J 

5827  Y-LGT  ( Ya/Ymn)*Sfy 

5830  PRINT  “PD  PA  10Q,*.Y,“YT* 

5833  NEXT  J 
5836  NEXT  I 

5839  PRINT  “PA  100.100  PU* 

5842  PRINT  “PA  0.-2  SR  1.5.2“ 

5845  Ii-LGT(X«in) 

5848  FOR  1-1  TO  Cx-M 
5851  Xa-X»in*1 0 ' C I - 1 ) 

5854  X-LGT(Xa/X»in)*Sf x 
5857  PRINT  “PA*:X.*,0:“ 

5860  IF  Ii>-0  THEN  PRINT  “CP  -2.  -2:LB10.‘PR  -2.2:LB*:IiJ“* 

5863  IF  li<0  THEN  PRINT  “CP  -2 . -2;LB1 0 ;PR  0,2;LB“;Ii:“* 

5866  Ii-Ii+1 

5869  NEXT  I 

5872  PRINT  “PU  PA  0.0* 

5875  Ii-LGT( Y«xn) 

5878  Y10-10 

5881  FOR  1-1  TO  Cy+1 

5884  Ya-Ymin*t  0  * ( 1-1 ) 

5887  Y-LGT(Ya/Yaun)*Sfy 
5890  PRINT  “PA  0.*;Y.“* 

5893  PRINT  “CP  -4,-.2S:LB10;PR  -2.2 :LB" Hi ; ““ 

5896  Ii-Ix-1 
5899  NEXT  I 
5902  BEEP 

5905  INPUT  “WANT  USE  DEFAULT  LABELS  ( 1  -Y  ,0-N)7* , Idl 
5908  IF  Idl<>1  THEN 
5911  BEEP 

5914  INPUT  “ENTER  X -LABEL" .X labels 
5917  BEEP 

5920  INPUT  “ENTER  Y-LABEL" , YlabeiS 

5923  END  IF 

5926  IF  Qpo<2  THEN 

5929  PRINT  “SR  1,2:PU  PA  40.-14:" 

5932  PRINT  "LBIT'.PR  -1.6,3  PD  PR  1.2.0  PUtPR  .5,-4;LB«o:PR  .5.1;* 

5935  PRINT  *LB-T;PR  .5.-1 ;LB3at:PR  .5,1;* 

5938  IF  Iun-Q  THEN 
S941  PRINT  “LB)  (K)“ 

5944  ELSE 

5947  PRINT  “LB)  (F)“ 

5950  END  IF 
5953  END  IF 
5956  IF  Oeo-2  THEN 
5959  IF  Iun-0  THEN 

5962  PRINT  “SR  J.S,2;PU  PA  40,-14;LBq  (W/m;SR  1.1.S;PR  0.5.i:LB2:SR  t.5.2;PR 
0.S.-1 ;LB) “ 

5965  ELSE 

5968  PRINT  “SR  !.5.2:PU  PA  34.-14;LBq  (Btu/hr:PR  .5, .5:LB. :PR  .5. -.5;“ 

5971  PRINT  “LBft-.PR  .5,1  ;SR  1,1.S:LB2;SR  1  .5,2:PR  .5.-1  ;LS ) ; " 

5974  END  IF 
5977  END  IF 
5980  IF  Qpo-0  THEM 
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5983  IF  Iun-0  THEN 

5986  PRINT  *SR  I . 5 . 2 : PU  PA  -12.40:01  O.tILBq  (W/m:PR  -I.O.S:SR  1.1.S:LB2:SR  I 
.S.2:PR  t. .5:L8) “ 

5989  ELSE 

5992  PRINT  “SR  1.5.2:PU  PA  -12.32:01  O.ULBq  (Btu/hr:PR  -.5. .5:LB. ;PR  .5. .5:* 
599S  PRINT  “LBftiSR  1.1.5;PR  -1,.S:LB2:PR  1,.5;SR  1,S.2:LB>* 

5998  END  IF 
6001  END  IF 
6004  IF  0po>0  THEN 
6007  IF  Iun-0  THEN 

6010!  PRINT  “SR  1.5.2;PU  PA  -12.38:DI  0,t:LBh  (W/m;PR  -1,.S;SR  1.1.S;LB2:SR  1. 
S • 2 • PR  • S •  5  * " 

6013  PRINT  “SR  1.2.2.4:PU  PA  -12.37:01  0,l:LBh;PR  1.0.5:L3o:PR  -1,0.S;LB  ( W/n 


6016  PRINT  “PR  -t,.S;SR  1,1.5;LB2:SR  1.5.2:PR  .S..S:LB.;PR  .S,0:LBK)“ 

6019  ELSE 

6022  PRINT  “SR  1.5.2:PU  PA  -12.28:01  0.1;LBh  <Btu/hr:PR  -.5.  .S:LB.  ,'PR  .5.. 5:' 

6025  PRINT  “LBftjPR  -1,.S:SR  1.1.5:LB2:SR  1.S.2:PR  .S,.5;LB.:PR  .5,.S:LBF)“ 

6028  END  IF 

6031  END  IF 

6034  IF  Idl-0  THEN 

6037  PRINT  “SR  1.5.2:PU  PA  50. -16  CP*:-LEN<XIabelS)/2:*0:LB*:XlabelS:** 

6040  PRINT  “PA  -14.50  CP  0 , “ ; -LEN(YLabelS) /2*5/6: “DI  0.1 ;LB*:YlabeiS:“" 

6043  PRINT  "CP  0,0  01* 

6046  END  IF 
6049  Ipn-0 
6052  X11-1.E+6 

6055  Xui  — 1.E+6 
60S8  Ien-0 
6061  Ifn-0 

6064  Ijom-1 
6067  BEEP 

6070  INPUT  “WANT  TO  PLOT  DATA  FROM  A  FILE  ( 1-r,0-N)?*,Ok 
6073  IF  Ok -1  THEN 
6076  BEEP 

6079  INPUT  “ENTER  THE  NAME  OF  THE  DATA  FILE’.D  files 
6082  ASSIGN  §File  TO  D_fileS 
6085  BEEP 

6088  BEEP 

6091  INPUT  “ENTER  THE  BEGINNING  RUN  NUMBER*. Md 
6094  BEEP 

6097  INPUT  “ENTER  THE  NUM8ER  OF  X-Y  PAIRS  STORED* .Npairs 
6100!  BEEP 

6103!  INPUT  “CONNECT  DATA  WITH  LINE  ( 1-Y,0-N)?*.Ici 

6106  BEEP 

6109  PRINTER  IS  1 

6112  PRINT  USING  “4X. ‘“Select  a  symbol:*** 

6115  PRINT  USING  *6X.**1  Star  2  Plus  sign*** 

6118  PRINT  USING  *SX.“*3  Circle  4  Square'** 

6121  PRINT  USING  ’6X.“*S  Rombus"** 

6124  PRINT  USING  “SX,**S  Right-side-up  triangle*** 

6127  PRINT  USING  *SX,“*7  Up-side-down  triangle*** 

6130  INPUT  Sym 

6133  PRINTER  IS  705 
6136  PRINT  “PU  DI" 

6139  IF  Sym- 1  THEN  PRINT  “SM*“ 

6142  IF  Sym-2  THEN  PRINT  *SM+* 

6145  IF  Sym-3  THEN  PRINT  “SMo" 

6148  IF  Md> 1  THEN 
61S1  FOR  1-1  TO  (Md-1) 


i 

il 


ENTER  -File: Ya.Xa 
NEXT  I 
END  IF 

FOR  1*1  TO  Noairs 

ENTER  3File:Ya.Xa 

IF  1-1  THEN  Ql-Ya 

IF  I*Npairs  THEN  Q2*Ya 

IF  0oo*1  THEN  Ya-Ya/Xa 

IF  Opo-2  THEN 

Q-Ya 

Ya-Ya/Xa 

Xa-Q 

END  IF 

IF  Xa<Xil  THEN  Xll-Xa 
IF  Xa>Xul  THEN  Xul-Xa 
IF  Iun-1  THEN 
IF  0po<2  THEN  Xa-Xa*1.8 
IF  0po>0  THEN  Ya-Ya*.17S1 
IF  Opo-O  THEN  Ya-Ya*.317. 

IF  Opo-2  THEN  Xa-Xa*.317 
END  IF 

X-LGT(Xa/X«in)*Sfx 

Y-LGT<Ya/Y»in)*Sfy 

Kj-0 

CALL  Syi»b(  X.  Y  ,Sy«,  Icl  ,K  j ) 

GOTO  S2S8 

IF  Sya^3  THEN  PRINT  *SN* 

IF  Sy«<4  THEN  PRINT  'SR  1.4, 2. 4* 

IF  Icl-0  THEN 
PRINT  “PA“ , X , Y , *  * 

ELSE 

PRINT  “PA* ,X, Y, “P0“ 

END  IF 

IF  Sv»>3  THEN  PRINT  *SR  1.2.1.S" 

IF  Syt»-4  THEN  PRINT  *UC2 . 4 . 39 . 0 .  -8 . -4 . 0 . 0 .8 . 4 . 0 :  * 

IF  Sy*-5  THEN  PRINT  “UC3. 0 .99 . -3 . -6 . -3 . 6 . 3 .5 .3 . -5: * 

IF  Sye-6  THEN  PRINT  *UC0 ,5.3,99 .3 , -8. -6 , 0 .3 .3; * 

IF  Sy«-7  THEN  PRINT  “UCO . -5 .3, 99. -3. 3 ,S, 0 , -3 . -8: * 

NEXT  I 
PRINT  *P1J* 

BEEP 

INPUT  *UANT  TO  LABEL?  ( 1-Y.O-N) *,Ilab 
IF  Ilab-1  THEN 
PRINT  “SP0:SP2" 

BEEP 

IF  < lab-0  THEN 
Xlao-5 
Y lab-95 

INPUT  “ENTER  INITIAL  X.Y  LOCATIONS* .Xlab. Ylab 

Xtt-Xlab-5 

Ytt-Ylab-8 

PRINT  *SR  1,1.5’ 

BEEP 

INPUT  “SELECT  LABEL  TYPE  (0-FILE  NAME.  1-TUBE  TYPEJ'.Ifntt 
IF  Ifntt-0  THEN 


6319 

PRINT  “3M;PA* , Xtt , Ytt, “LB 

7.  Heat 

File 

6322 

Ytt-Ytt-3 

6325 

PRINT  “3M:PA* ,Xtt , Ytt . “LB 

7.  Heat 

F  i  le 

6323 

ELSE 

6331 

PRINT  “3M:PA" .Xtt. Ytt. "LB 

7.  Hea  t 

T  ube 

Oil  Flux  Type 


6334  m-rtt-3 
6337  PRINT  "PA* .Xtt. Ytt. "LB 
6340  ENO  IF 

6343  IF  Sye-I  THEN  PRINT  "SH*“ 

6346  IF  Sye-2  THEN  PRINT  "SM+" 

6343  IF  Sy»-3  THEN  PRINT  "SMo" 

6352  Klab-1 
6355  ENO  IF 
6358  K/-1 

6361  CALL  Sywb(Xlab.Ylab.Sy»,Icl ,K; ) 

6364  PRINT  "SR  1.1.S:SM" 

6367  IF  Sye<4  THEN  PRINT  "PR  2,0' 

6370  BEEP 

6373  INPUT  "ENTER  BOP*. Bop 
6376  IF  Ifntt-1  THEN 
6379  BEEP 

6382  INPUT  "ENTER  TUBE  TYPE*. T type* 

6385  ENO  IF 

6388  IF  8op<J0  THEN  PRINT  "PR  2.0:LB* :3op,* "" 

6391  IF  Bop>9  THEN  PRINT  "PR  .5 , 0 JLB" ;8op ;  ** 

6394  Ihf-0 

6397  IF  Q! >02  THEN  Ihf-t 

6400  IF  Ihf-0  THEN  PRINT  "PR  4.0:LBInc* 

6403  IF  Ihf-I  THEN  PRINT  "PR  4,0:LB0ec" 

6406  IF  Ifntt-0  THEN 

6409  PRINT  "PR  2.0;LB':0_f ileS;** 

6412  ELSE 

64 IS  PRINT  "PR  2.0:LB";Ttyp«S;** 

6418  END  IF 

6421  PRINT  "SP0-.SP1  :SR  1.5.2" 

6424  Ylab-Ylab-5 
6427  ENO  IF 
6430  BEEP 

6433  ASSIGN  $Ftle  TO  * 

6436  Xll-XU/1.2 
6439  Xui-Xul*1 .2 
6442!  GOTO  8040 
6445  END  IF 
6448  PRINT  "PU  SM" 

645 t  BEEP 

64S4  INPUT  "WANT  TO  PLOT  A  POLYNOMIAL  <1 -Y, 0-N) 7* ,Go 

6457  IF  Go  on-1  THEN 

6460  BEEP 

6463  PRINTER  IS  1 

6466  PRINT  USING  "4X. “Select  line  type:""" 

6469  PRINT  USING  "6X.’"0  Solid  line"'" 

6472  PRINT  USING  *6X,*"1  Dashed""" 

6475  PRINT  USING  *6X.*"2..,5  Longer  line  -  dash""" 

6478  INPUT  Ipn 

6481  PRINTER  IS  70S 
6484  BEEP 

6487  INPUT  "SELECT  (O-LIN. 1-LOGtX. Y1 ) ",  Hog 

6490  Iprnt-1 

6493  CALL  Poly 

6496  IF  Ifn-0  THEN 

6499  N jo  in- 1 

6502  BEEP 

6505  INPUT  "ENTER  NUMBER  OF  FILES  TO  JOINT" ,N join 
6508  END  IF 

65 1 1  I  join-0 

6514  IF  Ifn<Njoin  THEN  Ijom-I 


65!  7 
6520 
6523 
6526 
6523 
6532 
6535 
6538 
6541 
6544 
6547 
6550 
6553 
6556 
6559 
6562 
6565 
6568 
6571 
6574 
6577 
6580 
6583 
6586 
6589 
6592 
6595 
6598 
6601 
6604 
6607 
6610 
6613 
6616 
6619 
6622 
6625 
6628 
6631 
6634 
6637 
6640 
6643 
6646 
6649 
6652 
6655 
6658 
6661 
6664 
6667 
6670 
6673 
6676 
6679 
6682 
6685 
6688 
6691 
6694 


IF  Ifn-0  OR  Ijom-1  THEN 

FOR  Ij-0  TO  3 

8s ( I  j )  -8s(  I; ) ♦B(  I j ) 

NEXT  I; 

If n«If n*1 
END  IF 

IF  Njoin-Ifn  THEN 
FOR  Ij-0  TO  3 
B< Ij)-8s( I;)/Njoin 
BstljJ-O 
NEXT  I; 

Ifn-0 

ELSE 

GOTO  6067 
END  IF 
BEEP 

INPUT  ‘ENTER  Y  LOWER  AND  UPPER  LIMITS' , YI1, Yui 
FOR  Xx-0  TO  Cx  STEP  Cx/200 
Xa-X«in*1 0  *Xx 

IF  Xa<Xll  OR  Xa>Xul  THEN  6655 

Icn-Icn+1 

Pu-0 

IF  Ipn-1  THEN  Idf-Icn  NOD  2 

IF  Ipn-2  THEN  Idf-Icn  MOD  4 

IF  Ipn-3  THEN  Idf-Icn  MOD  8 

IF  Ipn-4  THEN  Idf-Icn  MOD  16 

IF  Ipn-5  THEN  Idf-Icn  MOD  28 

IF  Idf-1  THEN  Pu-1 

IF  Opo-O  THEN  Ya-FNPoly(Xa) 

IF  Opo-2  AND  I log-0  THEN  Ya-Xa/FNPoly (Xa) 

IF  Opo-2  ANO  I log- 1  THEN  Ya-FNPoly(Xa) 

IF  0po-1  THEN  Ya-FNPoly(Xa) 

IF  Ya<  Yam  THEN  6655 
IF  Ya< Y 1 1  OR  Ya>Yul  THEN  6655 
IF  Iun-1  THEN 
IF  0po<2  THEN  Xa-Xa*1  .8 
IF  0po>0  THEN  Ya-Ya*. 1 761 
IF  Opo-O  THEN  Ya-Ya*. 31 7 
IF  Opo-2  THEN  Xa-Xa*.317 
END  IF 

Y-LGTf  Ya/Ymin)*Sfy 

X-LGT(Xa/Xmin)*Sf x 

IF  Y<0  THEN  Y-0 

IF  Y> 1 00  THEN  GOTO  6655 

IF  Pu-0  THEN  PRINT  "PA" , X . Y . *PD" 

IF  Pu-1  THEN  PRINT  "PA" ,X, Y . "PU" 

NEXT  Xx 
PRINT  "PU" 

GOTO  6067 
END  IF 
BEEP 

INPUT  "WANT  TO  PLOT  REILLY'S  SMOOTH-TUBE  DATA?  ( 1-Y.0-N) " , Irly 

IF  Qpo-0  OR  Qpo-1  THEN 

XI 1-3 

Xui-20 

END  IF 

IF  Opo-2  THEN 
X 1 1  - 10000 
Xui-100000 
END  IF 


'Cs.'C 


SS97  IF  Irly-1  THEN 
6700  Til *20 
6703  Yul-70 
6706  BEEP 

6709  INPUT  “ENTER  LOWER  ANO  UPPER  Y-LIMITS  FOR  PLOTTING" . Y 1 1 , Yul 

6712  FOR  Xx-0  TO  Cx  STEP  Cx/200 

6715  Xa-X*nv»1 0  ‘Xx 

6718  IF  Xa<Xl I  OR  Xa>Xul  THEN  6757 

6721  Xl-LOG(Xa) 

6724  IF  Opo-O  THEN  YL— 2 .S402837E- 1 +X 1*< 4 . 972015 1 -XI*2 .5 134787E- 1 ) 
6727  IF  Opo-1  THEN  YI— 2 .S402837E- 1 +X 1*< 3 .9720 1 5 1 -X 1*2 .5 1 34787E-1 ) 
6730  IF  Opo-2  THEN  Yl— 3.7073801E-t+Xl*<8.72S9190E-1-Xl*6.8826842E-3> 
6733  Ya-EXP(YI) 

6736  Y-LGT(Ya/Ymn)*Sfy 
6739  X-LGT  ( Xa/Xmn)  *Sf  x 
6742  Ipu-0 

6745  IF  Y< Yil  THEN  Ipu-1 

6748  IF  Y>Yul  THEN  GOTO  6757 

6751  IF  Ipu-0  THEN  PRINT  “PA“  .X.Y, “PD“ 

6754  IF  Ipu-1  THEN  PRINT  “PA*  ,X.Y,“P1J* 

6757  NEXT  Xx 
6760  PRINT  “PU“ 

6763  END  IF 
S766  BEEP 

6769  INPUT  “WANT  TO  PLOT  ROHSENOW  CORRELATION?  ( 1 -Y.O-N) *  .Irohs 

6772  IF  Irohs- 1  THEN 

6775  Yll-15 

6778  Yui-80 

6781  BEEP 

6784  INPUT  “ENTER  Tsat  (Deg  C)*.Tsat 
6787  Cat -.0040 
6790  BEEP 

6793  INPUT  "ENTER  Csf  (DEF-0.004) ‘ .Csf 
6736  Tf-Taat+2 

6799  FOR  Xx-0  TO  Cx  STEP  Cx/200 
6802  Xa-Xgun*10‘Xx 
6805  IF  Xa<Xl 1  OR  Xa>Xul  THEN  6889 
6808  Xl-LOG(Xa) 

6811  IF  0po<2  THEN  Tf-Tsat*Xa/2 
6814  Rho-FNRho( Tf 1 
6817  K-FNK(Tf) 

6820  Mu-FNMu(Tf) 

6823  Cp-FNCp(Tf) 

6826  Hfg-FNHf g( Tsat ) 

6329  Ni-Hu/Rho 
6832  Pr-Cp*Mu/K 

6835  0»ega-Cst*Hfg/Cp*( ( .014/(9.31 *Rho) ) * .S/(Mu*Hfg> ) * ( 1 ,/3) *Pr  ‘1.7 

6838  IF  Opo-O  THEN  Ya* ( Xa/Omegal ‘ 3 

6841  IF  Opo-1  THEN  Ya-( Xa/Owega) '3/Xa 

6844  IF  Opo-2  THEN  Ya-Xa ‘ ( 2./3J /Omega 

6847  IF  Opo-2  THEN 

6850  Tf c-Tsat+Xa/ Ya* .5 

6853  IF  ABS< Tf-T fc) > .01  THEN 

6856  Tf-(Tf*Tfc)*.5 

6859  GOTO  6814 

6862  END  IF 

6865  ENO  IF 

6868  Y-LGT (Ya/Ysin)*Sfy 

6871  X-LGT(Xa/Xrain)*Sf x 

6874  Ipu-0 
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7060 

7063 

7066 

7063 

7072 

7075 

7078 

7 08 1 

7084 

7087 

7090 

7093 

7096 

7099 

7102 

7105 

7108 

7111 

7114 

7117 

7120 

7123 

7126 

7129 

7132 

7135 

7138 

7141 

7144 

7147 

7150 

71  S3 

7156 

7159 

7162 

7165 

7168 

7171 

7174 

7177 

E-10) 

7180 

7183 

7186 

7189 

7192 

7195 

7198 

7201 

7204 

7207 

7210 

7213 

7216 

7219 

7222 

7225 

7228 

7231 

7224 


PRINT 

PRINT  USING  "10X.""  Fms  Tin  Tev 

IF  K-0  THEN 
PRINTER  IS  701 
PRINT 

PRINT  ’  Month,  data  and  time:  *:Date$ 

IF  Itt-0  THEN  PRINT  USING  “1  OX.* “Tube  Typa: 

IF  Itt-1  THEM  PRINT  USING  *10X.**Tuba  Type: 

IF  Itt-2  THEN  PRINT  USING  *10X.**Tuba  Type: 

PRINT 

PRINT  USING  - 1 0X. * “  Fms  Tin  Tev 


Tdrop" * ’ 


Wieiand  Smooth"” 
High  Flux"** 
TurDo-8"*" 


PRINT  USING  *10X.**  Fms  Tin  Tev  Vw42  Tdrop"** 

PRINTER  IS  1 
K-t 

EN0  IF 
BEEP 

INPUT  'ENTER  FLOWMETER  READING" .Fms 
OUTPUT  709: “AR  AFO  AL4  VR1* 

FOR  L-0  TO  4 
OUTPUT  709: *AS  SA* 

IF  L>0  AND  L<4  THEN  7141 
S-0 

FOR  1-0  TO  9 
ENTER  709 :£ 

S-S+E 
NEXT  I 

IF  L-0  THEN  Emt (0) -ABS (S/10) 

IF  L-4  THEN  Emt ( 1 )-ABS(S/1 0 ) 

NEXT  L 

OUTPUT  709: “AR  AF20  AL20  VR1 * 

OUTPUT  709: *AS  SA* 

Etp-0 

FOR  1-0  TO  9 
ENTER  709:Et 
Etp-Etp+Et 
NEXT  I 
Etp-Etp/1 0 
T in-FNTvsv( Emt ( 1  ) ) 

Tev-FNTvsv<  Emt ( 0  )  ) 

Grad-37. 3853-t.  J04388*T in 

Mdot-3.9657E-3+Fms*(3.619S5E-3-Fms*<8.32006E-6-Fms*( 1 .23688E-7-Fm*r*4.31897 
)  ) 

Uw-Mdot/ ( 1 00Q*PI*Di ’2) *4 
Tdrop-£tp*1 .£+6/ ( 1 0*Grad ) 

PRINT  USING  *1  OX ,3( DD . DD . 4X) , IX . Z.DD . 4X ,MZ. 4D" ;Fms. Tin, Tev. Um' 2, Tdrop 
BEEP 

INPUT  ’WANT  TO  ACCEPT  THIS  DATA  SET?  ( 1 -T ,0-N)  *  ,01c 
J-J*1 

IF  Ok-O  THEN 
J-J-1 
GOTO  7102 
ELSE 

OUTPUT  df  i le:Fms.Erat (*) ,Etp 
PRINTER  IS  701 

PRINT  USING  *  1 0X,3( DD.DD.4X ) , 1 X ,Z.DD , 4X ,MZ . 4D* ;Fms. Tin. Tev, Vw ‘2. Tdrop 

PRINTER  IS  1 

BEEP 

INPUT  “WILL  THERE  BE  ANOTHER  DATA  SET?  ( I -Y.O-N) * ,Go_on 
IF  Go  on-1  THEN  7102 
END  IF 

ASSIGN  4F i 1 a  TO  * 


7237  PRINTER  IS  701 

7240  PRINT  .  .  .  ,  _ 

7243  PRINT  USING  *  1 0X . ’ “NOTE :  **, ZZ .**  data  3«ts  are  storea  in  file  .15A  ,J.F 
i  lei 

7246  PRINTER  IS  1 
7249  SUBEND 
72S2  SUB  Uoort 
72SS  PRINTER  IS  1 
72S8  8EEP 

7261  INPUT  ‘Enter  Uo  File  Naae’.FileS 
7264  BEEP 

7 267  INPUT  "Nuaber  of  Data  Runs’.Nrun 

7270  INPUT  ‘Do  You  Want  a  Plot  Fi ie?<  1  -Y . 0-N ) * , Ip  lot 

7273  BEEP 

7276  IF  Iplot-I  THEN 

7279  INPUT  ‘Give  Plot  File  Naee* ,P_f i leS 

7282  CREATE  BOAT  P_fileS.4 

728S  ASSIGN  SPiot  TO  P_f ileS 

7288  END  IF 

7291  PRINTER  IS  701 

7294  PRINT 

7297  PRINT 

7300  PRINT  USING  *10X.**  Water  Yel  Uo’** 

7303  ASSIGN  SFiie  TO  FileS 

7306  IF  Iplot-1  THEN  ASSIGN  afilel  TO  PJ’iieS 

7309  FOR  1-1  TO  Nrun 

7312  ENTER  9File:Yu.Uo 

7315  IF  Ip  lot—  1  THEN  OUTPUT  «File1:Vu.Uo 

7318  PRINT  USING  *1SX.D.D0.6X.rtZ.3DE* iVy.Uo 

7321  NEXT  I 

7324  ASSIGN  afile  TO  ♦ 

7327  ASSIGN  *File1  TO  *  ....  „ 

7330  PRINT  USING  *10X.*‘N0TE:  **.ZZ.**  data  sets  are  stored  in  file  ‘‘.I5A  ;Nru 
n.FileS 


7333  IF  Ip  lot—  1  THEN  ,  ,  .. 

7336  PRINT  USING  "10X,**N0TE:  **.22.**  X-T  Pairs  are  stored  in  file  ‘*,I5A  ;Nru 

n,P_f i leS 

7339  END  IF 

7342  PRINTER  IS  1 

7345  SUBEND 


7348 

7351 

7354 

7357- 

7360 

7363 

7366 

7363 

7372 

7375 

7378 

7381 

7384 

7387 

7390 

7393 

7396 

7399 

7402 

7405 


SUB  Select 
CON  /Ido/  Idp 
BEEP 

PRINTER  IS  1 


PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 


*4X. "’Select  option:1 


PRINT  USING 
PRINT  USING  *6X 
PRINT  USING  "SX . 

PRINT  USING 
PRINT  USING 
INPUT  Idp 

IF  Idp-0  THEN  CALL  flam 
IF  Idp-1  THEN  CALL  Plot 


Taking  data  or  re-process i ng 
Plotting  data  on  Log-Log_”* 
Plotting  data  on  Linear1”* 
Make  cross-plot  coe+'ft  file** 
Re-circuiate  water**’ 

Purge’** 

T-Orop  correction’** 

Print  Uo  File*"" 

8  Modify  X-Y  file*** 

9  Move'** 

•sx.**10  Comb*** 

*SX, **T 1  Fixup’** 


’SX.** 
"SX.” 
*6X,  "* 
*6X . " * 
*6X. *  * 
’SX.** 

*sx.** 

*sx.** 
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old  data’** 


7408  IF  Idp-2  THEN  CALL  PI  in 
741 1  IF  Idp-3  THEN  CALL  Coet 
7414  IF  Idp-4  THEN  CALL  Main 
7417  IF  Idp-5  THEN  CALL  Purg 
7420  IF  Idp-S  THEN  CALL  Tdcn 
7423  IF  Idp-7  THEN  CALL  Uoprt 
7426  IF  Idp-8  THEN  CALL  Xymod 
7429  IF  Idp-9  THEN  CALL  Move 
7432  IF  Idp-10  THEN  CALL  Comb 
743S  IF  Idp-11  THEN  CALL  Fixup 
7438  SUBENO 
7441  SUB  Xymod 
7444  PRINTER  IS  t 
7447  BEEP 

7450  INPUT  "ENTER  FILE  NAME". File* 

7453  ASSIGN  SFilel  TO  FileS 
7456  BEEP 

7459  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS" ,Np 
7462  BEEP 

7465  INPUT  "ENTER  NEW  FILE  NAME’,File2S 
7468  CREATE  BDAT  Flie2S.S 
7471  ASSIGN  $File2  TO  FUe2S 
7474  BEEP 

7477  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS  TO  BE  DELETED* .Ndei 
7480  IF  NdeI-0  THEN  7492 
7483  FOR  1-1  TO  Ndei 
7486  BEEP 

7489  INPUT  "ENTER  DATA  SET  NUMBER  TO  BE  DELETED ", Nd ( I ) 

7492  NEXT  I 

7495  FOR  J-1  TO  Np 

7498  ENTER  *FUe1  :X.Y 

7501  FOR  1-1  TO  Ndei 

7504  IF  Nd(I)-J  THEN  7S16 

7507  NEXT  I 

7510  OUTPUT  3File2:X.Y 

7513  PRINT  J.X.  Y 

7516  NEXT  J 

7519  PRINTER  IS  701 

7522  ASSIGN  AFilel  TO  * 

7525  ASSIGN  .4Ftie2  TO  * 

7528  SUBENO 
7531  SUB  Move 
7534!  FILE  NAME:  MOVE 
7537! 

7540  DIM  Bop(SG) ,A<66) .8(66) .CC66) ,D<66) ,E(66) .F<66).G<66) ,H(66) ,J<66) .K(66)  ,L( 
66) . M(66) 

7543  DIM  ToidSf  66) [ 14] ,N(66) , Vr (66 ) , Ir ( 66) 

7546  BEEP 

7549  INPUT  "OLD  FILE  TO  MOVE" ,D2_f i leS 

7552  ASSIGN  4FUe2  TO  D2_fileS 

7555  ENTER  4Fi  le2:Nrun. Dates, Ldtcl ,Ldtc2.Itt 

7558  FOR  1-1  TO  Nrun 

7561  ENTER  »F!  ,e2:Boo( I) .ToidSI I) 

7564  ENTER  JFi*e2:rt< I) ,B( I) ,C( I) ,D< I) ,E( I) ,F< I) ,G( I> ,H< I) , J( I) ,K( I) ,L< I)  ,M<  I>  ,N 
(I) 

7567  ENTER  »Fi le2 : Vr ( I ) , Ir ( I) 

7570  NEXT  I 

7573  ASSIGN  »Fiie2  TO  ♦ 

7576  BEEP 

7579  INPUT  "SHIFT  DISK  AND  HIT  CONTINUE* .  01c 
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j 


7532  3£cp 

7535  INPUT  -INPUT  3DAT  SIZE'. Sue 
7538  CREATE  3DAT  D2_f i leS .Si te 
7591  ASSIGN  if  i  le I  TO  D2_fileS 
7594  OUTPUT  *Fi lei ;Nrun. Date*. Ldtcl .Ldtc2. Itt 
7597  FOR  1*1  TO  Nrun 
7600  OUTPUT  »File1;8op(I).ToldS(I) 

7603  OUTPUT  if  i  lei  :A(  I )  .3(  I)  ,C(  I )  .D(  I )  .E(  I)  .F(I )  ,G(  I)  ,H(  I) .  J(I)  .KlI )  ,L(  I)  .mi ) 
N(  I) 

7606  OUTPUT  3Fi iel :Vr ( I) , Ir( I) 

7609  NEXT  I 

7612  ASSIGN  if lie!  TO  ♦ 

7615!  RENAME  'TEST'  TO  D2  file* 

7613  BEEP  2000.. 2 
7621  BEEP  4000.. 2 
7624  BEEP  4000.. 2 

7627  PRINT  'DATA  FILE  MOVED' 

7630  SUBENO 
7633  SUB  Como 
7636!  FILE  NAME:  COMB 
7639! 

7642  DIM  Emt (12) 

7645  BEEP 

7648  INPUT  “OLD  FILE  TO  FIXUP*. D2  file* 

7651  ASSIGN  9file2  TO  D2  file* 

7654  D1  fileS-'TEST* 

7657  CREATE  BDAT  01  file*. 30 

7660  ASSIGN  if  i  lei  To  Q1_fileS 

7663  ENTER  Afile2:Nrun. Date*. Ldtcl .Ldtc2. Itt 

7666  IF  K-Q  THEN  OUTPUT  ifi lei ;Nrun, Date*. Ldtcl ,Ldtc2. Itt 

7669  FOR  1*1  TO  Nrun 

7672  ENTER  ifile2:3oo.ToldS.Emf (*) ,Vr .Ir 
7675  OUTPUT  *f; lei :8op . Told* .Emf (*) ,Vr . Ir 
7678  NEXT  I 

7681  ASSIGN  if  l Ie2  TO  * 

7684!  RENAME  'TEST*  TO  Q2_file* 

7687  BEEP  2000.. 2 
7690  BEEP  4000.. 2 
7693  BEEP  4000.. 2 
7696  BEEP 

7699  INPUT  -UANT  TO  ADD  ANOTHER  FILE  ( 1 -Y . 0-N ) ?*  . Ok a 
7702  IF  Oka-1  THEN 
7705  K-1 

7708  BEEP 

7711  INPUT  -GIVE  NEW  FILE  NAME"  ,Nf i leS 
7714  ASSIGN  4Ftle2  TO  Nf^eS 
7717  GOTO  7663 
7720  END  IF 

7723  ASSIGN  ifile2  TO  ♦ 

7726  SUBEND 

7729  SUB  Fixup 

7732!  FILE:  FIXUP 

7735!  DATE:  Feoruary  18.  1986 

7738! 

7741  DIM  Emf (12) 

774a  BEEP 

7747  INPUT  "OLD  FILE  TO  FIXUP' . D2_f i ie* 

7750  ASSIGN  if  t le2  TO  Q2_fileS 
7753  D!  f i leS-'TEST’ 

7756  CREATE  BDAT  D1_f  i  l““i  .20 
7759  ASSIGN  ifilel  TO  OI_fileS 
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7762  ENTER  «F  i  ia2 :Nrun . Dates  .Ldtc l  . i_dtc2 .  It  t 
7755  Nemo*  2  7 

7763  OUTPUT  4F  i  le  I  1  Nr'jno  .Da  teS  .Ldtc  1  .Ldtc2 , 1 1 1 
7771  FOR  1-1  TO  Nrun 

7774  ENTER  »F  t  le2  :Sod.  ToldS  .E.nt  (  *)  .Vr  .  Ir 
7777  IF  1-1  THEN  7783 

7780  OUTPUT  sF i le I :3oe . ToldS.Emt ( *) . Vr . Ir 
7783  NEXT  I 

7786  ASSIGN  »Fiie2  TO  ♦ 

7739  ASSIGN  SFilal  TO  * 

7792!  RENAME  'TEST'  TO  D2  file* 

7795  BEEP  2000.. 2 
7798  3EEP  4000.. 2 
7801  BEEP  4000,. 2 
7804  SUBENO 


i  m '  if1  i  tf  u  f  up  up  i  ■  j 


APPENDIX  C 

EXAMPLES  OF  REPRESENTATIVE  DATA  RUNS 
The  two  following  data  runs  (TXE222  and  GTB185)  are 
representative  samples  of  the  data  taken  over  the  course  of 


i 

1 

J 


the  Investigation.  The  first  data  run  for  the  Thermoexce 1 -E 
surface  was  for  a  decreasing  heat-flux  condition  and  an  oil 
concentration  of  three  percent.  The  second  data  run  was  for 
the  GEUA-T  surface  was  also  for  a  decreasing  heat-flux 
condition  but  with  an  oil  concentration  of  one  percent. 


Month,  da  to  and  tiM  :  OS: 06: 11  :*1 :20 

NOTE:  Proqras  nas#  :  DRP4 
Disk  nusOsr  •  13 
Old  fils  nas#:  TXE222 


This  data 

sst  taksn  on  :  02: 18: 

13:42:58 

Tubs  Nuaosr:  10 

Ooto  Sot 

Nusosr 

-  i  Suite  an  t  - 

3.0 

02: 

18:13:48:44 

TC  No: 

1 

2  3  4 

5 

6 

7  8 

T  mmo  : 

3.71 

9.45  9.00  8.48 

9.70 

8.12 

10.21  8.59 

Tm 

T1  iad 

Tliqd2  Tvoor  Psot 

Tsuoo 

a. 7. 

2.23 

2.17  2.06  -1.27 

-18.9 

Thotoo 

Htuds 

Ode 

S.S4I 

l.4SSE‘04  9.SI9E‘04 

Doto  Sot 

Nusbsr 

•  2  Bulk  Oil  Z  - 

3.0 

02: 

i8:i3:49:S7 

TC  No: 

1 

2  3  4 

5 

6 

7  8 

Tmo  : 

8. SI 

9.01  8.42  8.00 

9.58 

3.29 

10.33  8.51 

Two 

Tl  iqa 

Tltad2  Tvoor  Psot 

Tsuoo 

8.55 

2.24 

2.22  2.09  -1.73 

-18.3 

Thotoo 

Htuos 

Qdo 

S.3I3 

1.S08E--04  9.53IE‘04 

Doto  Sot 

Nuaosr 

■  3  Built  an  z  • 

3.0 

02:18: 13:50:51 

TC  No: 

1 

2  3  4 

5 

6 

7  3 

Tooo  : 

9.68 

8.38  8.33  7.89 

9.84 

8.07 

10.2 Z  8.64 

Two 

TUdd 

Tliad2  Tvoor  Psot 

Tsuoo 

8.53 

2.26 

2.19  2.07  -1.74 

-18.3 

Thotoo 

Htubs 

Qdo 

5.305 

1  .SUE- 

►04  9.52SE*04 

Ooto  Sot 

Nuaosr 

-  4  Built  an  7.  * 

3.0 

02: 

18:14:02:37 

TC  No: 

1 

2  3  4 

5 

6 

7  8 

Toon  : 

6.93 

7.05  5.71  6.48 

7.23 

5.24 

7.63  6.56 

Two 

Tl  idd 

Tltod2  Tvoor  Psot 

Tsuoo 

5.56 

2.2S 

2.21  2.14  -1.73 

-18.  S 

ThotoO 

Htuds 

Qdo 

4.433 

1.41 4£.04  S.269E‘04 

Ooto  Sot 

Nusosr 

•  5  Bulk  an  x  - 

3.0 

02: 

18:14:03:13 

TC  No: 

1 

2  3  4 

5 

6 

7  8 

T.oo  : 

6.32 

7.18  5.30  S. 56 

7.35 

6.25 

7.55  5.52 

Two 

Tl  idd 

Tl:od2  Tvoor  Psot 

Tsuoo 

5.71 

2.28 

2.25  2.18  -1.59 

-18.4 

ThotoO 

Htuos 

Qdo 

4.444 

l.408E‘04  S.2S8£‘04 

Doto  Sot 

Huadsr 

-  5  Bulk  Oil  1  • 

3.0 

02: 

18:14: 13:S6 

TC  No: 

1 

2  3  4 

5 

5 

7  8 

T.oo  : 

s.so 

5.55  5.41  5.26 

5.58 

4.30 

5.34  5.06 

Two 

Tl  iqd 

Tliod2  Tv  tor  Psot 

Tsuoo 

5.28 

2.32 

2.28  2.35  -1.58 

-16.3 

Thotoo 

Htuos 

Qdo 

2.381 

I .3046 

‘04  3.388£‘Q4 

7 


Data  Set  Nuaoar  ■  7  Bulk  Oil  7.  •  3.J 

TC  No:  !  2  3  a  5 

Taae  :  S.4S  S.SO  5.37  5.31  S.6I 

Tua  Tliqd  Tliqd2  Tvaer  Peat  Tium 

5.26  2.23  2.26  2.35  -1.68  -16.2 

Thatao  Htuoa  Ode 

2.387  1 .2996*04  3.8816*04 

Data  Sat  Nuaoar  *  3  Sulk  Oil  Z  •  3.0 

TC  No:  l  2  3  a  S 

Taae  :  a.2S  4.33  4.22  4.!6  4.49 

Tua  Tliqd  Tliqd2  Tvaer  Plat  Tsuae 

4.20  2.24  2.16  2.64  -1.77  -15.8 

Thatao  Htuoa  Ode 

1.898  1 . 1 766*04  2.2496*04 


Oata  Sat  Nuaoer  -  9  Bulk  Oil  7.  -  3.0 

TC  No:  I  2  3  4  5 

Taae  :  4.33  4. an  4.26  4.23  4.52 

Tua  Tiled  Tliqd2  Tvaer  Peat  Tauae 

4.25  2.28  2.23  2.71  -1.71  -15.8 

Thatao  Htuoa  Ode 
2.004  1.1636*04  2.3446*04 

Oata  Sat  NuaOar  *  10  Bulk  Oil  Z  *  3.0 

TC  No:  1  2  3  4  5 

Taao  :  3.68  3.68  3. SB  3.S4  3.30 

Tua  Tiled  Tliqd2  Tvaer  Peat  Tauae 

3.S9  2.33  2.27  2.35  -1.66  -16.4 

Thatao  Htuoa  Ode 

1.293  1.1376*04  1.4696*04 


Oata  Sat  NuaOar  ■  11  Bulk  011  2  *  3.0 

TC  No:  I  2  3  4  S 

Taao  :  3.S6  3.61  3.SS  3.51  3.75 

Tua  Tliqd  Tliqd2  Tvaer  Paat  Tauae 

3.53  2.28  2.21  2.34  -1.72  -16.4 

Thetab  Htube  Ode 

I. 285  I . 1 42E*04  1.4676*04 


02:18:14:14:29 
6  7  8 

4.88  5.82  5.02 


02: tS: 1 4:20:33 
6  7  9 

4.07  4.58  4.06 


02:18:14:21:05 
6  7  8 

4.10  4.65  4.11 


02:18.-14:26:17 

6  7  a 

3. SI  3.3S  3.47 


02:18:14:26:48 
6  7  8 

3.44  3.80  3.39 


Oata  Sat  Nuabtr  •  12  Bulk  Gil  2  • 

3.0  02:18:14:33 

:  1 1 

TC  No: 

1 

2  3  4 

5  6  7 

3 

Taae  : 

3.06 

3.07  3.02  3.00 

3.12  2.38  3.14 

2.30 

Tua 

Tliqd 

Tliqd2  Tvaer  Peat 

Teuae 

3.01 

2.43 

2.11  3.43  -1.69 

-17.1 

Thatab 

Htuoa 

Ode 

.739 

1.1206*04  8.2836*03 

Oata  Sat  NuaOar  -  13  Bulk  Oil  Z  - 

3.0  02:18:14:33 

:43 

TC  No: 

1 

2  3  4 

5  6  7 

9 

Taae  : 

3.00 

3.01  2.94  2.92 

3.05  2.90  3.08 

2.81 

Tua 

Tliqd 

Tliqd2  Tvaer  Peat 

Teuae 

2.34 

2.36 

2.05  3.41  -1.77 

-17.2 

Thatao 

Htuoa 

Ode 

■  .738 

1.1366*04  8.3866*03 

Data  Sat  Nua6«r  •  14  Sulk  Oil  7.  - 

3.0  02:18:14:40 

:43 

TC  No: 

1 

2  3  4 

5  6  7 

8 

Taao  : 

2.70 

2.74  2.73  2.71 

2.30  2.71  2.79 

2.61 

Tua 

Tl  iod 

Tliqd2  Tvaer  Peat 

Teuae 

2.71 

2.4! 

2.03  3.79  -1.75 

-17.3 

Thatao 

Htuba 

Qdp 

.488 

1 .07SE* 

04  5.2576*03 

Oata  Sat 

Nuabar 

-  IS  Bulk  3‘1  7.  -  3 

.0  02:18:14:41:' 

19 

TC  No: 

1 

2.66 

2  3  4 

5  6  7 

3 

Taao  : 

2.70  2.65  2.67  2 

.74  2.66  2.71 

2.54 

Tua  Tliqd 

Tliqd2  Tvaer  Peat 

Teuae 

2.65 

2.39 

2.00  3.78  -1.77 

-18.0 

Thatao 

Htube 

Ode 

.458 

1 . 1 456* 

04  5.2446*03 
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Month.  Jot*  and  tiao  :US:06:  It  :a7:  44 


MOTE:  Proqraa  naao  :  D8P4 
Dlak  nuaeor  •  <0 
Old  fiio  naao:  STBISS 
Thia  data  aot  takon  on  :  It :2S: 11:35:33 
Tub*  Nuaeor:  9 


Data  Sat  Nuaeor  •  I  Sulk  0:1  Z  -  1.3 

TC  No:  t  2  3  4  S 

Tooa  :  9.27  8.61  7.83  7.37  7.83 

Too  TUod  Tliad2  Twaer  Put 


1 1 : 25: IS: 46: 12 
6  7 

7.82  9.8S  7 


Thotae  HtuUo  Ode 

5.728  t.S82£*04  9.062E*04 


Data  Sot  Nuaeor 


3  Bulk  Oil  Z 


II  :2S: IS:S4:0S 


Tea 

Tllod  TUod2 

Tv*pr  Put 

Tauao 

6.66 

2.2!  2.22 

2.77  -1.75 

-17.0 

Thotab 

Htubm  Odp 

4.448 

t.348€*04  5.39SE-Q4 

Data  Sot 

Nuaeor  -  4 

Bulk  Oil  7.  - 

1.0  1 t :2S: !5:S4:S5 

Tliqd2  Tvaor  Paat 


Tliod2  Tvaor  Paat 


Thotae 

Htufcp 

Ode 

3.S83 

I . Q36E-04  3.7I2E-04 

Oata  Sot 

Nuobor 

•  6  Bulk  Oil  Z  - 

u.a 

ll  :2S: 16:04:22 

TC  No: 

1 

2  3  4 

s 

6  7  1 

Taao  : 

6.11 

6.23  S.84  5.84 

5.80 

5.77  6.42  5 

Tua 

TUod 

Tliqd2  Tvaor  Paat 

Tauao 

5.85 

2.27 

2.26  3.22  -1.69 

-17.4 

Thotab 

Htubo 

Ode 

3.588 

1. 0346*04  3.7T2E-04 

Oata  Sot  Nuaeor 

-  7  Bulk  Oil  X  - 

1.0 

1 1 : 25: 16: 10:  IS 

TC  No: 

t 

2  3  4 

5 

6  7  1 

Toaa  : 

5.41 

S.S2  5.19  5.32 

5.25 

5.16  5.50  5 

Tua 

T1  iqd 

TUqd2  Tvaor  Poat 

Tauao 

5.24 

2.24 

2.20  3.61  -1.74 

-17.7 

Thotab 

Htuto* 

Ode 

3.020 

7.320E*03  2.21 1E*04 

Data  Sot 

Nuab«r 

-  8  Bulk  OH  Z  • 

1 .0 

1 1:25:16:12:1 1 

TC  No: 

1 

2  3  4 

5 

6  7  8 

Taoe  : 

5.2S 

S.32  4.38  5.13 

5.06 

4.95  5.29  4. 

Tua 

Ti  lod 

Tllqd2  Tvaor  Paat 

Tauao 

Thotab 

Htub« 

Ode 

2.663 

S.lueHKJ 

1 .362E*04 

Oata  Sot 
rr  . 

N«i«to«r  - 

10  Sulk  Oil  Z  - 

1.0 

11:25: 16: 18:04 

Oata  Sat  Nuaber  •  It  Bulk  Gil  X  *  1.0 

TC  No:  I  2  3  «  5 

Teeo  :  4.75  4.88  4.42  4.34  4.75 

Tua  Tliod  Tliqd2  Tvaor  Paat  Tauae 

4.SS  2.30  2.29  3. IS  -1.66  -IS. 3 

Thatao  Htuoa  Ode 

2.3S0  3.273£»03  7.72SE»03 


Oata  Sat 

Nu*o«r 

•  12 

Bulk  GU  X  - 

1.0 

TC  No: 

1 

2 

3  4 

5 

Taae  : 

4.76 

4.38 

4.42  4.34 

4.74 

Tea 

TUfld 

Tliod2 

Tvaor  Paat 

T«u«o 

4.66 

2.30 

2.29 

3.13  -1.66 

-16.0 

TKatad 

Htuo# 

Ode 

2.360 

3.264EHJ3  7. 

702E.03 

Oata  Sat 

Nu*o«r 

-  13 

Bulk  Oil  X  - 

1.0 

c 

12  3  4 

4.45  4.54  4.09  4. SO 

Tliod  Tliod2  Tvaor  Paat 

2.27  2.27  2. 35  -t .69 

Htuoa  Ode 
2.361E-03  4.944E-03 


TC  No: 

Taao  : 

Tua 

4.37 
The  tab 
2.094 

Oata  Sat  Nuacar  *  14  Bulk  Gil  2  ■ 

TC  No: 

"Tea 

4.39 
Thacao 
2.096 


12  3  4 

4.46  4.55  4.12  4.61 

Tliod  Tlioo2  Tvaor  Paat 

2.30  2.29  2.94  -1.66 

Htuba  Ode 

2.352E-03  4.931£*03 


4.39 

Tauae 

-14.9 


1.0 

5 

4.43 

Tauae 

-14.9 


1 1 :2S: >6:24:54 
6  7  3 

4.61  4.62  4.58 


11:25:16:25:38 
6  7  8 

4.61  4.60  4.58 


1 1 :2S: 16:31 :21 
6  7  9 

4.34  4.36  4.28 


1 1 :2S: 16:31 : 53 
6  7  8 

4.38  4.38  4.32 


Oata  Sat 

Nuaoar  *  15  Bulk  Oil  X  * 

1.0 

TC  No: 

1  2  3 

4 

S 

Taae  : 

a.  19  4.29  3.84 

4.38 

4.07 

Tua 

Tliod  Tliqd2  Tvaor 

Paat 

Tauae 

4.09 

2.23  2.23  2.83 

-1.73 

-14.8 

Thetao 

Htuoa  Ode 

1.863 

I.S82E.03  2.346E.03 

Data  Sat 

Nuaber  *  16  Bulk  Oil  X  - 

1.0 

TC  No: 

1  2  3 

4 

5 

Taae  : 

4.20  4.30  3.83 

4.36 

4.07 

Tua 

Tliqd  Tliqd2  Tvaor 

Paat 

Tauae 

4.09 

2.23  2.22  2.82 

-1.74 

-14.8 

Thatao 

Htuba  Qdo 

1.862 

1.S82E-03  2.946E*03 

Data  Sat  Nuaoar  *  17  Bulk  Oil  X  • 

1.0 

TC  No: 

1  2  3 

4 

5 

Taae  : 

3.97  4.04  3.67 

4.17 

3.77 

Tua 

Tliod  Tliqd2  Tvaor 

Psat 

Tauae 

3.85 

2.19  2.19  2.71 

-1.78 

-15.2 

The  tab 

Htuba  Qde 

1.660 

1.I46€*03  I.901E.03 

Oata  Sat 

Nuaber  -  18  Bulk  Oil  X  * 

1.0 

TC  No: 

1  2  3 

4 

5 

Taae  : 

3.97  4.04  3.68 

4.19 

3.78 

Tua 

Tliqd  Tliqd2  Tvaor 

Paat 

Tauae 

3. as 

2.19  2.19  2.71 

-1.77 

-15.2 

The tab 

Htuba  Qdo 

1 .660 

1 . 1 47E*03  1.904E-03 

Data  Sat 

Nuaber  •  19  Bulk  Oil  X  - 

1.0 

TC  No: 

1  2  3 

4 

5 

Taw  : 

3.69  3.73  3.50 

3.77 

3.40 

Tua 

Tliqd  Tliqd2  Tvaor 

Paat 

Tauae 

3.51 

2.19  2.19  2.58 

-1.78 

-15.4 

The  tab 

Htuba  Qdp 

1.326 

8.573E*02  1 . I37E*03 

Oata  Sat 

Nuaoar  '  20  Bulk  Oil  X  • 

1 .0 

TC  No: 

1  2  3 

4 

S 

Two  : 

3.77  3.81  3.57 

3.84 

3.46 

Tua 

Tliqd  Tliqd2  Tvaor 

Paat 

Tauae 

3.S8 

2.25  2.25  2.62 

-1 .71 

-15.7 

The  tab 

Htuba  Qde 

1.335 

9.S30E*02  1 . I52E»03 

l 1 :2S: IS:37:S3 
6  7  8 

4.06  4.00  3.99 


il:2S: <6:38:28 
6  7  8 

4.05  3.38  3.99 


H:2S:l6:4«:S6 
5  7  8 

3.81  3.67  3.74 


11:25:16:45:33 
6  7  8 

3.82  3.66  3.72 


11:25:16:51:19 

s  7  a 

3.46  3.25  3.32 


1 1 :25: 16:53: 16 
6  7  8 

3.54  3.32  3.38 
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APPENDIX  D 


UNCERTAINTY  ANALYSIS 


The  uncertainty  of  the  heat- transfer  coefficient  at  60 
kW/m2  and  8  kW/m2  for  data  runs  of  the  GEWA-T  tube  (GTB185 
with  one  percent  oil  concentration)  and  of  the  Thermoexcel- 
E  (TXE222  with  three  percent  oil  concentration)  are 
analyzed  below.  The  method  of  analysis  is  based  on  Kline- 
McClintock  [Ref.  20]  method  of  uncertainty  analysis. 

The  heat-transf er  coefficient  is: 

h  =  <qc)/(T«0  —  T| »  t )  (D.l) 

and 


-  T« , 


=  Tm  i  -CQc  In  (D,/D,  )]/(2  IT  k  L)  -  T, 


(D.  2) 


where 


h 

qc 

Tm  o 
T*  *  T 
T„  , 
Qc 
D, 

Dj 

k 

L 


heat- transfer  coefficient 

heat  flux  corrected  for  end  loss 

average  outer  tube  wall  temperature 

saturation  temperature 

average  inner  tube  wall  temperature 

heat  input  corrected  for  end  losses 

tube  internal  diameter 

tube  external  diameter 

tube  wall  thermal  conductivity 

length  of  tube's  heated  surface 


« 
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To  begin  the  analysis,  a  dummy  variable  is  assigned  to  the 
conduction  term  as  follows: 

DELT  =  C Qc  In  (D,/D,  )]/(  2  7Tk  L)  (D.3) 

In  accordance  with  Kline  and  McClintock.  the 

uncertainty  of  the  heat-transfer  coefficient  is: 


The  uncertainty  for  the  conduction  term  is  estimated  since 
the  logarithmic  error  is  negligible  when  compared  to  the 
other  error  of  the  analysis.  The  uncertainty  of  the 
conduction  term  is  as  follows: 


Sdelt 

DELT 


where  Qc  =  qc  7T  D?  L 

The  uncertainty  of  this  term  is  as  follows: 


5qc  = 

Qc 


ft)!  *  ft)1  *  (ft’ 


(D.  6) 


I  <  3 

(D.  7) 


Table  4  contains  a  listing  of  the  computed  terms  of 
equations  D. 1  through  D.7  for  selected  data  runs  for  the 
GEWA-T  and  the  Thermoexce 1 -E  tubes  at  specified  heat-flux 
setting  of  60  kW/m!  and  8  kW/m: .  The  algorithm  in  the  data- 
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TABLE  4 

.  UNCERTAINTY  ANALYSIS 


File  Name 
Heat  Flux 
Oi  1  * 

GTB185 
60  kW/m2  8 

1% 

kW/m2 

1% 

TXE222 
60  kW/m2 

3% 

8  kW/m2 

3% 

49a.. 

Qc 

0.015 

0.015 

0.015 

0.015 

4k. 

k 

0.03 

0.03 

0.03 

4l 

L 

0.0025 

0.0025 

0.0025 

0.0025 

DELT 

0.486 

0.064 

0.486 

0.064 

6  DELT 

DELT 

0.  069 

0.520 

0.069 

0.526 

T«  i  (  C ) 

6.845 

4.678 

6.904 

2.964 

Tm  t  C  C) 

0.49 

0.065 

0.  49 

0.065 

Tb.t  (  C) 

2.215 

2.295 

2.  265 

2.  205 

T.,t  <K) 

0.  1 

0.  1 

0.  1 

0.  1 

ILu. 

T„0  -T.  a  T 

0.  Ill 

0.028 

0.111 

0.088 

6  DELT 

Two  ~T  SAT 

O.  008 

0.014 

0.008 

0.046 

....  4.Ta..0...T _ 

Tuo  ~Ts  A  T 

0.  023 

0.042 

0.023 

0.  135 

s9c 

qc 

0.  02 

0.02 

0.02 

0.02 

5h 

h 

0.  115 

0.056 

0.  115 

0.  169 

h  ( W/m2  K  ) 


13480 


3273 


14140 


11200 


reduction  program  used  this  approach  and  assumed  that  all 
the  eight  wall  thermocouples  were  used  to  calculate  the 
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